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Slowly driven elastic interfaces, such as domain walls in dirty magnets, contact lines wetting a non- 
homogenous substrate, or cracks in brittle disordered material proceed via intermittent motion, called 
avalanches. Here we develop a field-theoretic treatment to calculate, from first principles, the space-time statis- 
tics of instantaneous velocities within an avalanche. For elastic interfaces at (or above) their (internal) upper 
critical dimension d > due (due = 2,4 respectively for long-ranged and short-ranged elasticity) we show that 
the field theory for the center of mass reduces to the motion of a point particle in a random-force landscape, 
which is itself a random walk (ABBM model). Furthermore, the full spatial dependence of the velocity corre- 
lations is described by the Brownian-force model (BFM) where each point of the interface sees an independent 
Brownian-force landscape. Both ABBM and BFM can be solved exactly in any dimension d (for monotonous 
driving) by summing tree graphs, equivalent to solving a (non-linear) instanton equation. We focus on the 
limit of slow uniform driving. This tree approximation is the mean-field theory (MFT) for realistic interfaces 
in short-ranged disorder, up to the renormalization of two parameters at d = due. We calculate a number of 
observables of direct experimental interest: Both for the center of mass, and for a given Fourier mode g, we 
obtain various correlations and probability distribution functions (PDF's) of the velocity inside an avalanche, as 
well as the avalanche shape and its fluctuations (second shape). Within MFT we find that velocity correlations 
at non-zero q are asymmetric under time reversal. Next we calculate, beyond MFT, i.e. including loop correc- 
tions, the 1-time PDF of the center-of-mass velocity u for dimension d < due. The singularity at small velocity 
'P(u) ~ l/it^ is substantially reduced from a — l(MFT)toa — 1 — |(4 — d) + ... (short-ranged elasticity) and 
a = 1 — |(2 — d) + ... (long-ranged elasticity). We show how the dynamical theory recovers the avalanche- size 
distribution, and how the instanton relates to the response to an infinitesimal step in the force. 



I. INTRODUCTION 

Elastic interfaces driven through a disordered medium have 
been proposed as efficient mesoscopic models for a number 
of different physical systems and situations, such as the mo- 
tion of domain walls in soft magnets d-lit], fluid contact lines 
on a rough surface ll9l-fTTll. or strike-slip faults in geophysics 
lfl2l - [l5ll . Their response to external driving is not smooth, but 
exhibits discontinuous and collective jumps called avalanches 
which extend over a broad range of space and time scales. 
Physically, these are detected e.g. as pulses of Barkhausen 
noise in magnets iH 0, [l6l - [l8ll . slip instabilities leading to 
earthquakes on geological faults llsl [l2[ [l9l - i22ll . or in frac- 
ture experiments I23I-33JI. While the microscopic details of 
the dynamics are specific to each system, an important ques- 
tion is whether the large-scale features are universal ||34| . The 
most prominent example are the exponents of the power-law 
distribution of avalanche sizes P{S) ~ (for earthquakes, 
the well-known Gutenberg-Richter distribution ll35l - B7ll ) and 
durations, which are believed to be universal. Beyond scaling 
exponents, the question of whether the shape of an avalanche 
is universal is of great current interest 03811 . Understanding 
whether and how universaUty arises, and obtaining quantita- 
tive predictions for avalanche statistics beyond phenomeno- 
logical models are some of the main challenges in the field. 

Historically, the elastic interface model has allowed for an- 
alytical progress thanks to a powerful method, the Functional 
Renormalization group (FRG). This method was first devel- 
oped to calculate either the static (equilibrium) deformations 
of an interface pinned by a random potential (e.g. the rough- 
ness exponent) [i39h421 . or the critical dynamics at the depin- 
ning transition which occurs when applying an external force 
/ > /c ll4ll I43I - I50I1 . These results are obtained in an expan- 



sion in the internal spatial dimension d of the interface, around 
the upper critical dimension due, in a loop expansion. Despite 
these successes the study of avalanches in elastic systems has 
remained centered on toy models fllsl [13I1 or on scaling ar- 
guments and numerics 10, |45l I5ll - l56ll 7 Several other impor- 
tant models have been used to describe avalanches, such as 
the random-field Ising model ll57l - l59ll and discrete automata 
known as sandpile models, for which analytical results exist 
ll60l - l67ll . However, exact results on the avalanche statistics are 
notably hard to obtain. 

One simplifying feature of the interface model in its basic 
version, i.e. with over-damped dynamics, is that it satisfies 
the no-crossing rule, or Middleton theorem, which guarantees 
only forward motion after a finite time, and uniqueness of the 
sliding state ll68l - l70ll . This allows to define unambiguously, 
at fixed driving velocity v, a quasi-static limit w = 0+ which 
we have studied with high precision both from numerics and 
using the FRG, testing the agreement up to two-loop accu- 
racy iItiIi . Recently, we have developed FRG methods Il72l - l76ll 
to calculate the statistics of avalanches for elastic interfaces, 
both in a static, and quasi-static framework, obtaining e.g. the 
distribution P{S) of their size, i.e. the total area swept dur- 
ing an avalanche. Initially our calculation focused on static 
avalanches, i.e. switches in the ground state. However, thanks 
to Middleton's theorem, it can be extended to quasi-static driv- 
ing: Since the system visits a unique sequence of metastable 
states, we define quasi-static avalanches in a stationary regime 
(for V = 0+) as jumps from one metastable state to the next. 
The avalanche size S depends only on the initial and final con- 
figuration, and is a property of the quasi-static limit. We found 
[73, 77] that to 1-loop accuracy P{S) is the same as for de- 
pinning as for the statics, although we expect them to differ at 
2-loop order 
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In this paper we extend our study to the dynamics inside an 
avalanche; we calculate the probability distribution of the in- 
stantaneous velocity during an avalanche. Although we focus 
on the small-driving-velocity limit, it is a truly dynamical cal- 
culation. To properly define the avalanche statistics, we found 
it important to separate two very different velocity scales: (i) 
the small driving velocity v, which allows to separate different 
avalanches and to define a stationary regime; (ii) the motion 
inside an avalanche, which is much faster than the driving ve- 
locity V, and independent of it for small v. It is this fast motion 
that we study here. 

To this aim, we consider the following over-damped equa- 
tion of motion, which reads, in its simplest form (for short- 
ranged elasticity of the interface), 

77o?i(x,t) = c'V'^u{x,t) + F{u{x,t),x) + m'^[w{t)-u{x,t)] . 

(1) 

Here and below, we denote indifferently by u{x, t) or 
dtu{x, t) the local interface velocity. The time-dependent 
scalar function t), x describes the displacement of 
a d-dimensional interface in a d + 1-dimensional system. The 
quenched random force x) can be taken as a Gaussian 
random variable, short-ranged in x-direction, but with arbi- 
trary correlations in u-direction, 

x)F{u', x') = 5'^{x - a;')Ao(w - u') . (2) 

In most applications, the disorder Aq{u) is a short-ranged 
function. The interface is driven and confined by a parabohc 
well of curvature m^, which advances according to 

w{t) = vt . (3) 

This model, and this type of driving, is of experimental rel- 
evance for the systems mentioned above. In some cases, it 
requires an extension of the elastic kernel to non-local elastic- 
ity, which amounts to replacing in Eq. ([T]i, in Fourier space, 

cq^ + ^ €{q) = g{q)-^ . (4) 

The combination e{q)\uq\'^ is the energy associated to the 
mode q, which includes the elastic energy plus the coupling to 
the quadratic well. We have defined its inverse g{q)^ i.e. the 
(static) propagator, which we use extensively below. One ex- 
ample is e{q) ~ c{q^ + ji'^)''^'^, or more complicated kernels, 
and we always denote g{q = 0) = 1/m^ and e{q) ~ q'* at 
large q. For a contact line, m is related to the inverse capillary 
length /i (usually called k), set by surface tension and grav- 
ity ITSll and 7 = 1. For a magnet, m is set by the so-called 
demagnetizing field S 01 and 7 = 1 in some situations 
dominated by dipolar forces, while 7 = 2 in others. In frac- 
ture experiments, e.g. when breaking apart two plates which 
have been sintered together ll23l - l27ll . is proportional to the 
inverse thickness of the plates, and usually 7 = 1. 

A toy model to describe the avalanche dynamics which re- 
sults from Eq. ([T]i has been proposed by Alessandro, Beatrice, 
Bertotti and Montorsi (ABBM) 10,11], and further developed 
in lisi Issi l79l - [8lll . It approximates the motion of the domain 
wall, i.e. a system with many degrees of freedom, by the mo- 
tion of a point, at position u{t), which satisfies the equation of 



motion 

T]u{t) = F{u{t)) - [u{t) - w{t)] . (5) 

In 0], the random pinning force F{u) acting on this point was 
postulated to be a Gaussian with the correlations of a random 
walk, 

[F{ui)-Fiu2)f = 2a\ui~U2\ , (6) 

where cr > characterizes the disorder strength. One of 
the motivations for this assumption was that the model be- 
comes solvable. Although a crude description, it was used 
extensively to compare with Barkhausen-noise experiments 
on magnets, with success in some cases (systems with long- 
ranged elasticity) and failures in others IsI-ItI [111 HH. The 
most natural interpretation is that u{t) may represent the av- 
erage height of the interface, u{t) — J d'^x u{x,t), and 
that the ABBM model gives a mean-field description of the 
elastic interface. The random force F{u) is then interpreted 
as an effective random force, sum of the local pinning forces 
in some correlation volume. This is in agreement with the 
remark IHI^ that for infinite-range interactions the effective 
disorder is indeed correlated as in (|6]l. Thus this view has been 
taken for granted for a while. However, until now, there was 
no derivation from first principles starting from the realistic 
microscopic model of an elastic interface. 

In this article, we go beyond this simple toy-model descrip- 
tion of avalanches, and consider the motion of an elastic inter- 
face given by Eq. ([T]). We use the dynamical field theory and 
methods from the functional renormalization group (ERG). 
Let us recall that the upper critical dimension is due = 27 
in general, hence d^c = 4 for short-ranged elasticity, and 
due = 2 for the most common long-ranged elasticity, i.e. mag- 
netic systems with dipolar forces, the contact line or fracture. 
In this article, we will show: 

(i) In the small driving-velocity limit, all correlation func- 
tions (in time and space) of the instantaneous velocity 
u{x, t) can be computed (to lowest order in v) in a di- 
mensional expansion around d^c- This is done by com- 
puting averages of exponentials of the velocities (gen- 
erating functions), whose 0{v) contribution allows to 
extract the full probability distribution of the velocity 
field u{x, t) during an avalanche. 

(ii) At the upper critical dimension d = d^c, and in the 
small-m limit, the velocity field in an avalanche has 
the same space-time statistics as the Brownian Force 
Model (BEM) with renormalized parameters -q ?/„i 
and a — > am- The BEM is a model for an interface de- 
scribed by ([Hi where F(ti, x) are Brownian motions in 
u, of variance a, uncorrected in x. It is a generaliza- 
tion of the ABBM model to a set of elastically coupled 
ABBM models. For the BEM the generating functions 
of the velocity are obtained exactly in any dimension 
d by summing only tree graphs. Furthermore one can 
consider that the "tree theory" is the correct mean-field 
theory and describes the system for d > d^c, with full 
universality at d — d^c and small m. 
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(iii) The ABBM model (O with the force-force correlator 
(|6]l correctly describes the avalanche motion of the cen- 
ter of mass of the interface for d > c?uc in the limit 
V = 0^ . Universality arises for d = due and small m, 
with a dependence of the effective parameters rj ^ rjm 
and a — > am that we computed. 

(iv) Even for d ~ due the original ABBM model is not suf- 
ficient to describe the velocity correlations of different 
points on the interface, or the statistics of Fourier modes 
q 0- The latter can however be obtained from the tree 
theory (i.e. the BFM) which we show to be equivalent to 
solving a non-linear instanton equation. From this we 
obtain e.g. the avalanche shape at finite qatd = d^c- 

(v) Finally, for d < due the velocity field in an avalanche 
has universal statistics not given by the BFM, nor, for 
the center of mass, by the ABBM model. It can be 
obtained within an e = due — d expansion. We show 
that the one-time center-of-mass velocity distribution 
diverges at small velocity not as P{u) ~ l/w, but with 
a modified exponent 

Piu) ^ 4 . (7) 

For short-ranged elasticity the exponent is (with e = 
4- d): 

2 

a = l- -e + 0(e^) non-periodic, RF (8) 

a = 1 — -g + 0{e^) periodic . (9) 

For long-ranged elasticity (7 = 1), the exponent is 
(with e = 2 - d): 

a = l- -e + 0(e^) non-periodic, RF (10) 
9 
2 

a = 1 — -e + ©(e") periodic. (11) 
3 

A short report of some of our results has already appeared 
as a Letter [82j|. The present study is the starting point of a 
calculation of the avalanche shape and duration to order 0(e) 

Since the methods used here (based on the dynamical MSR 
path integral) are quite different from the usual Fokker-Planck 
approach to solve the ABBM model iQ-Ht], our study also pro- 
vides a new way to solve the ABBM model. In particular, we 
find that generating functions can be obtained from the solu- 
tion of the non-linear instanton equation. This new connection 
has been exploited and extended in ifsill to derive new results 
for the ABBM model (and elastically coupled ABBM models) 
for finite v > Q and for a non-stationary avalanche dynamics. 

One should emphasize that the methods introduced in the 
present work strongly rely on the Middleton theorem. Al- 
though specific results are obtained for an over-damped dy- 
namics, the present methods can be extended to any dynamics 
which satisfies the Middleton theorem. As an example, we 



have recently studied the ABBM model in presence of retar- 
dation issil . A much greater challenge for the future would 
be to extend these methods to models where the no-passing 
rule does not apply, such as models with inertia or relaxation 
which have been proposed, e.g. to study earthquake dynam- 
ics ifssll . There the very existence of a quasi-static limit is 
much less clear, and may depend on details of the dynamics. 
Some steps in that directions have been taken in ll86ll . Finally, 
let us also mention related studies of static avalanches in spin 
glasses using Replica Symmetry Breaking iIstIIssIi . and in the 
Random-field Ising model lls^ . 

The outline of this article is as follows: 

In section |ll] we introduce the interface model, define im- 
portant observables, and explain our strategy for their calcu- 
lation. We also review the expected scaling relations for the 
avalanche statistics. 

In sectionHni we construct the theory at tree level. We start 
with calculating the moments of the instantaneous velocity 
in subsection IIII Al before introducing in subsection IIIIBI a 
non-linear equation, which we call the instanton equation, to 
efficiently resum them. In subsection IIII CI we calculate the 
joint probability distribution for the center-of-mass velocity at 
one and several times. From that we extract various velocity 
probability distributions, and calculate the average shape of 
an avalanche, as well as its variance which we call the second 
shape. In subsection IIII Dl we show that the solution of the 
instanton equation encodes the response to a small step in the 
applied force. In subsection IIII Fl we recover the quasi-static 
avalanche-size distribution. In Section lTlI Gl we discuss the re- 
lation between the tree theory and the mean-field theory: We 
show that the tree theory is equivalent to (i.e. is exact for) the 
Brownian force model, and, for the center of mass only, to 
the ABBM model. We also show that the so-called improved 
tree theory, i.e. the tree theory with renormalized values for 
the disorder and the friction parameters, is the correct mean- 
field limit (for d = due) of the underlying field theory to be 
discussed in the following section |lVl Our approach is based 
on the Langevin equation and on the MSR dynamical action; 
alternatively one can use a Fokker-Planck description, as is ex- 
plained in subsection llll G 41 It is this latter description which 
was introduced by ABBM |0, Ht] for a particle, but whose use 
seems to be restricted to the latter In subsection IIII HI we ob- 
tain a number of results beyond the center-of-mass motion, 
such as the local averaged shape following a local step in the 
force, as well as the spatial and time dependence of the second 
shape. 

In section UV] we study the loop corrections, for d < due. 
We explain the general framework in subsection II V Al before 
introducing a simplified theory in section HV Bl containing all 
the needed ingredients for the one-loop calculation. The latter 
is solved perturbatively in subsection IIV CI We then discuss 
in detail the 1-loop, i.e. 0{e), corrections to the velocity dis- 
tribution in subsection llV Dl We derive the necessary counter- 
terms in subsection IIV Fl The extension to long-ranged elas- 
ticity is detailed in subsection llV Gl 

The above theory was developed in terms of the velocity 
ii as the dynamical variable. In section |V] we discuss how to 
perform the same calculations using the more standard theory 
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in terms of the position u. While this is more involved, it 
avoids certain technical problems which may be present in the 
velocity theory, and confirms the validity of the latter 

Several technical issues are presented in appendiceslAltolPl 



II. MODEL, OBSERVABLES AND PROGRAM 

A. The bare model 

We consider an elastic interface of internal dimension d, 
with no overhangs, parameterized by a time-dependent real 
valued displacement (or height) field u{x, t) = u^t G M, with 
X <E M''. It evolves in presence of a random pinning force 
x) according to the simplest possible overdamped equa- 
tion of motion, 

VodtUxt^ / {g^^)xx'{wx't-Ux't) + F{uxt,x) . (12) 

J x' 

Here ryo is the bare friction coefficient and {g^'^)xx' is the elas- 
tic matrix, with propagator g^x' = Qx-x' and g{q) = = 
/j. e^'^^gx in Fourier space and we define the (squared) mass 
rn^ = g^J^Q- Every where we denote equivalently := J d'^x 
and /( / dt. The interface is driven by an external 
quadratic potential centered at position Wxt- The total external 
force acting on the interface is noted 




with ft = m'^wt for spatially uniform driving. Equivalently, 
for inhomogeneous driving, Wxt denotes the reference inter- 
face position in the absence of disorder and in the limit of 
very slow driving (hence this notation is useful in the statics 
and the quasi-statics). We focus on the case of local or short 
range elasticity g^^ := + m?, with an elastic constant set 
to unity by choice of units. We will however also give the re- 
sults for more general non-local elasticity, see the discussion 
after Eq. Q. We focus on a uniform driving at fixed velocity 
V, Wxt = wt ~ vt. This leads to Eq. ([T]i in the introduction. 

The pinning force is chosen as indicated in Eq. (|2|i, where 
Ao(w) is the microscopic (bare) disorder correlator and 
denotes disorder averages. For realistic disorder the bare dis- 
order correlator is smooth. Note that for the bare model, we 
always assume (unless stated otherwise) a small-scale cutoff 
in X, either a lattice spacing a, or that Ao(u) decays on a 
finite correlation length r[. This insures the existence of a 
Larkin scale 190], which produces a small-scale cutoff for 
avalanches. We denote 5*0 the small-scale cutoff on their size. 

The above model exhibits two important properties: Due 
to statistical translational invariance of the disorder and its 
(5-correlations in internal space, the model possesses the so- 
called statistical tilt symmetry (STS) which guarantees that the 
elasticity Qg is uncorrected by fluctuations (loop corrections), 
see e.g. IITSll for notations and some definitions in this section. 
The second important property of the model is the Middleton 



theorem': If the driving force ni^wt is an increasing function 
of time, Wt > (positive driving), and if velocities are all 
pos itive at t = 0, Ux.t=o > 0, then they remain so at all times 
16811 . In particular, for a finite interface (of size L), submitted 
to positive driving, all velocities become positive after a finite 
driving distance, and the memory of the initial condition is 
erased. 



B. Quasi-static observables 

In this paper we focus on the stationary state of the model 
with fixed driving velocity wt = vt, hence iixt = v. We 
focus on the small-velocity limit v — 0+, i.e. on the vicinity 
of the quasi-static depinning transition. At a qualitative level, 
it is expected that because of disorder, at scales larger than the 
Larkin length Lc, the interface is rough at all times, i.e. self- 
affine [uxt — Ux'tY ^ Ix— a;' p^, with the roughness exponent 
C = Cdop of the depinning transition ll9ll - l93ll . Because of the 
mass term, the interface flattens for scales jx — > Lm, 
{uxt — Ux'tY ^ with Lm ^ 1/m for local elasticity. We 
are interested in the universality which arises in the small -to 
limit, i.e. for L,„ ^ Lc- 

It is also expected that on scales larger than the Larkin 
scale, the motion is not smooth but proceeds by avalanches, 
i.e. the system jumps from one rough metastable state to the 
next one. Thanks to the Middleton theorem there is a well- 
defined quasi-static limit, i.e. a function Ux{w) such that for 
w = 0+ one has Uxt = Ux{wt) where Wt = vt is the position 
of the center of the quadratic well. The sequence of visited 
states is unique. The quasi-static process Ux{w) was defined 
in 19411 and studied numerically in li7lL [tsIi , see also for 
an experimental realization. Note that the process Ux{w) is 
different from u^J^*'{w) defined in the statics fTSl ] which de- 
scribes shocks, i.e. switches in the ground state-. However, 
there are close analogies, hence similarities in notations in this 
section and in Ref. 17311 . The quasi-static process jumps at a 
set of discrete locations Wi, i.e. 

Uxiw) ^^^S^x^iw - Wi) . (14) 

i 

We also consider the motion of the center of mass of the inter- 
face, denoted 

ut := L--" f Uxt . (15) 

J X 

For V = 0+, it converges to the quasi-static process ut — 



' For a model discrete in x, this is the case if (g~^)xx' ^ for ^ 
Then iixt > if fxt > and Uxt^ > at some initial time t;. 

2 ^stat jg jjjg minimum-energy configuration for a given w. In contrast, 
for a particle m(ui) is the smallest root of the equation m^w = F(u) — 
•m?u and, similarly, for an interface Ux {w) is the metastable state with the 
smallest Uxiw) for all x. 
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u{wt) for the center of mass, denoted 



u{w) = L 



-d 



Ux{w) 



(16) 



Here 5*^ is the size of the i-th avalanche. In the statics, the 
statistics of these shocks was studied in Ref. iItsIi . Here one 
can also define their size density (per unit w) as 



p{S) = poP{S) =Y.6iS~ S^Siw - w{) . (17) 



The probability distribution P{S) of the size is normalized to 
unity. Since one can show that lf8lll94ll 



m [w — uiw)] ~ fc{fn) 



(18) 



the critical force at fixed m, it implies u{w) w w, hence the 
process follows the center of the well, although with a delay. 
This shows that the total density po P^r unit w is related to the 
average size as 



(5) 



(19) 



where here and below = / ASP{S)f{S) denotes the 

(normalized) average of f{S). Note that the existence of a 
short-scale cutoff (and a Larkin scale) guarantees that po is 
finite, although itjnay diverge if these cutoff scales go to zero. 

As shown in ivs]] there is an exact relation between the sec- 
ond moment of the avalanche-size distribution and the cusp in 
the renormalized disorder correlator. 



(5^) _ -A'(0+) 
2(5) ~ to4 



(20) 



It defines the avalanche-size scale S„i, which behaves as 
Sm ^ m^('^+'»^ at small m. The definition of the renor- 
malized disorder correlator A(u) is recalled below and its 
salient property is that it is non-analytic, even if the bare dis- 
order is smooth. This relation holds in any dimension, for 
statics and quasi-statics, i.e. depinning (with, accordingly dif- 
ferent values for A'(0+) and the roughness exponents). The 
only assumption is that all motion takes place in shocks or 
avalanches, as in (fT4] i. which usually holds for small enough 
771 (see IstIi for a case where the contribution from the smooth 
part of u{w) is calculated explicitly). 

The convergence to the quasi-statics in the small- t; hmit oc- 
curs on time scales t,,, 5w/v where 5w ^ Wi+i — Wi is the 
typical avalanche separation, is called the waiting time (un- 
til the next avalanche). On the other hand, the motion inside 
an avalanche occurs on the so-called duration time scale 



< Sw/i 



(21) 



where z is the dynamical exponent at depinning. In this pa- 
per we always assume v small enough so that the order of 
scales is as given by Eq. ( |2TI ). i.e., the avalanche duration 
is much smaller than the waiting time between avalanches, 
so that successive avalanches are well separated. In practice. 
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FIG. 1: Schematic plot of the instantaneous velocity as a function 
of vt for different v. The area under the curve is the avalanche size 
hence is constant as u — > 0^. The quasi-static avalanche positions 
Wi are indicated. 



when L ^ Lm and at least for SR elasticity, it may be suffi- 
cient to ask that successive avalanches occurring in the same 
region of space be well separated, i.e. that (|2TI) holds when 
Sw is the typical separation of avalanches in the same region 
of space. The condition ( |2TI ) is equivalent to the condition 
Ljn ^ S,v, where is the correlation length near the depin- 
ning transition ll43l - l45ll . 



C. Dynamical observables 

Our aim is to obtain information about the dynamics in an 
avalanche. For simplicity we will first consider the 77-times 
(instantaneous) velocity cumulants ut^...Ut„ for the center 
of mass, and discuss space dependence later The important 
property about avalanches, and non-smooth motion in general, 
is that in the limit v ^ 0'^ 



ut,...ut„' ^ vfih, ...,tn) + Oiv^) (22) 
ii;;Zii;:^vf{h,...,tn) + Oiv^) . (23) 

This means that cumulants and moments are 0{v), and have 
the same leading time dependence. This is very different from 
a smooth motion, for which they would be 0(u"). Here we 
are considering times much shorter that the waiting-time scale 
6w/v, hence a single avalanche. The result (l22l) can be un- 
derstood as follows: The above cumulants are non-negligible 
only when all times are inside the same avalanche. When that 
occurs, the velocities are 0(u°), with a magnitude studied be- 
low. Let us suppose that the separation of the times ti is of the 
order of T. The above cumulants are thus dominated by the 
probability that exactly one avalanche occurs in a time inter- 
val of duration T (with T ^ Aw/v). This probability is in 
terms of the total avalanche density po 

Prob(one avalanche in[-T/2, T/2]) = povT < 1 . (24) 
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More precisely, one can establish the sum rule 

[-T/2,T/2]" 

(25) 

which is valid as long as pqvT <C 1. It comes from the fact 
that the total displacement L"^ J dt iit during the avalanche i 
is equal to its size, 



S, =L'^ j di in 



(26) 



It is clear from the above that the difference between moments 
and cumulants is at most of 0(?j^). The sum rule (l24b thus 
connects dynamical quantities to quasi-static ones. It provides 
a valuable consistency check for our dynamical calculations. 



D. Strategy 

Let us now summarize our strategy. We will calculate the 
velocity cumulants from perturbation theory in an expansion 
in the disorder Naively this expansion is in the bare dis- 
order Ao(u). To lowest order the n-times cumulant ( l22b is 
0(Aq^^) and, as we will see below, is obtained from tree 
graphs in the graphical representation of perturbation theory. 
For each n, 1-loop graphs only occur at the next order, i.e. 
0(Aq ), and so on for higher-loop graphs. Hence we start by 
examining the perturbation theory at tree level in the next sec- 
tion. We compute explicitly the lowest moments, and then 
show that there exists a much more powerful method, based 
on a simplified field theory, which allows to sum all tree dia- 
grams and compute directly the Laplace transform of the joint 
probability distribution P{utn ...jiit^) of the velocities at n 
times. 

In practice it is in fact more accurate to work with the renor- 
malized disorder We recall that the renormalized disorder 
correlator A(w) is defined in the quasi-static theory from the 
center-of-mass fluctuations as 



[u{w) — w][u(w') — w'] =L '^A{w — w') 



(27) 



The function A(u) depends on m, with A(u) = Aq{u) for 
m — > oo. At small m it takes the universal scaling form 



1 



e/2 



h = 



1 



(<Z2 + 1)2 



(28) 
(29) 



It is given here for SR elasticity. The rescaled correlator A(u') 

converges to a (non-analytic) fixed-point form A(u)) > 

K*{'w) = 0{e). Here C, is the roughness exponent at depin- 
ning. The rescaled correlator K{w) obeys, as a function of 
m, a FRG flow equation which was obtained at the depin- 
ning transition, together with its fixed points, to two loops in 
ll4lll48ll and checked in ItiIi where A (it) was measured from 
numerics. 



Since it is the renormalized disorder, which is small, we 
then reexpress the perturbative expressions of the velocity cu- 
mulants in terms of /S.{w) directly. Thus we generate an ex- 
pansion in powers of e for these quantities. The leading order 
is determined solely by tree graphs in the renormalized dis- 
order A (each cumulant being of order A"~^) and is valid 
for d = (iuc (to some extent it is also valid for d > c?uc, see 
the discussion below). This leads to the tree-level result for 
the velocity probabilities. Corrections to the tree-level result 
are obtained in the next section by adding the contribution of 
one-loop diagrams, i.e. the next order in A* = 0{e). 

In the remainder of the paper we will switch to the comov- 
ing frame, unless explicitly indicated. Hence we define for 
w = vt 

Uxt ^vt + ii^t , (30) 
where Uxt satisfies the equation of motion: 

vdtUxt = ^luxt + F{vt + Uxt, x) - iri^Uxt - m' ■ (31) 
Below we will denote u by u for simplicity. 

E. Expected scaling forms for avalanche statistics 

1. Size distribution 

The size distribution is by now the best known one. Let 
us first recall our previous results lf72l ItsIi for the avalanche- 
size distribution in the small-m limit, i.e. Sm ^ <5'o, where 
Sq is the microscopic cutoff, and Sm the scale of the large 
avalanches, given by Eq. ( l20b . For S S> 5*0, the size distribu- 
tion ^(5*) takes the form 



P,i,e(5) = P{S) = ^p{S/S^) 



(32) 



Depending on the dimension d, p{s) takes different forms: (i) 

for d = d,,r 



1 



20F 



S-3/2g-./4 



(ii) for d < due. 



p{s) 



— —s exp 
2^ \ 



B 



(33) 



(34) 



to first order in 0(e = duc-^), where ^-1, C 0(e) 

are given in iItsIi . The exponent i5 = 1 + -^[e — C) + O(e^) 
and the avalanche exponent 



(35) 



were agree to first order in e with the Narayan-Fisher (NF) 
conjecture SSUl, which relates the avalanche-size exponent 
and the roughness exponent via 



T = 2- 
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(36) 
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Here 7 = 2 for SR elasticity, 7 = 1 for LR elasticity^ and 
due = 27. This conjecture agrees well with numerics for 
d — 1,2,3 lf72l[75ll . both for the statics and quasi-statics (with 
the respective values for Q, but it is not known if it is exact 
(see the discussion in section VIII-A of ItsII ). It was pro- 
posed by NF for depinning only, but recently we have found 
a general argument for the statics as well, based on droplet 
considerations [82, ,8811 . 

Here we will recover the above results, within a dynamical 
calculation, to tree level in d = d^c, and to one loop, 0(e), for 
d < d,,r. 



2. Duration distribution 

Assuming one can define unambiguously the duration T of 
an avalanche (see the discussion below in Section lTlIG3b the 
duration exponent a is defined through the small-T behavior 
of the duration distribution 



1 



-Pduration(T) = ^/(T/To) , 



(37) 



where Tq is a large-time cutoff, and /(O) a constant. This form 
has been conjectured in various articles, see e.g. A simple 
scaling argument relates a to r and z, the dynamical exponent. 
One writes S - and T - L^ hence S - T^d+Q/z^ 

Then 



-fcluration (^) ^ -^size ^^S* ' 



dT 



T" 



(38) 



with 



(r-l)(d + C) 



(39) 



If we use in addition the NF conjecture ( [36b one obtains 

a = 1 H , 



(40) 



a relation which was conjectured previously, see e.g. H. It is 
not known at present whether these conjectures are exact. The 
methods of the present paper allow to determine Pduration(7')- 
Here we obtain it to tree level, and in llssll to one loop. 



3. Velocity distribution 

Here we obtain the distribution of velocities in an avalanche 
in the form 



Viu) ^ - 



(41) 



We will obtain a = 1 at the mean-field level (tree theory), as 
in the ABBM model, and a < 1 for d < due- It turns out that 
our result for the exponent a is not straightforward to derive 
from scaling arguments. Hence it may be a new independent 
exponent. 



III. TREE-LEVEL THEORY 

In this section we implement the program explained above 
to lowest order, i.e. at tree level. Hence we construct the 
proper mean-field theory for the interface. We will use sys- 
tematically the notation A for the disorder vertices and 77 for 
the friction. Hence, if one substitutes Aq and 770 one gets the 
naive perturbation result, i.e. genuine tree graphs. If one con- 
siders A and 77 as the renormalized disorder correlator and 
friction, one obtains the result using the so-called "improved 
action'^Le. the limit for d ~ d^c of the effective action (see 
Refs. iFtsI F74I1 for more details on these definitions). This 
amounts to summing tree graphs plus those loop diagrams 
which renormalize friction or disorder at d = due- Some- 
times we will denote A — > A^ and rj 7]„i to remind that 
these quantities are m dependent. In simple terms, the results 
expressed in terms of A^ and rjm are numerically accurate at 
d = due, with the correct, and universal, dependence on m for 
small m. 

It is useful to recall here the result of iItsI for the generating 
function and avalanche-size distribution at tree level. 



(c'^s - 1) 



(S) 



25„ 



1- v^l - 4A5„ 



(42) 
(43) 



We have added the subscript 5* to distinguish from the nota- 
tion for the dynamical generating functions introduced below; 
let us also note that we use indistinguishably the three suffixes 
"tree", "MF", and "0", to indicate tree, i.e. mean-field quanti- 
ties. Eq. (1431 ) holds both for the statics and quasi-statics, and 
will be recovered below in the dynamical approach. 



A. Calculation of moments 

The equation of motion dsTl ) in the comoving frame can also 
be written as 



Uxt — j Rxt,x't' 
I x't' 



F{vt + Ux't',x') - r]v 



where 



Rxt,x't' = Rx-x'.t-t' = {vdtStf + yxx'^^tt')^ 



(44) 



(45) 



is the bare response function with Rxt = Jq'^'^'^^Rqt- In 
Fourier space it reads 



1. First moment 



(46) 



^ The exponent 7 is often called /i in the literature, see e.g. 0|. 
In Section lTvl the notation a is used for a different quantity, see Eq. j2991 . 



We start with the first moment, which defines the critical 
force fc ~ fci'iTT', v)- Taking the disorder average of (l44b we 
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have 



fc ■■= F{vt + Uj:t,x) 

This yields the exact equation 



(47) 
(48) 



[ Rx-cc'.t-t' [F{vt + u,,t, ,x')-h], (49) 

Jx't' 



from which we now compute the cumulants to leading order 
in perturbation theory. 



2. Second moment 
To lowest order in A one finds from Eq. ( |49] l that 

Uxitl'^X2t2^ = / Rxi-x' ,tl~t'Rx2-x' ,t2-t" ^t"y) ■ 

Jx't't" 

(50) 

From this we obtain the cumulant of the center-of-mass veloc- 
ity. 



^ , 1 

ut^ut^ =L dt^dt2^ 



X / e 1 

I Sl<tl,S2<t2 
,,2 



(tl-Si) —(t2-S2 



A(v(si -S2)) 



V 



-L-^^ j e-^(^^+^^)A"(«(ti-i2-ri+T2)). 

Tl>0,T2>0 

(51) 

Let us now consider the limit of w 0+, and assume that 
A(u) has a cusp, i.e. 

A"(w) = 2A'(0+)(5(u) + A"(0) + 0{\u\) . (52) 

Then we find that 



ut,ut2^ = -2i-'*A'(0+) — 



e 1 



■(2Ti-ti+t2) 



T Jti>0 



T J 

ri>0,T2>0 



+Oiv^) 
Hence we obtain 



(53) 



■\t2-tl\ 



(54) 



Note that the cusp is crucial to get non-smooth, avalanche mo- 
tion: Since ii = v, the term of order v in the above equation 
is possible only since the manifold moves with velocity ii of 
order one (i.e. independent of v) for a time of order 1/v. In 
the absence of a cusp, ii ^ v, and the second cumulant of the 
velocity is O(u^) indicating a smooth motion. To this order. 



the typical time scale of an avalanche is read off from the 
exponential in the first line of Eq. ( |54T ). as 



1] 



(55) 



In the improved action, rj will be renormalized to 77 = -qm, as 
is discussed below. 

Using that the size of an avalanche is 5* = L'' iit, we can 
now integrate over the time difference to obtain 

dt = -2vL'^ ^'^^P . (56) 



Using Eq. ( fT9l ), i.e. po = L"^/ {S), this exact relation agrees 
with the general sum rule for n = 2, provided Eq. ( |20l ) holds. 
This is indeed an exact relation obtained both in the statics 
and in the quasi-static limit in ll73l[82ll : it relates the cusp to 
the second moment of the avalanche-size distribution. 

In order to simplify the notations for the calculation of 
higher cumulants, we now switch to dimensionless units. 
They amount to replacing 



X — > x/ni, 
and A'(0+) - 



L L/m, 
4-^'A'(0+ 



m 



t^tTm, V^v/t„i (57) 

. In effect this is equivalent to 



setting 7] = m~ = 1. 

We now reproduce the above result, introducing a graphical 
representation which will be useful for the calculation of the 
higher cumulants. Let us consider Eq. (l50l ) integrated over 
space and rewrite it graphically as 



L'^ut^ut2'' = 



t2 



(58) 



Here the dashed line represents the disorder vertex A which is 
bilocal in time and the fuU lines are response functions (|46l) . 
here taken at zero momentum q = 0. (For details on this 
standard graphical representation see e.g. I^H.) The second 
velocity cumulant thus reads 

tl t2 

L'^iit.ut," = dtA2 t 1 ■ (59) 



Hence the time derivatives act on the external legs. We now 
use the fact that the response function depends only on the 
time difference, i.e., 

dtiRq,ti-Si = —ds^Rq^ti-Si , (60) 

where here and below we denote Rt Rq=o.t = 0{t)e^^ 
in our dimensionless units. Hence, by partial integration, we 
can move both time derivatives to act on the disorder ver- 
tex as ds^ds2 which produces the term — u^A"(u(si — S2)) 
as in Eq. ( 1531 ). To lowest order in v this can be replaced by 
— 2wA'(0"'')(5(si — S2), hence the two internal times are iden- 
tified. This can be represented as 



^2i;A'(0+) 



-(2si-ti-t2) 



Si <min(ti ,(2) 



= -'yA'(0+)c-l*i~*=^l , (61) 
recovering the above result (|54T ) to lowest order in v. 
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Third moment 



We are now ready to compute the third cumulant. Here and 
below we label external times by ti and internal times by Sj 
(black dots). To lowest order in the disorder, one finds from 
Eq. gill; 



6Sym 



(62) 

where Sym denotes symmetrization w.rt. the external times 
ti . Hence one has 



1' 



(63) 

The first thing one could do is to perform the dt^ derivative, 
using partial integrations 



= / i?t3_,,,a,,A'(s3-s4) 

J S4, 

= -2A'(0+) j Rt,-sJis3~Si) 



(64) 



Note that we have safely replaced A'('i;(s3 — 54)) by A'(s3 — 
S4) since we anticipate that to lowest order we will need only 
A'(w) = A'(0+)sgn(u) + 0{u). Note that there is no bound- 
ary term if time integrals are performed from [—00, 00] and 
the theta function is included in R. By this procedure, the 
term A (53 — 54) will have exactly two derivatives. However, 
to be able to proceed further, it is better to consider 9*3 
simultaneously, while symmetrizing at the same time leading 
instead to (passing always one external derivative onto each 
disorder vertex-end): 



-A'(0+^ 



dt^R{t2 ~ S3)R{t3 - S4)'5(s3 - Si) 
+ R{t2 - Si)dt3R{h - S3)(5(s3 - S4) 



= -A'(0+) [dt., + dt,) [R{t2 - s3)R{h - S3) 
= A'(0+)a,3 \R{h - s3)R{h - S3) 



(65) 



Integration by part w.rt. S3 is then possible, and together 
with taking dt^ on the left branch and using time translational 



invariance of i?. 



and Rti-si respectively leads to two 



derivatives on the lower vertex A(si — S2). 

In summary, we find that the surplus external derivatives 
can always be passed down in the tree, so that at the end each 



vertex receives exactly two derivatives. This means that we 
can rewrite (1631 ) as 



L^^.ut.ut, =6v A' (0+)^ Sym f f 

J Si J S2 




where the points are intermediate times, and the arrows stan- 
dard response functions. We now have to compute this new 
diagram, with the huge simplification that vertices are now lo- 
cal in time and which apart from the vertices contains only 
response functions. 

We also note that the single v factor comes from the lower 
vertex: This can be interpreted as the point in space and time, 
where an avalanche is triggered with rate v. 

Let us now complete the integration over internal times. To 
this aim, let us fix the smallest internal time si, and integrate 
over S2- 




— Rti-si 



R, 



t2—S2 Rt3—S2 Rs2—Si 



g-[max(t2,t3)-Sl] _ g-(t2-Si)-(t3-Si) 



— Rti-Si 

xQ{si < mm{t2,t3)) . 
Integrating once more gives 

Iv 2v ^3 
'2 



(67) 




_ _g2 min(ii ,t2,i3) — max(t2,t3) 



^3 min(ii ,t2,*3) — *l— 12— *3 



(68) 



Finally, after symmetrization it simplifies into 



6 Sym 



jjmin(ti,t2,t3) — max(ii,t2,t3) 




(69) 



Hence, assuming that the external times are ordered as ti < 
t2 < ts we obtain our final result for the third velocity cumu- 
lant as 



(70) 



= 2wA'(0+)^e"(l*i~*"l+l*i"*^l+l*2"*-'*l)/^ . 

Note that the final expression is simple, while the starting one 
was quite non-trivial. 

We can now check that the sum rule ( |25] ) is satisfied. Indeed 



-2d 



= 12i;A'(0+)2 



UtiUt2Ut^ 



0=tl<t2<t3 



12wA'(0 



(71) 



recovering the result of 117311 . and which can be obtained by 
expanding ( |43] ) for the third moment of the avalanche-size dis- 
tribution. 
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4. Fourth moment 



The sum rule gives 



The higher moments can be computed using the same 
method, as the same simplifying features can be generalized. {S ) 

The result for the fourth cumulant is, supposing the times are (5) 
ordered as ti < t2 < ^3 < ^a'- 



-24uA'(0+)^Sym 




(73) 



ut.ut^ut^uu = 120w|A'(G+)|3 



t2-t3 .t4 



(74) 



which coincides with the result for the fourth moment of 
(72) Eq. m- 



5. Fifth moment 
Finally, we give the fifth moment 



= i'A'(0+)''5!Sym 



7- 4fl~ = : '■ = 

r 5^ 4v 3v 2v ^ I IV 3> 

<l 

02 +2 

'3 

4 •4 

?;A'(0+)''[8C*^"*'^ _^4gtl+t2-t3-t5 _^8gtl+t2-t4-t5 ^ 4(,tl+t3-t4-t5] 

I 






(75) 



We check that 



{S* 



(S) 



= L 



id 



ut^iit^ut^uuuu 5!x 14uA'(0+)* 



(76) 

coincides with the result for the fifth moment of (l43l l. 

The above results suggest that there is an underlying sim- 
plification at the level of tree diagrams of the original field 
theory, which is non-local in time, into a field theory which is 
local in time. We now show how the latter arises. 



B. Generating function and instanton equation: Simplified 
(tree) field theory 

Since here we want to study the temporal and spatial statis- 
tics of the instantaneous velocity field, we define the fol- 
lowing generating functional of a (possibly space- and time- 
dependent) source field X^t, 



G[\] :== G^t Ax*(f+ux*) _ 



(77) 



We remind that we are working in the comoving frame, i.e. 
V + iixt is the velocity of the manifold in the laboratory frame. 
The functional G[A] encodes all possible information. In par- 
ticular, all moments can be recovered by differentiation w.rt. 



the source. In this article we focus on the small driving- 
velocity limit. In view of the results of the previous sections, 
it will be sufficient to compute the generating function. 



Z[\] := L-'^dyG[X] 



v=0+ 



(78) 



which contains the leading 0{v) dependence of all moments 
in the limit of small velocity w = 0+. 

It turns out that, within the tree level theory, it is possible to 
compute these generating functions and obtain all cumulants 
at once, as well as the velocity distribution. We now show 
how this simplification occurs. 

We start not from the equation of motion ([T]), but from its 
time derivative in the comoving frame^ 

i7jdt-Vl+m^)u^t^dtF{vt + u^t,x) + f^t-m^v. (79) 

For completeness we wrote it for arbitrary driving f^t ~ 
{m? — V\)wxt, however we will mostly specialize to uniform 
driving, i.e. w^t = v, fxt = m^v, in which case the last term 
is zero. We denote indifferently time derivatives by ii or dtu, 
and for now we use the original (microscopic) units. Again, 



' Below, when indicated, we will alternatively use this equation in the labo- 
ratory frame, which amounts to setting v = Oin Eq. il9\ . 
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one has to set 77 770, A Aq for a derivation starting from 
the bare model, or the renormalized parameters if one deals 
with the improved action. 

We now average over disorder (and initial conditions) using 
the MSR dynamical action S associated to the equation of 
motion (|79] i: 

S = Sa+ Sdis (80) 
So= [ u^tivdt - ^1 + m^)u^t (81) 

Jxt 

5dis = / UxtUxt'dtdt'A{v{t-t') + u.j,t-Ua;t,). (82) 

^ Jxtt' 

Note that this is the dynamical action associated to the veloc- 
ity theory, i.e. in terms of Uxt and iixt to be distinguished from 
the one usually considered, associated to the position theory, 
in terms of Uxt and Uxt, to be discussed below. 

The generating function (iTTb can then be written as 



G[A] = / V[u]V[u]e 



5a = 5 - / Xxt{v + Uxt) 

I xt 



(83) 



(84) 



with G[0] = 1 and Z[0] = 0, since the dynamical partition 
function is normalized to unity. We can rewrite for the time- 
derivatives appearing in Eq. 



dtdt'A{v{t -t')+ Uxt - Uxt') 

= {V + Uxt)dt'A'{v{t - t') + Uxt - Uxt') 

= {v + Uxt)A'{0+)dt'Sgiiit - t') + ... . (85) 

Here we have used that sgn(w(i — t') + Uxt — Uxt') = 
sgn(t— i'), i.e. the motion for w > is monotonously forward, 
as guaranteed by the Middleton theorem ll68l . The neglected 
terms in Eq. dSSl l are higher derivatives of A(m)|„^o+- As we 
discuss below at length, they contribute only to 0{e = 4 — d) 
to Z[}], hence they can be neglected at tree level. This is 
consistent with our findings in the previous section that only 
A'(0+) appears at tree level. Hence we can rewrite the dis- 
order part Sdis of the dynamical action, which is a priori non- 
local in time, as S^is ~ + where 



= A'(0+) / uxtuxtiv + iixt) 



(86) 



is an action local in time. Furthermore we recognize the cu- 
bic vertex which generates the simple graphs obtained in the 
previous sections by a systematic perturbation expansion. The 
action 



5* 



ctrcc 
" '-'dis 



(87) 



is the so-called tree-level, or mean-field, action. Note that if 
we use the improved action, it then includes the loop correc- 
tions to 77 and A, and yields the correct result for d = d^c = 4, 
making the dependence in m explicit as 77 — >^ 77,„ and A — > 
A,„, see the discussion below and in Ref. iFtsIi . Note that due 
to the STS symmetry mentioned above, m^, the elastic coeffi- 
cient in front of W'^Uxt, and v are not corrected. 



We can now study algebraically the tree approximation 



•u=0+ 

G*'™[A] = J 2?[7i]2?[7i]e-'5"°° 



^tree 



Xxt{v + Uxt) 



(88) 
(89) 
(90) 



Note that the highly non-linear action dSTl ) ( [82b has been re- 
duced to a much simpler cubic theory. Cubic theories among 
others describe branching processes, such as the Reggeon field 
theory jgsll for directed percolation. The present theory how- 
ever is simpler, and can be reduced to a non-linear equation as 
we now explain. 

Remarkably, considering (|88] |. one notices that iixt appears 
in 5]^'°° only linearly, i.e. in the form J^^UxtOxt[u, X]. It 
can thus be integrated out, leading to a (5-function constraint^ 
Tixt ^{Oxt[u, A]). Hence in the tree-level theory the field iixt 
is not fluctuating, but given by the non-linear equation 



ivdt 



nn?)uxt 



A'(0+)7i2^ + A,t = 



(91) 



This equation is the saddle-point equation w.rt. ii of 5^''°° in 
presence of a source, and is satisfied exactly. We also call 
it the instanton equation. We denote u^^ the solution of this 
equation for a given source field Xxt with ii^^^ = 0. After 
integration over iixt, we thus obtain from Eqs. (|87] | to ( |90l i: 



G[A] = 



(92) 



Z[X]=L-' / [A.,-A'(0+)(7l^,f 

Jxt 

^-L-" [ {Tjdt + Vl~m')u' 

Jxt 



xt 



m'L-" I ul. 



(93) 



Here we have used the saddle-point equation (l9Tl i and, in the 
last equality, assumed that u^^ (resp. V xU^) vanishes at large 
t (resp. x). This is insured if the source vanishes at infinity 
which we assume in the following. Note that since Z[X] is 
independent of the velocity, Eq. (|92] l gives the full dependence 
at finite v, a fact which is exploited and studied in detail in Ref. 



In summary we find that the calculation of Z[X\, i.e. of all 
cumulants of the velocity field, is equivalent to solving the 
non-linear equation (|9T] ). The solution v/^i can be constructed 
perturbatively in an expansion in powers of the source Xxt- To 
lowest order 



XxtR 



xt-^xt^x't' 



0{X^) 



(94) 



x,t 



* Equivalently one can view iixt as a response field associated to tlie equa- 
tion Oxt [«, A] = 0. 
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where Rxt,x't' is the usual bare response function (l45T l. Inte- 
grating Eq. (1991 ) or ( |94] |. one finds 



Z(A) = L 



(95) 



which is consistent with iixt = (v is uncorrected). Pursuing 
to O(A^) and higher orders, one recovers the velocity cumu- 
lants obtained in the previous sections, and in addition obtains 
their full spatial dependence. Instead of working perturba- 
tively, we obtain and analyze in the next subsection the (joint) 
probability distributions of the velocity at one (and several) 
times, focusing on the simplest observable, the center-of-mass 
velocity iit- 

Let us note that the simplified (tree) theory defined above 
does not contain all tree graphs. There are other tree graphs 
involving A"(0) and higher derivatives, as e.g. the following 
configurations of order w^. 



+ 



+ 



(96) 



While they are similar to those in Eq. (1581 ). different classes of 
trees appear starting at the fourth moment, as e.g. 



(97) 



These diagrams are characterized by the fact that they have 
two (or more) roots (lowest vertices), and are of order (or 
higher). The full tree theory is studied in section|V]and can be 
reduced to two non-linear saddle-point equations. However 
since these additional tree graphs lead to contributions which 
are of higher order in v, to study a single avalanche in the 
small-iJ limit, they are not needed. 

Finally, it is important to stress that the above simplified 
tree theory corresponds to the problem of an elastic manifold 
in a random-force landscape made out of uncorrected Brow- 
nian motions, for which it is exact for monotonous driving. 
This is the BFM, discussed in Section lTlIGI 



C. Joint probability distributions for the center-of-mass 
velocity 

To analyze the results, it is convenient to use dimensionless 
equations, hence replacing x x/m, and t Tmt. In mean 
field Tm = ri/m^, A Xxm/Sm, u^t -> Uxt/im^Sm), and 
V — >■ vvm, where f„i = Smm'^/T„i, L -> L/m. We start 
by using these units and, whenever indicated, switch back to 
dimension-full units in discussing the final results. We also 
keep the factor of L'^ in the beginning, but later on we find it 
convenient to suppress it. That amounts to a further change of 
units asv^ vvm with Vm = {mL)^'^Vm whenever indicated 
below. 



1. 1-time center-of-mass velocity distribution 

The center-of-mass velocity distribution is obtained by 
choosing a uniform A^^f ~ At. The 1-time probability is ob- 
tained from the inverse Laplace transform of Z{\), choosing 
At \5{t). 



Z{\) = L-''9„ei"^(''+«)|, 



=0+ 



(98) 



Here ii = Ut=i), and the tilde on Z{\) reminds us that we use 
dimensionless units. The saddle-point equation ( |9T| l admits a 
spatially uniform solution Uxt — uu thus we need to solve 

{dt - l)ut + u1^ -\5{t) . (99) 

The boundary condition is ut at i = ±oo, leading to 

A 



Ut 



A + (l- A)e-* 
This gives the generating function 



6{-t) 



z{\)^ jut = -Hi-\) . 



(100) 



(101) 



We now want to infer from this the 1 -time velocity distribution 
in an avalanche. Before doing so, let us restore dimension-full 
units. We assume that in the limit u = 0+ there are times when 
the velocity is exactly zero, i.e. v+ii = (since we use the co- 
moving frame) and times (when an avalanche is proceeding) 
when the velocity is non-zero. This picture is confirmed by 
results below ^. Hence the 1-time velocity probability (at say 
time t = 0) must take the form 



P{u) = (1 - pa)5{v + U)+ PaViil) 



(102) 



Here pa is the probability that t ~ belongs to an avalanche, 
and Viii) is the conditional probability of velocity, given that 
t = belongs to an avalanche. Both V and P are normalized 
to unity. One notes the two (always) exact relations {u)p ~ 
and Pa{v + u)-p — v. It is easy to see that 



Pa = Pov{t) 



(103) 



The mean duration of an avalanche is (t) = r,; where 

Na is the total number of avalanches and the duration of the 
i-th avalanche^. Now from Eq. ( 1102b one has 



(104) 



' This gives the universal regime for u 3> f o ■ For velocities smaller than 
the cutoff vq one expects a dependence on the details of the dynamics. 

^ Note that we are imphcitly working to lowest order in ii, at small v. Hence 
the fact that Pa increases linearly with v. while (r) remains constant, does 
not conflict with the requirement that pa < 1, since we study here the 
regime of small pa. At larger v. avalanches will merge, and formula | |1031 
ceases to be valid. 
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Taking a derivative w.r.t. v, one obtains to leading order in 

i; = 0+ 



Z{\) 



The identity 

Z{\) = -ln(l - A) = 



— — Z{m '^VmX) 



- 1 



dx 



(105) 



(106) 



allows to perform the inverse Laplace transform^ of Eq. (IIOII ). 
We thus obtain, in the slow-driving limit, the distribution of 
the instantaneous velocity of the center of mass for vq <C u 
(where wq is a small-velocity cutoff) as 



V{u) 



1 



1 

— e 

X 



(107) 



We have defined i'm ~ {mL)^'^Vm = L~'^Sm/Tm- 
This agrees with the above exact relation which becomes 
Po{t){u)'p = 1 in the limit of v = 0+. One notes that the 
distribution of small velocities diverges with a non-integrable 
l/ii weight. Since 7'(u) should be normalized to unity, the 
ensuing logarithmic divergence requires a small-velocity cut- 
off vq. This leads to the additional relation 



P0{T)Vm 



In 



"0 



(108) 



Hence we already anticipate that the average avalanche dura- 
tion will exhibit a logarithmic dependence on the small-scale 
cutoff, as confirmed below. Let us note that the rescaled func- 
tion p{x) is not a bona-fide probability, rather it is normal- 
ized such that /dxxp(x) = 1. Finally let us comment on 
the typical scale of the center-of-mass velocity, Vm- Since 
iit ~ L^'^ jj. Uxt we find that the scaling variable x enter- 
ing p{x) is the ratio of the instantaneous increase in the total 
area swept by the interface, iixt, divided by its typical value 
Sm/Tm (hence it does not contain the factor of L~'^). 

Let us indicate here for completeness the 1-time instan- 
ton solution in dimension-full units, as well as the generating 
function: 



G(A) 



1 



-u. 



dimlcss 
t 



ril?Sr)_ 
vm^L'^ f dt ut 



{t/ Tjji , XSni I Tjyi ) 



'-L''Z{\S,„/t,„) 



(109) 
(110) 



We recall that ( 11071 ). and all formulae concerning the center- 
of-mass velocity distribution, assume that the driving velocity 
V is small enough at fixed L so that only a single avalanche 
occurs, Pq ^ 1; hence v scales as ~ L~'^. If L goes to infin- 
ity first, at fixed small v, multiple avalanches occur along the 



'in practice one performs the Laplace inversion on Z'{\) which yields 
u'P{u), thus has no singularity at ii = 0. 



0.5 




Ja 


- 


1 1 


- 


/ / ^'^ 




^ 0.1 


'''''''' 







FIG. 2: Solution ut of the instanton equation J99t as a function of t 
for a source A(t) = \\i(t — t\) + A25(t — t-z), with t\ — 0, and 
t2 = —1 < ti. The function ut has the following properties: (i) It 
has the form u{t) — ^^-t_i on any interval where the source X(t) 
vanishes, (ii) It is zero for for t > tiby causality, (ii) It jumps by Ai 
(here -0.5) at t = ti and by A2 (here -0.155) at i = t2- 



interface. For small enough v they occur at far away locations 
(distances 1 /m) and are statistically independent. In that 
case the center-of-mass velocity distribution can be computed 
from convolutions of the distribution (1107b . It tends to a Gaus- 
sian distribution for large L and fixed v. The present results 
thus describes mesoscopic fluctuations. 



2. Exact result for the p-time generating function 

We now obtain the generating function for the p-time dis- 
tribution of the center-of-mass velocity. 



v=0+ 



(111) 



by solving Eq. ( |99] l in presence of the source At = 
Y!,^=i ^ ^j)- In this subsection we order the times as 

< ti, although in the following 



t 



P+i 



-00 < tp < 



subsections we will choose the opposite order 
The solution reads 



Ut 



p e{t,+i < t < t,)u^- 

^ (1- V 

3 — ^ 3 3 



(112) 



with tp^i = —00 < tp < ■ ■ ■ < ti, u^- = Xj + and 
u,+ = 0. Integration of (I112l l leads to Zp Zp{Xi, Ap) 



with 



ZJp 

We used the definition 



^ln(l - Zj+l,jU^-) 



l_e-l*--*.l 



(113) 



(114) 
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hence in this section = 1 — e*' *j with i > j. To generate 
Zp one can construct a recursion relation for the argument of 
the logarithm. From the above, one finds 



n,+i = A,n,+i?,n,_i 



A. 



Zj+2j + lAj + l 



1 



with Ho = 1 and Hi = 1 — Z21A1, so that 



Inn, 



^ 1 ^1 ; 

here tp+i is set to —00. This leads to 

Z2 = — ln(l — Ai — A2 + A1A2Z21) 
Z3 = - In (^1 - Ai - A2 - A3 

+ ^i^jZij — AiA2A3Z322;2i 



(115) 
(116) 

(117) 



(118) 



(119) 



(120) 



By inspection of the higher-order results, we arrive at the fol- 
lowing conjecture for tp < ... < ti 



Zr, 



In 1-J]A, (121) 



13=2 l<ii<i2<...<iq<p j = l 

Note that this expression corrects a misprint in an earlier ver- 
sion of the result, Eq. (17) in jsJl- This can also be written 
as 



p-i 

= - In (1 - ^(-l)'=tr(iVAf'=) 



A;=0 



-ln(l-tr(iV(l + A/)~^)) 



A, 



U > J) 



(122) 

(123) 
(124) 



The functions Zp possess an interesting factorization property, 
which we demonstrate on the simplest example Z3; Suppose 
that we choose A2 = = — (i_Ai)e*i-'->+Ai ' 

finds that itt = in the interval < t < t2. This leads to 



A2 = - 



(l-Ai)e*i-t2+Ai 



^3(Ai,A2,A3) ^ =Z(A3)Z(Ai,A2) , 

(125) 

which we have checked explicitly. It implies that the observ- 
able e^^"'2+^i"ti for this particular relation between A2 and 
Ai is strictly statistically independent from the velocity at any 
time in its past. It would be interesting to investigate further 
the consequences of this property. 



3. 2-time probability 

Here we consider the joint velocity distributions at two 
times, and choose ti < t2 (from now one we choose the nota- 
tions of times in the more natural order ti < fi+i). We expect 
that in the limit w 0+ the 2-time probability takes the form 
(with iij iit.y. 

P(ui, U2) = (1 - gi - (72 - qi2)S{v + ui)6{v + U2) 

+ q25{v + Ul)V2{u2) + qi2'P{ui,U2) 

+qiS{v + U2)Vi{ui) . (126) 

The four terms, in the order of their appearance, are plotted on 
Fig. [3] The expression qi2 ~ vqi2 is the probability that both 
ti and t2 belong to an avalanche (case (iii) of Fig.[3]l. In the 
small-?; limit we are studying here, it must then be the same 
avalanche, and qi2 must be proportional to v. The quantity 
U2) is the normalized velocity distribution, conditioned 
to that event, qi ~ vq[ (resp. q2 = vq'2) are the probabili- 
ties that ti (resp. t2) belongs to an avalanche but not t2 (resp. 
ti), and7'i(?ii) (resp. 7^2 ("i)) the distribution conditioned to 
that event, (cases (ii) and (iv) of Fig. [3]). The first term in the 
decomposition ( 1126b ensures that the probability is correctly 
normalized. 

Integrating over 112, one recovers the single-time distribution; 
hence comparing with Eq. (11021 ) we have 

Pa = qi + qi2 = q2 + qi2 , (127) 

PaViui) ^ qiViiiii) + qi2 J (0121^(111, U2) , (128) 

and similarly for iii. Hence, qi = q2. From the definition 
( lllll ) of Z2 = Z2{Xi, A2) and Eq. (1126b we now have 



v=0+ 



Z2 = a„e^i("+"i)+^2(i'+ii2) _ 1 
= q[ y"d?ii7'i(?ii)(e^^"^ -1) 
+92 / du27'2(^2)(e^^"^-l) 



q'i2 / duidu2V{ui,U2){e 



1 ■ (129) 



We remind that here and below (until stated otherwise) we 
have suppressed all factors of L'^. The latter are restored be- 
low, when going to the result in dimension-full units'". Note 
that the symmetry of Z2(Ai, A2) in its arguments further im- 
plies that Vi{u) = V2{u) and that V{ui,U2) is also a sym- 
metric function of its arguments. Hence there is no way to tell 
the arrow of time from the velocity distribution of the center 
of mass at the mean-field level. Below we will however show 
that an asymmetry in time arises for finite Fourier modes, or 



Units of the center-of-mass velocity are then Vm which does contain the 
factor L~ , see the remark at the beginning of section lmcl 
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(i) 



(ii) 



(iii) 



(iv) 



FIG. 3: The four cases in Eq. ( I126t : both times outside the avalanche (i), only t2 inside the avalanche (ii), both times inside the avalanhe (iii), 
only ti inside the avalanche (iv). 



local velocities, already at the mean-field level. As a conse- 
quence, it will also arise for the center of mass at one-loop 
order JsH, i.e. for d < due- 

Taking now one derivative w.rt. Ai of ( |129t , one obtains 
from the formula (|119l l for Z2 via Laplace inversion the com- 
bination 

ill [q['Pi{ui)S{u2) + q[2'P{ui,U2)] 



lt: 



^dx,Z2i\uX2) 

1 + S2Z 



1 + SI+S2+ S1S2Z 

1 



"l(l + °2) 

LT~_^ e i+"2 



(130) 



We denote z := = 1 — e 1*^ We now use the general 
result 



LT 



-1 „d+ 



\/^/i(2V^)e'^-^" (131) 
V u 



with d = —u\jz, a = — z)/z'^, b — 1/z, and Ii the 
Bessel-/ function. This yields the smooth part, in dimension- 
less units, as q[2'P{ui,U2) = p{ui,U2) with 



,. . s -"i+"2 Z ( r^^\J\ - z\ 

P2[Ul,U2)=e " j== Ii[2^/UiU2 

Z\'UlU2 \ ^ J 



In dimensionfuU units 



(132) 

(133) 
(134) 



units) 



ai,(e^i("+"i) - l)(e^2(i'+«2) _ 1) 

= ^2(Ai,A2)-Zi(Ai)-Zi(A2) 

, / 1 - Ai - A2 + ZA1A2 

= — m 



v=0+ 



(1-Ai)(l-A2) 
duidu2p(ui,U2)(e^i"i - l)(e^^"^ - 1) (136) 



which is indeed satisfied by the function ( 11321 ). 

The (5-function piece in ( |131t allows to obtain (72 "7^2 (^2) in 
( 11301 ) (in dimensionfuU units) as 



(137) 



Normalization leads to qi ~ {v/vm)^n{zVm/vQ), in agree- 
ment with the results ( 11351 ), ( 1103b . (11081 ) and the sum rule 
( 11271 ). Note that ( 11371 ) can be obtained directly from Laplace 
inversion (in dimensionless units) of liniA2-i.-oo c'ai ^2 = 
z/(l — zXi) since that limit selects" the d{u2) piece in ( 1130b : 
equivalently, the first terms in ( 1 129b are 



q[ / dwi7'i(?ii)(e-^i"i - 1) = -ln(l - zAi) 



(138) 



Finally (1128b follows from the trivial identity Z2(Ai,0) 
^i(Ai). 



4. Avalanche duration 

The distribution of avalanche durations can be obtained by 
several methods. Let us recall that avalanche durations are 



Since V{ui, U2) is normalized to unity, integrating Eq. (1133b 
over both variables, one obtains the probability that both ti 
and t2 belong to an avalanche, 

qi2^vq[2^ ^\ii{l/z) . (135) 

For consistency we can check that the combination which in- 
volves only q'^.2'T^{'^i ; ^2) leads to a relation (in dimensionless 



Recall that the Laplace transform /(A) = LTu^_;\=s/(u) : = 
/dile^"7(?i) satisfies: (i) /(A) = 1 for f{u) = 5{u), (ii) /(A) - 
/(O) = -ln(l - A) for /(«) = sl± and (iii) /(A) = jj^^^^ for 
/(«) = ti°'e~", a > —1. Second, the behavior of f(u) at ii near zero 
is related to the behavior at A — > —00 of /(A): if the limit of A — > —00 
in /(A) exists, and is non-vanishing, it picks out the term ~ 5(u). The 
term /(0+) is extracted, in the same limit, from the term ^ l/(— A) in a 
large-A expansion. 
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well-defined as time intervals where the velocity is strictly 
positive. Consider then the probability that there exists an 
avalanche starting in [ti,ti+dti] and ending in [t2,^2 + di2]- 
On the one hand, this is equal to 



P(ti,t2)dtidt2 



duration 



(t = - ti)drd 



(139) 



where PdurationCT) is the probability distribution of avalanche 
durations. On the other hand it also equals 



dti dt2 9*1 9*2 512 



(140) 



where qi2 computed above is the probability that ti and t2 
belong to the same avalanche. From Eqs. (11351 ) and (11141) we 
obtain the distribution of durations as 



Pduration(^) 



(1-z) 



(141) 



where we recall z = Z21 = l — e 1*^ and in dimensionfuU 
units 



Pd 



uration 



(r) 



1 



r/T 



1 1 



(142) 



This probability distribution has a power-law divergence for 
small durations t T,n, 



duration 



(r) 



1 



(143) 



i.e. there are many short avalanches. We assume a micro- 
scopic cutoff time tq. The mean duration exhibits a diver- 
gence, i.e. 



1 



■ In ( ^ 

To 



(144) 



as a function of tq. However, higher moments are well- 
defined (i.e. independent of short scales). The expression 
(1144b is in good agreement with our previous result (1108b if 
one assumes ln( — ) « ln(— ). 

There are several other ways to obtain the duration 
distribution. First one notes that performing the limit 
limA2^.-oo dt-i constrains the avalanche to end at i2, and sim- 
ilarly — limAi-).-oo dtx constrains it to start at t\. Hence, in 
dimensionless units one recovers 



lim 

Ai , A2 — ^~ 



~dt^dt^Z2 = Pduration(i2 " ^l) ■ (145) 

It also yields another method to obtain from ( 11401 ), writing 

q'l2 = / dsi / dS2 Pduration(s2 - Si) 



lim 

Ai , A2 — — oo 

-ln(z), 



Z2(Al, A2,tl - t2) - Z2(Al, A2,oo] 



(146) 



inserting Eq. (11451 ) (second line) and Z2 from Eq. (11 191 ), in 
agreement with Eq. (1135b in dimensionless units. 

Another way to obtain the duration is as follows: We note 
that when the avalanche starts and ends, the velocity must van- 
ish. Hence the duration distribution can be recovered from 
7^(0+, 0+) which should be proportional to the probability 
that an avalanche starts at ti and ends at t2- We can indeed 
check on our result (1132b , ( 11331 ) that 



'7i2^(0+,0^ 



POT„ 



duration 



(r - t2 - h) ; (147) 



hence this is true, up to a normalization. We note that this 
term can also be obtained as the coefficient of 1/ (Ai A2) in an 
expansion of Z2 at large (negative) Aj. 

To study the temporal avalanche statistics, it turns out to be 
more efficient to use two properties simultaneously: (i) iii = 
outside the avalanche, an event whose probability can be 
selected by taking the limit A^ —00; (ii) taking a on 
the generating function multiplies by u*., hence is non-zero 
only if ti belongs to the avalanche. Using these properties 
we will now show how to generate the p-times distribution 
of velocities inside an avalanche conditioned to start and end 
at some given times. In particular, we recover the duration 
distribution, from the normalization, and we compute shape 
functions, which are of high interest in view of experiments. 



5. 1-time velocity distribution at fixed duration and mean 
avalanche-shape 

We start with the information contained in the joint 3-time 
distribution, which can be obtained from Z3 in dl 19b . Choos- 
ing again ti < t2 < ts, and generalizing the form (11261 ), we 
expect that the joint distribution contains a piece 



^^9l3,2'^13.2(u2)'5(wi)(5(-U3) , 



(148) 



where vq'-^^ 2 probability that ti and do not belong to 
an avalanche while t2 does, and 7^13.2 (1*2) is the velocity dis- 
tribution conditioned to this event. From the above remarks, 
to obtain this piece we need to inverse-Laplace transform 



lim 9a, ^3 = T—r 

AijAs— > — 00 — A2 

2212:32 Z21 Z32 



Hence we find in dimensionless units 

g;3,2^13,2(u2) = -e-""^ . 



U2 



(149) 
(150) 

(151) 



Integration over 112, in presence of a small-velocity cutoff vq, 
leads to 



'7i3,2 = -Infewo 



(152) 



Taking two time derivatives we recover the duration distribu- 
tion 

« « ' _ dt.bdt^b 

— OtiOtaQlsa — — ^duratioiAT — — tl) , 



(153) 
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FIG. 4: "Pulse-shape": The normalized velocity at time t in an 
avalanche of duration T for T <C Tm (lower curve) to T S> Tm 
(upper curve). 



using that dt^dt^h ~ 0. We also find the distribution of the 
velocity at ^2 conditioned s.t. the avalanche starts at ti and 
ends at t^. 



P(m2|13) 



-^tl^t3[gl3,2^13,2(M2)] 
^'dmation(''' = H — ^l) 



This leads to 



P{U2\13) ^ U2b^ 



(154) 



(155) 



From this one obtains the shape function 

("2)13:= J dU2U2P{u2\l3) 

_ 4sinh(2^)sinh(2^(l-^)) 



sinh I 



(156) 



2t„ 



for a fixed avalanche duration r = ta — ii, denoting t = 
t2 — ti. We have restored all units in the last line. This form 
interpolates from a parabola for small r <C t™ to a flat shape 
for the longest avalanches (see Fig. |4|i. The result holds for 
an interface at or above its upper critical dimension, which 
previously was used 13 8.1 on the basis of the ABBM model. 

An alternative approach is to obtain ^3(0+, U2, 0+) from 
Z3(Ai, A2, A3). As discussed above, one needs to extract the 
coefficient of l/(AiA3) in the large Ai, A3 expansion of Z^. 
Hence we first need to calculate 



^2|13 



(A2) 



Aid 



A?d~ 



lim lim -^Z3(Ai,A2,A3) 

A3— ^ — 00 aA3 Ai— > — 00 clAi 

A2 (Z3I + Z21 (^32 - 1) - ^32) -2:31 + 1 

(231 - A22;2lZ32) ^ 
-I 2 

1 b 



2 sinh I 



6- A2 



(157) 



dimensionless units) gives 

LTI^^_^^Z2|i3(A2) = 



X P(U2|13) 



4sinh2((^i^) 

= P(M2|13)Pduration(il - ^s) , 

where P(u2|13) is given in Eq. ( 11551 ), and Pduration in 
Eq. (fT42l) . 



6. 2-time velocity distribution at fixed duration and fluctuations of 
the shape of an avalanche: The "second shape" 

We now derive the 2-time velocity distribution at fixed 
avalanche duration. For that we consider the term 
(5('iti)(5(it4)gi4^23'Pi4. 23(^*27 ""3) in the joint 4-time distribu- 
tion (with ti < t2 < < t4) which can be obtained from 
Z4. We recall that 

pr • 11^^ -dtA^ [gM,23^14,23(u2,M3)] 

P{u2, U3|14) = fTTTT? ^^^^^ 

is the 2-time velocity distribution at fixed avalanche dura- 
tion T = ^4 — ti. We expect, and will check below, that 
-9ti dtiq'u.23 = Pduration(r = ti - ti), i.e. comparing with 
(11531 ), the number of intermediate points does not matter 

The simplest quantity to obtain is the 2-time shape function. 
Indeed multiplying ( 1158b by U2U3 and integrating, one finds 



{U2U3)l4: 



4 A2=0,A3 



It is easy to calculate from ( I121l i 



(159) 



lim 

2:212:43 



9a29a3^4|;,^^0^^3^0 



z 



(2:322:41 - 2:312:42) 



41 

sinh'(i(t2 - h)) sinh'(l(^4 - h)) 
sinh^(i(t4 - ti)) 



= 4 



(160) 



Taking two derivatives in ( |159l l one finds a complicated ex- 
pression for (m2U3)i4 which however simplifies greatly if one 
forms the cumulant combination and uses the above result for 
the shape. Then both results can be summarized, introducing 
the function h{t) := 4 sinh(f;/2), as (in dimensionless units): 



("2)14 

(u2W3)l4 



h{U - t2)h{t2 - h) 

("2^3)14 - (u2)i4(m3)i4 

1 fh {U-t3)h{t2-ti) 



(161) 



(162) 



6 is defined in Eq. (11501 ). The inverse Laplace transform (in 



Hence the fluctuation of the shape has a simple expression, 
and it would be nice to measure it in experiments. We call 
this the "second shape " since it gives more information about 
the avalanche statistics than the usual shape, the average of 
the velocity. The second shape tells about the variability, i.e. 
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fluctuations of the avalanche shape. For t2 = ^3 one recov- 
ers the relation {u^Y ~ between second cumulant and 
mean of the single time velocity distribution (11551 ). Note that 
the second cumulant always starts quadratically in time near 
the edges. It is quite remarkable that the dimensionless ratio 



{u{t2?)u 



14 



(163) 



is independent of ^i, ^4- This is an important signature of 
the mean-field theory which should be studied in experiments. 
On figure |5] we have plotted 



C(t,T) := 



mu{-t)Y 



m) (H-t)) 



(164) 



ti = -T/2,t4=T/2 



It measures the correlations between the left and right part of 
the avalanche. 

One can go further and obtain the full 2-time distribution. 
For this one notes that the function q'l^ 23^14,23 is obtained 
(in dimensionless units) by Laplace inversion as (i = 2, 3) 

<7'i4,23"2W3:Pi4,23(w2,W3) 



LT-i,^Y lim d,,d,,zA 

\ Ai,A4^ — 00 / 



lt: 



Z2lZi3{z3iZi2 - Z32Zil) 



[zil + 332:312:43 + 52221(242 + 53232243)]^ 

(165) 



We have used the result (I121l i. The normalization is obtained 
by integrating'- the above ds2 ds3 leading to 



gi4,23 = «9l4,23 



V , 2:42231 
— In 

Vrn 241^32 



(166) 



This is the probability that there is an avalanche starting in the 
interval [ti, t2] and ending in the interval [t3, t4]. Indeed one 
can check for consistency that integrating the duration distri- 
bution ( 1142b we obtain 



t't2 

vpo / dt' 



dtp, 



duration 



= (167) 

Vm Z41Z32 



Laplace inversion of (11651 ) w.r.t S2 yields an expression equal 
to minus the derivative —dt of (1131b . with other values for 
a = U2a' , b,d = —U2d'. Finally we find 



9l4, 23^14, 23 (m2,U3) 

with 



V U2U3 V / 



-bv.3 



d' 



Z31 
221232 
231242 



■242 
232^43 



232241 



1 



221232243 



4sinh^( 



(168) 

(169) 
(170) 



There seems to be a non-commutation of limits, lience we need to take first 
tlie large- A limit. 




0.1 0.2 0.3 0.4 

FIG. 5: The velocity correlation C{t, T) of Eq. l fT64l l for T = 1 



t 

0.5 T 



This leads to the final expression 

914,23^14, 23 ("2, U3) 



(171) 



1 



VU2U3 



Xe~( l-^'l-'^ "'"e'3-t2 -1 l_et2-'3 c*4 -"ts - 1 ) 

The 2-time velocity distribution at fixed avalanche duration 
T = <4 — is then obtained as 

P(w2,W3|14) = jp— 7 . (172) 



-19149^4, 



23 



This leads to the result 

P{U2,U3\U) 



(173) 



V"2M3 



\/u2U3 



^ 2sinh(^^) ^\sinh(ii^) 
sinh2(*i^) 
^ 4 sinh^ ( ) sinh^ ( ) 

One can check its normalization using the useful formula 



oo />oo 

dx / dy^/xy Ii{2ay^)e 



bx—cy 



{be- a'^y 



while derivatives w.rt. b and c allow to recover shape cumu- 
lants such as (11621 ). For instance one finds the third cumulant 
of the shape as 

{ulusY ^ {ulu3) - {ul){'il-3) - 2{U2){U2U3) + 2(m2)^(u3> 

1 h{t2 - 



2 h{t4 - ti)'- 



-Hh ~ t2)h{t3 ~ uy 



(174) 



This procedure can be pursued to obtain higher p-time distri- 
butions at fixed avalanche duration. We will stop here, and 
just point out that one can check explicitly that (11211) satisfies 



lim dx^dx^Z^^O 



Ai .Ag-^OO 



(175) 
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consistent with the fact that there is only a single avalanche 
to this order, since a non-zero value would require that <2 and 
t4 are in two separate avalanches, since the limit A3 — > —00 
selects U3 = 0. 



D. Interpretation of the instanton solution: response to a small 
step in the force 

Here we examine the question of what is the physical mean- 
ing of the instanton solution u^j? We show that it encodes the 
(linear) response to a small (infinitesimal) step in the applied 
force at x, t, equivalently a small kick in the driving veloc- 
ity. The inverse Laplace transform of u^^. is then related to the 
change in the probability distribution of ii due to this kick. 

First note that the action in presence of the source A, noted 
^tree ([9o] | ig guch that iixt does not fluctuate. This means 
that all cumulants of ii and u involving at least 2 response 
fields vanish. In other words, in any expectation value the 
field Uxt can be replaced by u^^. Hence from Eq 



Ux't' = {Ux't')s\ 



1 /- 



G[X] 



1 (5e-''x* '^"^'^^"'""^'^ 



G[A] SU 



1 Se-I^^t 



G[X] 



gx 



SWx"t' 



(176) 



By definition of the response field, since Ux>t couples to 
fx't' = J^„ 9x']c"''^x"t', see Eqs. ( fT2l) and ( |79] |. it is the re- 
sponse to a change in the driving from Wxt = vt ^ vt + 6wxt, 
and more precisely to an infinitesimal kick Swxt = Sw S{x — 
x') S{t — t') in the velocity at position x' and time t' . Note that 
(1176b is independent of v, a fact which comes from the form 
(I92I 1 and is a peculiarity of the tree theory (at fixed r/ and <t). 

Taking a derivative w.r.t. A at A = 0, and comparing with 
(I94I 1 yields the property that for the tree theory the exact re- 
sponse function TZxt.x't' (in the velocity theory) is uncorrected 
by disorder. 



7^ 



xt,x't' 



[Ux't'Uxt/S^'^" 



Rx 



:= [Ux't'UxtjSo ) 

(177) 

as clearly the cubic vertex (l86b cannot lead to corrections of 
the response. This is in agreement with the fact that the ef- 
fective action T = S for this theory as discussed in detail in 
Isill . Note that Eq. ( 1177b is a non-trivial property for w = 0+, 
since then, in most realizations of the disorder, the particle is 
not moving and under a kick it will experience only a small 
avalanche (of the order of the cutoff). 

Let us now use Eq. ( 1176b in the limit of w — > 0+, i.e. order 
in V, but to lowest order in the perturbation 6f, 



- 1 = / Sfx't'U^x't' + 0{(,5ff) . (178) 



We used that G[A] = 1 for tj = 0+. The instanton solution 
thus gives the statistics of the motion induced by the kick. 



For instance, let us apply Eq. (11781 ) to calculate the center-of- 
mass velocity it 1 = iit^ at time ti, choosing Xxt = X5{t~ti), 
given that there was an infinitesimal uniform kick 5wxt = 
5w5{t — to) at some time to < ti, on top of the f = 0+ 
stationary state. The instanton solution is uniform Uxtg = Uj), 
and precisely encodes that information 



(179) 



Note that Eq. ( |179b can be generalized to any source A(t), 
hence the instanton solution u^^^ gives the first order in 5w of 
the generating function of velocities at any later times; u^^^ 
does not depend on the sources at times smaller than to- 
Performing the inverse Laplace transform of the instanton 
solution w.rt. s := —XL'' gives 



lt: 



ow 



(180) 



This is the linear change of the velocity distribution at time ti 
as response to an infinitesimal kick at time to < ti- Using the 
explicit form for the instanton solution (llOOl ) and performing 
its Laplace inversion we find from ( 11791 ) (restoring all units): 



6{ui) ( 1 
6w e 



Sw 



1 



3(tl-t0)/T™ - 1 



+ - 



+ 0{5w^), (181) 



which is interpreted as follows; For v = 0+, at a given time 
tp , almost surely the particle has zero velocity. The infinites- 
imal kick at time to produces an avalanche (it gives a velocity 
= Sw/Tm) which most of the times dies out well before 
time ti (in a time ^ tq, the microscopic cutoff time). Excep- 
tionally rarely, however, and with probability 0{5w), this kick 
produces a larger avalanche, i.e. lasting a time of order t„i. 
Hence the result that the response function is unchanged by 
disorder is not trivial at all: For most realizations TmSuti I 
is very small; however for some realizations Siit-^ = 0{1) 
hence Siit^/Sw ^ \/6w. After averaging over disorder these 
rare events lead to the bare response function, which is 0(1). 

Let us now comment on stationary versus non-stationary 
avalanches. In previous sections, and most of the paper, we 
study avalanches in the steady state, obtained by time-uniform 
driving Wxt = vt (with small v). These can thus be called 
stationary avalanches. Adding a kick at time to leads to non- 
stationary driving. Indeed the avalanche generated by the kick 
appears non-stationary, i.e. P{ui) in ( 11811 ) is quite different 
from the 1-time distribution found in Eqs (11021 ). ( 1107b . It is 
time (i.e. ti) dependent, and for instance the average veloc- 
ity decays exponentially, iii = iHg-(ti-*o) Qjjg ^.^jj ^gj^ 
whether such non-stationary avalanches are qualitatively dif- 
ferent from the stationary ones. 

For an infinitesimal kick, this is not the case. Indeed, if one 
considers as in Section IllICI to lowest order in v the steady 
state, i.e. the distribution of probability of iii = ?i(ti), condi- 
tioned to an avalanche having started at to, one obtains exactly 
P{ui), as given in Eq. ( 11811 ): As usual, this conditional proba- 
bility is obtained as —dta(l'iPi{ui) using formula ( 1137b (ti, t2 
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there are to, ti here, respectively). This, in fact, is more gen- 
erally true: Namely an infinitesimal uniform kick at time to 
produces the same velocity statistics for t > to as condition- 
ing an avalanche in the steady state to start at time to- It can 
be shown at the mean-field level from the identity 



lim ( 

A— — oo 



-dt. 



+00 



dt u 



A(t) = XS{t - to) + fi{t) 



(182) 
(183) 



Here = for t < to, but ^{t) is arbitrary for t > to- The 
rh.s of ( |182t is related (via Laplace inversion) to the effect 
of the infinitesimal kick at time to on the joint distribution of 
the velocities at all later times, while the l.h.s. is related to the 
velocity distribution conditioned to the avalanche starting at to 
(the conditioning results from the operation — lim>_j._oo dtg 
as we learned in Section lTlICI ). The proof of this result, which 
is easy to obtain from the instanton equation, and more details 
on these properties will be given in Ii83il . 



E. Finite step in the force and arbitrary monotonous driving 

For completeness, let us discuss the case of a finite kick, 
studied in [84]. First one notes that one can generalize 
our method to arbitrary monotonous driving. Starting from 
Eq. (|79] l in the laboratory frame (i.e. setting v ~ 0), but with 
arbitrary driving fxt > 0, we follow the same steps as in Sec- 
tion UnB] to obtain for the generating function of velocities 



(184) 



Here a = — A'(0+). The Middleton theorem allows to re- 
strict the path integral to positive velocities Uxt > 0. Again, 
integrating over ii^t enforces the instanton equation to be sat- 
isfied. Inserting its solution thus eliminates all terms propor- 
tional to ii, such that we are left with llsill 



(185) 



As written, on an unbounded time domain, this formula holds 
if and only if all trajectories are forward for all times. It 
can thus be applied for w = 0+ and an infinitesimal kick 
fxt = 5 fxt > 0, recovering (1178b and (11791 ) by expanding 
to lowest order in 5f (and to order in i;). It also holds for 
any finite kick, and allows to study arbitrary non-stationary 
monotonous driving as done in detail in 118411 . For instance, 
one can prepare the system di t = to in the quasi-static 
Middleton state Ux{w): In the distant past one first drives 
monotonously with fxt > to erase the memory of the initial 
condition, then stops driving. The above formula implies 



(186) 



This can be used to study non-stationary avalanches obtained 
from the Middleton state at t = to, generated by applying 
a finite kick Sf = ra^Sw at time to- Interestingly, these 
avalanches can also be shown, within mean field, to be equiv- 
alent to those of the steady state, under conditioning of the 
velocity a t to to be equal to uto+ ~ 5w as will be dis- 
cussed in llo . Note however that these formulae do not say 
anything about non-monotonous driving as in the hysteresis 
loop, which remains to be investigated. They only pertain to 
avalanches in the Middleton state. 

Consider now an application to a spatially non-uniform 
kick at time to, of arbitrary finite strength fxt = SfxS{t — to)- 
It is interesting to note that any observable involving the cen- 
tor of mass at later times depends only on 5 fx, since the 
associated source Xxt = At is spatially uniform; hence the 
instanton solution is also spatially uniform, u^^ = u^- One 
consequence is that the probability that the avalanche which 
started at to has terminated before ti. 



P{T <ti)^ lim eJ^^^-< 

A— )■ — CX3 



_.At=A5(t-t]_) 



1 



0{5fl] 



(in dimensionless units) also depends only on 5 fx- This is 
because, although an avalanche has ended if and only if all 
Uxt = 0, thanks to Middleton's theorem this is equivalent to 
the center-of-mass velocity being zero. Hence we can use the 
uniform source At = \5{t — ti), leading to the above explicit 
expression, which we use below. 

As a last application, to be discussed again below, consider 
an arbitrary driving fxt > for t > to with the initial con- 
dition ( 1187b . Let us define a kick of finite duration t[ — to as 
a driving such that fxt > for to < t < tf and fxt = 
for t > tf. Consider a source Aa;t = 12^=1 ^ tj) with 
to < ti < ■ ■ ■ < tn- The solution of the instanton equation 
with such a source was studied in Section UlI CI One can 
check that in the limit of all Ai — )• — 00 the instanton solution 
takes a very simple form (in dimensionless units), namely 



Ed{tj^i <t < tj) 
1 — e*J^* 

Hence we obtain the joint probability 

Prob(uti = 0, Ut2 = 0, . . . , Ut„ = 0) 

AtlJxt _ /■*" 



exp 



(188) 



(189) 



We can learn a lot from this formula: First, for n = 1, we 
see that utj can vanish (strictly) only either when the driving 
has stopped strictly before ti, e.g. fxt = for alHo < ^/ < 
t < ti, or if it stops at ti, e.g. fxt ~ (t — ti)" with a > 



with initial condition 



Uxto 



Uxto = Ux{wto) 



(187) The time ordering there was opposite. 
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such that the integral remains finite. Hence a kick of finite 
duration produces only a single avalanche which lasts longer 
than t[ — to, more precisely, taking a derivative w.rt. ti, 



P{T = h) = 



4sinh2(ii^) 



(190) 



Then, for n > 1, formula ( 11891 ) allows to analyze the case of 
a succession of several kicks of finite duration. Because the 
joint probability takes the form of a product on each interval 
[ti,ti^i], it shows that for a given f^t the events u{ti) are 
statistically independent 

To conclude, let us note that the formula (11851 ) being more 
general, it also allows to study the properties of stationary 
avalanches in the steady state with constant driving Wxt = v 
(see e.g. Isill ). However formulae such as ( |189l l and ( 1190b 
do not readily apply (they would lead to divergent integrals). 
This is because one must perform the limit of infinite Laplace 
parameters Xi after integration over time, the physics of the 
single avalanche being restored for v = 0+ as explained in 
details in Section lTlICI 



F. Recovering the quasi-static avalanche-size distribution 

Here we show how to recover the quasi-static avalanche- 
size distribution, first within the stationary state at a constant 
small driving velocity v, by measuring for a finite time, and 
second in a non-stationary setting, by driving the system over 
a finite distance. The results for the avalanche-size distribu- 
tion in a finite time window are new and of experimental in- 
terest. Some results at the end about a finite driving are also 
new. 



1. Steady state: Limit of infinite time window 

Consider the center of mass, i.e. the total size S of an 
avalanche. In the limit of small v, in the comoving frame, 

the latter is S* = L'^ J^^^^dt{v + ut), where T is a time 
much larger than the typical single-avalanche duration, but 
much shorter than the waiting time between two consecutive 
avalanches. We want to compute 



(191) 



One would like to take T oo, and consider a static source 
Xt := X. The instanton equation then admits static solutions. 



The one of interest is 



4Acr 



2cr 



(193) 



The other root is not continuously related to w = at A = 0, 
and for this reason we reject it. The solution (11931 ) has to be 
injected into Eq. ( 1931 ). Due to the time integral in the latter, 
this leads to an infinite Z{X). Hence to recover the avalanche- 
size distribution from the dynamics in the setting of a constant 
driving, w{t) = vt, one must be more careful and consider T 
large, but not infinite. For instance, we may consider a square 
source 



At = X0{t2 ~ t)9{t ti) 



(194) 



with ti = -T/2 and = T/2. If T is large enough, the 
solution is expected to look like 

ut = , t>t2 (195) 
Ut = u{X) , <i < t < i2 (196) 
ut=0 , t<^ti. (197) 

One then finds, expanding ( l92l ) in small v, 



(198) 

We work here in the limit T ^ Tm, but povT ^ 1. On the 
other hand, we know that quasi-static avalanches obey 



^XL-'[u{w)^u(a)] J dSp{S){e^-'^ - l)w + 0{w^) 

= L'^Zs{X)w + 0{w'^) . (199) 
Here we denoted (instead of Z( A) as in Ref. iFtsIi ') 

ZsiX) = L-^ (e^^ - 1), = ^ ( {^'') l) (200) 

the generating function for quasi-static avalanche sizes. 
{■ ■ ■) p denotes the un-normalized average'^ w.rt. p and we 
have used ( fT9] l to transform it into a normalized average 
over P{S). Identifying w = vT and the total displacement 

u{'w) — u{Q) = ut/2 — we obtain 

Zg{X) = iTi^uiX) . (201) 
Hence we recover the tree result for the size distribution (17311 



^(A) = 



1 - yi - 4A5'„ 



(202) 



Ut=U 



'Tn^u + au^ = —A 



(192) 



It leads, upon inverse Laplace transformation, to P{S) given 
by Eqs. (|32] | and (l33T l. Note that the same procedure can be 
performed to recover the local avalanche-size distribution, by 
considering a time independent but space dependent solution 
of the instanton equation. One then recovers, for instance, the 
results obtained in section IX of |73fl. 



There is no contradiction witli tlie fact tliat for a single kick lit^ = 
implies ut2>ti = 0: Indeed, the probability of the second event is one if 
the driving vanishes on the interval 1 1 , t2 • 



Note however that the expression with p also holds for a continuum 
avalanche process with no cutoff. From il9\ it is normalized to the vol- 
ume (S)p = L'*, see jvl. 
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2. Steady state: Distribution of avalanche sizes during a finite 
time window 

To be complete, we now show that the solution of the in- 
stanton equation indeed has the form (11951 ) at large T, i.e. 
that the static fixed point is attractive. This also provides a 
novel physical observable for measurements restricted to a fi- 
nite time window. The effect of finite space windows has been 
studied before in the avalanche context in Isoll . while a gen- 
eral study of windows in scale invariant Gaussian signals can 
be found in S. 

We solve the instanton equation in dimensionless units for 
a square source \t ■= X6{t2 — t)6{t — ti). 



{dt - i)ut + u^t = -Ae(t2 - t)e{t - h) 



Its solution is 



ut = 
1 

"*^2 



ti <t <t2 



t > t2 



(203) 



(204) 



(pxiz) ~ tanh(z), C\ = arctanh( 
0a(^) = coth(z), C\ = arcoth 



V\/l -4A 

( ' 



Ut 



t<ti. 



A < 
A > 



The two branches depending on the sign of A are actually 
identical (by analytic continuation) since tanh(z + z7r/2) = 
coth z. We see on these solutions that the above fixed-point 
form (11951 ) indeed holds at large T. 

We now study the probability distribution of the total dis- 
placement during a time-window size T, i.e. of the observable 



U 



T/2 



dt{v -i- Ut) ■ 



(205) 



-T/2 



This quantity is clearly of experimental interest. (For simplic- 
ity we have suppressed all factors of L'', which can be restored 
at the end). It should interpolate between the distribution of 
the instantaneous velocity at short times, and the distribution 
of sizes of quasi-static avalanches at large times. To check 
this, we compute Z{\) = ut using Eq. ( 1204b . which leads 
to 



Z(A) = 



T + ln(l - 4A) 



In I (1 - 2A) sinh ( -\/l-4A 



T 



T 



+ Vl -4Acosh l^-Vl "4Aj J 

AT+— — + -A2(2A-l)r3 + 0(T^) . (206) 
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In the last line we have indicated the behaviour at small T. 
The series expansion in A, which gives the moments, is also 
instructive. 



Z{\) = \T + \^{T + e- 



l) 



+2[r-2 + e-^(2 + T)]A^ + 0(A'') . (207) 

It shows that U = vT exactly, as expected, and that at large T 
all moments grow linearly as 



UP = V [cpT + dp + 0(r"''e-^)] , 



(208) 



i.e. up to exponentially decaying terms, and with possible 
power-law prefactors. 

As in the preceding section, in the small-w limit the proba- 
bility distribution of U is expected to take the form 



PiU) = (1 - povT)SiU) + PovTViU) 



(209) 



Here povT is the probability that an avalanche has started 
inside the time window T. Note that if U is non-zero, the 
avalanche can have started anytime during the time window 
and may, or may not, have finished during that time. U thus 
contains information about the signal measured in a time win- 
dow without the necessity to determine when the avalanche 
starts or ends. 

Since U ~ povT{U) ~ vT from the above, (where and 
below (...) denotes moments w.rt. the distribution V{U)), us- 
ing Eq. ( fT9] l we obtain the remarkable property that the first 
moment of the distribution P{U), 



(U) = (S) 



lim (U) 

T— !-oo 



(210) 



is independent of T. The distribution V{U) can then be ob- 
tained via Laplace inversion. 



1 



UV{U) = LT^i^dxZiX) 



X=-s 



LT 



-1 



T(2s+ 1) 



4s (sT - 1) 



ViiTT ((2s + 1) tanh (iV4s + IT) + V4s + l) 
4s(T+ 1) + r 



(4s + 1) 



(211) 



For s = 0, this yields Eq. ( I210l l. The Laplace inversion is 
performed in Appendix|A] Here we give some general features 
and limiting behaviors. First note that for any finite T the 
apparent singularity at s = —1/4 is fictitious, since the LT is 
analytic there. The closest singularity is at si{T) < —1/4, 
and the leading exponential decay at large U is proportional 
to e'^iC^)^ where si(T) = -l/T at small T, and si = -1/4 
at large T. 

Examining Eq. (121 lb at large s ^ max(l, 1/T^) shows 
that the small-?/ behaviour at fixed T is independent of T, 
and given for U <C min(l, T^) by 



nu) 



{S) 



20F[/3/2 



(212) 
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The persistence of this strong divergence at small U, which 
requires a short-scale cutoff Uq ^ So, is consistent with the 
property (121 OK since demanding normalization to unity of 
ViU) leads to ([/) - y/U^, i.e. ([/) - VUqS^. in dimen- 
sionfull units. 

At large T one can set the tanh in Eq. (121 II ) to unity and 
obtain 



ViU)^{S)< 



-C//4 



(213) 



1 



1 - erf 



U 



The neglected terms ... give the subdominant exponentially 
decaying part in Eq. ( |208t . while the linear and constant parts 
(i.e. Cp and dp) are reproduced by this formula. It thus gives 
the leading correction to a measurement of the avalanche-size 
distribution if the time window is not large enough. Restoring 
units we find that these corrections are decaying quite slowly 
as 0{t„JT). They do exhibit a divergence - 1/{2TU) at 
small U, but which is too weak to correct the tail [/^^ with 
T = 3/2 which agrees with the distribution (|32] |, ( l33T l. 

We note that the above formulae (in Laplace) are reminis- 
cent, but different, from the ones leading to the joint distribu- 
tion of avalanche durations and sizes given in Il84ll . 



3. Avalanches size distribution in non-stationary driving 

In the first part of this section, we have considered what 
happens when measuring the avalanche-size distribution in the 
steady state obtained by constant driving wt — vt, during a fi- 
nite time. On the other hand, one may also consider what 
happens when the system is driven only over a finite distance 
6w, i.e. in a non-stationary setting. For this we recall the dis- 
cussion of arbitrary monotonous driving in Section IIII El and 
use formula (11861 ). We work in the laboratory frame and focus 
on the case where the system is prepared at rest in the Mid- 
dleton state, as described there and in Ref. llsil . i.e. Wt = Wtg 
for ti < t < to and ti — > — oo. The driving is turned back on 
at to. Hence at t = to one has Uxto = Ux{wtg), zero velocity 
Uxta = 0, and formula (1186b holds for t > to. Since the par- 
ticle has been at rest for a while for t < to we define the total 
avalanche size as 



dtiit = L'^{u+oo - uto) 



(214) 



To compute its distribution we can choose a source X^t = A, 
for t > to, independent of space and time. The advantage 
of this setting is that the instanton solution is then simply the 
constant solution, Uxt — m(A) for t > to, given by Eq. (11931 ). 
Hence one has, denoting wq = wtg '■ 



Sw := j Wf dt = 

ho 



Wo 



(215) 
(216) 



Note that at this stage we consider an arbitrary driving wt > 
for t > to, i.e. we only assume that Sw < oo. We have not 
assumed it to be slow or small. To fix ideas, two extreme 
examples are: 

• A kick Wxt = S'wS{t — to) 

• A constant driving during a finite window, Wt = v for 
to < t < ti and Wt = for t > ti, such that Sw = 
v{ti - to). 

Now we know, from Middleton's theorem, that 



u+oo := lim ut ~ u[wo + Sw) 

f oo 



(217) 



Hence we have found that 



- Sw 



(218) 

with Sm = cr/m'*, for arbitrary Sw. In the limit of small Sw, 
from the definition ( 11991 ) of Za{\) we recover again Za{X) = 
m^u{X). But we find more. By Laplace inversion one obtains 
the distribution of S, 



PsUS) 



L'^Sw 



s_ 

-e 4s™ 



(219) 



This is Eq. (33) of Ref. Il84ll where it was obtained for the kick 
and for a particle {d = 0), but as we see here, it is indepen- 
dent of the precise fonn of the driving, depending only on Sw. 
What is remarkable is that the probability (1219b is two things 
in one: 

(i) It is the distribution of size S = dt ut of the 
avalanche, produced by an arbitrary driving resulting in a total 
shift of the quadratic well of Sw = /^^ Wt- Since the driving 
velocity can be arbitrarily large this is a priori a non-trivial dy- 
namical observable. Note that for the kick one is guaranteed 
that there is a single avalanche, but if wt has a more compli- 
cated form then S may encompass several avalanches, sepa- 
rated by time regions where ii = Q, e.g. for a succession of 
several finite duration kicks, as discussed in Section lTlIEl 

(ii) It is also the distribution of 



S = L'^ [u{wo + Sw) - u{wo)] 

= / Ux{wo + Sw) - Ux{wo) 



(220) 



a quasi-static observable, which for finite Sw may also encom- 
pass several quasi-static avalanches, since e.g. {S) = U^Sw 
In the limit of small L'^'Sw <^ Sm one recovers the form (l32l 
of P{S) for a a single avalanche for S 3> Ssw where Ssw = 
{L'^Sw)'^ / Sm acts as a small-scale cutoff. The true single- 
avalanche limit however is reached only when Ssw ~ <S'o- 



In the limit where the microscopic cutoff So — ^ there are infinitely many 
small avalanches. 
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The fact that (i) and (ii) are the same is a simple, but remark- 
able, consequence of Middleton's theorem. The fact that the 
form for P{S) is given by ( 12191 ), and the property ( |218t , are 
a consequence of the simplified tree theory'^. As discussed 
below, its use is justified for d > d^c, and a priori only in the 
limit of slow driving w ^ Vm- The property (12181 ) is consis- 
tent with u{w) being a Levy process, i.e. a jump process made 
of statistically independent avalanches, each distributed with 
the single avalanche distribution P{S) from ( |32] |. The prop- 
erty recovered here is also present in the statics, i.e. for the 
process u^*''^^{w), in mean field, in the BFM and in the Burg- 
ers equation. It has the same P{S), as is discussed in detail in 

Q. 

Finally, a similar analysis can be performed for the proba- 
bility distribution of the local observable 



S't' 



(l)x[ux{wo + Sw) - Ux{wo)] 



(221) 



One must then solve the space-dependent instanton equation 
with a source Xxt ~ <t>{^)^ which is a hard problem. In 
the case (t){x) ~ 5{xi), i.e. a hyperplane in a rf-dimensional 
space, the time-independent instanton solution is known, see 
section IX of ItI: 



(222) 
(223) 



A(Z) = ^Z(Z-6)(Z-12). 



The Laplace inversion is involved, but a simple generalization 
of Eq. (220) in fTSl] . The same trick yields the (normalized) 
probability distribution of S'^ ~ S, 



P. 

Pwis) 
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(224) 
(225) 



Here (S"^) = wL'^ and for L'^Sw < Sm, i-e. w < 1, one 
recovers pm(s) wp{s) with p(s) = 2Ki{-^)/{tts), the 
(rescaled) single-avalanche size distribution obtained in iItsIi . 
These results are exact for the BFM (discussed below), an ap- 
plication being a single-site avalanche for a string (d = 1). 



Hence we first discuss these parameters and their universality. 
In a second stage, the tree theory is identified with the BFM 
and the ABBM model is recovered. 

We recall that the upper critical dimension is due = 27 for 
an arbitrary elastic kernel behaving as e{q) ~ q"^, i.e. d^c = 4 
for usual SR elasticity (7 = 2) and due = 2 for the most 
common LR elasticity (7 = 1). 



1. Improved tree theory and the parameters of the model 

We have shown above that to lowest order in perturbation 
theory in the bare disorder, all generating functions of the ve- 
locity, to first order in v, are given by the sum of tree graphs. 
Equivalently, they can be computed from the simplified tree 
action S'*'^'^'' defined in Eq. dSTl i. At the bare level, this action 
only contains three parameters 770, m and ctq = — Aq(0+). 
These bare parameters are corrected by disorder, and acquire 
a dependence on m, as we now discuss. 

Let us now use well-established results from the ERG ap- 
proach to the statics and dynamics of elastic interfaces. First, 
TO is uncorrected to any order in perturbation theory thanks 
to the STS symmetry, hence we can use everywhere the bare 
mass m. Second, perturbation theory converges for d > due 
(in a sense recalled in AppendixlBl). Third, at c? = due there are 
only two operators which become marginally relevant. The 
first one is the local part of the renormalized disorder, A(u), 
which actually is a function of u; so in principle there is an 
infinity of marginally relevant directions. However, as far as 
single avalanches are concerned, we only need A'(0"''): It is 
shown in Section |IV] that the higher derivatives lead to loop 
corrections, i.e. are important only for d < due- The second 
parameter is the renormalized friction 77. Both parameters, 
A'(0''') and 77, receive logarithmically divergent corrections 
in d = due from 1-loop diagrams. These are cut off by the 
mass m and can be resummed using the FRG flow equation to 
1-loop order 

Let us now determine the renormaUzed parameters at the 
upper critical dimension d = due- Define £ := ln(A/?7i), 
where A is a small-scale UV cutoff; at d = due, for SR elas- 
ticity, set 



G. Mean-field theory for avalanches: The Brownian-force 
model and its ABMM limit 

We are now ready to discuss the correct mean-field theory 
for the avalanche motion of elastic interfaces in the limit v ^• 
0+, and to identify its universal properties in the limit of small 

TO. 

In a nutshell, the mean-field theory is the tree theory, with 
however a renormalization of two parameters of the model. 



A(u) = 8tt'^A{u) = 8Tr^A{ur'^')e~^+^^' . (226) 

Then the FRG flow equation for A(u) is (B.14) in ill. As 
m 0, the rescaled correlator tends to a fixed point A(m) 
A*(u), which is the same one obtained to first order in a e- 
expansion i.e. A*(u) = lime_>o A*(u)/e. Similarly, see e.g. 
114811 . one obtains 

9, In 77 = -A"(0+) = -A"{0+)r^ . (227) 



The full free theory with an arbitrary A(«) does not satisfy property | |2181 . 



Hence, the two parameters of the model acquire a universal 
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dependence on to, in the limit of to — ^ 

77^?7,„~?7o[1ii(A/to)]'' (228) 

Both zi and Ci are defined by the 1 -loop result for the dynamic 
and roughness exponents, 

z = 2- A*"(0+) = 2 + zie + 0(e2) (229) 
C = Cie + 0(e2) (230) 



done is to replace the original equation of motion (|79] l in a 
disorder described by the gaussian force correlator Ao(u) by 
a disorder described by a (renormalized) correlator A(u) = 
A(0) + A'(0+)|m|, since we have neglected all higher-order 
derivatives A*^") (0''"); the latter become important only upon 
considering loop corrections to the velocity distributions. This 
means that this (simplified) tree theory describes exactly an 
elastic manifold in a Brownian force landscape F{x,u) with 
Gaussian correlations, 



(Ci-l)/3 = -2/9 for non-periodic F{x,u)F{x' ,u') = (5''(x - a;') A(0) - ct| 



with Qi ~ 1/3 and zi 
SR disorder 

The above formulae extend to LR elasticity by changing 
everywhere above m — defined below in (14421 ). and the 
factor Stt^ — > Cd=d..,.'v (see its definition and detailed dis- 
cussion in section X of IITSll ) with 6*2,1 = 27r, the fixed point 
A*(u) being unchanged. 

We can now make a precise statement, based on the effec- 
tive action V of the theory. For its definition see ll97ll. and 
in the context of FRG e.g. iHJtI 98] (statics) and HI] 
(dynamics), summarized in llSlfl . Appendix A. It is a gen- 
eral property of F that all connected correlations of the theory 
(here of the velocity field) are tree diagrams in F: The vertices 
of the trees are vertices not of the original action S, but ver- 
tices of F, i.e. renormalized vertices, which contain all loop 
diagrams. 

When d due and in the limit of to 0, the effective ac- 
tion F becomes simpler and its limit is the so-called improved 
action. This is discussed in Appendix|Bj where we show how 
the irrelevant operators become negligible for d ~ due, when 
properly scaled. For instance, the higher time derivatives in 
the equation of motion, or higher disorder cumulants, become 
negligible, and one can focus on ?/ and A(m) only. 

If in addition one considers positive driving only, fxt > 0, 
then for d = due the effective action of the velocity theory is 



(231) 



5, 



tree I r;, A' (0+)' 



I.e. the tree action with the renormalized 
parameters a — > (j,n and 77 ?y,„ . It sums tree graphs except 
for the renormalization of a and rj, which contain loop cor- 
rections. This remains true for d > due, where cr and jy flow 
to non-universal limits as to — > 0, as discussed in Appendix 
IB] Note that the statement we make here is only for v = 0+: 
Since we have not analyzed the FRG flow at non-zero v, we 
focus on the limit of small driving. This also means a small 
step in the force, i.e. a small kick, in the non-stationary setting 
discussed in Section lTlIEl 

For d < due the behavior is universal but different from 
mean-field, and is analyzed in SectionlTVl 



2. Brownian force model (BFM) or elastically coupled ABB M 
models and universality 

The mean-field tree-level theory has a very simple inter- 
pretation. It is clear from Section IIIIBI that what has been 



Since Sm = dm/'m'^, this corrects a misprint in Eq. (108) of 17311 . 



where a = — A'(0"''). Such a landscape is constructed in a 
spatially discretized version, by considering that for each x, 
F{x, u) performs a Brownian motion (BM) as a function of 
u, and that these BMs are mutually independent for different 
X. Furthermore, they are stationary Brownian motions, hence 
they are constructed by considering e.g. a much larger peri- 
odic system in the u direction. An elastic manifold of internal 
dimension d in such a landscape is called the Brownian force 
model (BFM) ll74ll . The statics of this model was studied in 
ll74ll . As we discuss below, a non-stationary BM version can 
also be considered. 

Hence, from the previous paragraph we conclude that the 
full statistics of the velocity field in an avalanche for an inter- 
face at d > due identifies in the small-TO, small-z; Hmit with 
that of the BFM, with parameters a (t,„, r] ^ rj^- This 
BFM can also be described as a set of ABBM models for each 
Uxt with an elastic coupling g^^, between them. 

A crucial property of the BFM is that the dynamics of the 
center of mass of the elastic manifold is described by the 
ABBM model i.e. by equations Q and Q. Intuitively 

it is easy to understand why: To compute center-of-mass ob- 
servables in perturbation theory we need to consider all graphs 
with external momenta set to zero, q = 0. However, since we 
have summed only tree graphs, it implies that all propagators 
are evaluated at q ~ 0. Hence, apart from the (non-trivial) 
renormalization of the parameters of the model, in effect, the 
avalanche dynamics of the center of mass iit for w = 0+ is 
described by the ABBM model, i.e. a single point driven in 
a long-range correlated random-force landscape, F{u), with 
Brownian statistics. It amounts to suppressing the space de- 
pendence in Eq. ( |79l ), hence corresponds in our general model 
to the special case d = and Ao(0) — Ao(u) — a\u\. 

Let us now connect our previous results, obtained directly 
for the center of mass of the interface, to the standard analy- 
sis of the ABBM model. Then we will revisit the BFM, and 
finally calculate observables beyond the center of mass, re- 
quiring the full power of the BFM. 



3. Center-of-mass obsen'ables and ABBM model 

Let US recall the original solution lH, H] of the ABBM 
model, based on a Fokker-Planck approach (see more details 
in ifsill). The equation of motion (|79] l for the instantaneous ve- 
locity in the laboratory frame v = ut of a particle in a Brown- 
ian landscape (suppressing internal degrees of freedom x) can 
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be written as a stochastic equation 

rjdv = m^(w — v)dt + dF 



(232) 



where dF^ — 2a\/dt . The associated Fokker-Planck equa- 
tion for the probabiHty distribution Q = Q(v, t|vi, 0) of the 
velocity at time t, given velocity vi at time t = is 



Tjdt 



-d,ivQ)+m^{v-v)Q 



7] 



(233) 



It satisfies Q(v2, 0+|vi, 0) = (5(v2 — Vi). It is normalized to 
unity at all times upon integration over the final velocity v, 
thus it is the propagator of the system. For w > 0, it evolves at 
large times to the stationary (zero current) distribution Qq := 
limt_j.oo Q with 



Qo(v) 



- - 



V /Vrr 



(234) 



Here Vm = Sm/T„i, Sm = a/m^ and r,„ = rj/m^. Note 
that here we study a point particle, hence the velocity scale is 
f„i; if we study the center of mass of an interface, it is to be 
replaced by Vm as discussed in Section lTlIG II 

One notes that taking w — > 0+ and forgetting the normal- 
ization, Qq converges to the single-time velocity distribution 
obtained above in Eq. ( 1107b by a completely different method. 
There, the normalization was fixed from considerations of a 
small-scale cutoff. Similarly, in the limit f 0+, one finds 
that the propagator takes the form 



Q(v,t|vi,0) 



1 ^ 



V t 



Vl 







(235) 



with 



Q{v2,t\vi,Q) = Vie"" 



P2{VI,V2) H 1 



Vl 



~5{V2) 



(236) 

The term ^2(^1, ^^2)7 given by Eq. ( 1132b . is indeed a solution 
of ( |233] | with (5(v2,0+|vi,0) = (5(v2 - vi). We note that the 
piece ~ 5{v2), which corresponds to avalanches which have 
already terminated at time t, is necessary for Q to conserve 
probability, i.e. such that d?;2Q(w2, ^l^i, 0) = 1 for all t. 
Since Q is a conditional probability, we can also consider the 
joint distribution of velocities. 



(3(w2,i|wi,0)pi(wi) = Q{v2,t\vi,Q) 



1 

— ( 

Vl 



(237) 



We find that it reproduces the 2-time probabilities given in 
Eqs. (11331 ) and ( 11371 ). More details about the ABBM propa- 
gator and how it behaves in the w — > 0+ limit can be found in 
AppendixiDl 

By using the dynamical field theory of interfaces, we have 
in this paper obtained a novel, and completely independent 
way to solve the ABBM model. Indeed, our method is even 
more powerful, since it allows to treat interfaces and spa- 
tial degrees of freedom, and it can be extended beyond the 



tree level, as will be discussed in the following sections. Al- 
ready its consequences for the ABBM model itself are quite 
interesting: By allowing to compute directly Laplace trans- 
forms through the instanton equation ( |99l ). it provides a useful 
complementary method to the Fokker-Planck approach. For 
avalanche observables it is quite efficient, as was shown in the 
previous sections and Ref. ||84|l . For other observables (such 
as [/ = /^^^2 '^^ non-locality in time makes it very hard 
to obtain the result via the Fokker-Planck method. On the 
other hand, one advantage of the Fokker-Planck approach is 
that since v(t) is a Markov process, the n-time velocity prob- 
ability can be written in a factorized form as 



^ n— 1 

<l'l...n'P(.V,i,...,Un) = —£-"1 Q(Mj + itj + i|Wjij) , 

^1 , 
J = l 

(238) 

where „ is the probability that all n times belong to an 
avalanche. Curiously, it is not easy to recover that property 
immediately from our general expression for Z„ . In Appendix 
|E]we check it exphcitly for n = 3. 

Let us note that since the ABBM model is the zero- 
dimensional limit of the equation of motion (|79b of an inter- 
face, the dynamical-action method can be applied. Hence we 
just found that, for the ABBM model at t; = O"*", the tree 
approximation is exact. In the field theory for the velocity 
it means that the effective action F equals the bare action 5", 
and there are no loop corrections. Hence A' (it) = Ag(u) = 
— (7 sgn(u) is an exact ERG fixed point with scaling exponent 

= 4 — d, as already noted in the statics in iItsIi . Crucial for 
this remarkable property is that the force landscape is a Brow- 
nian, and even in = 0, this is not valid for any other, e.g. 
shorter-ranged, force landscape. These properties and a direct 
solution of the ABBM model at any v are discussed in llsill . 

A word of caution should be said about the notion of the 
duration of an avalanche. In the present tree-level mean-field 
theory (and similarly in the ABBM model) avalanche dura- 
tions can be defined unambiguously for a continuum version 
where the small scale cutoff Sq — > 0, and accordingly the 
avalanche density po 00, as the velocity ii exactly vanishes 
at some time for tj = 0+. In that version there is an infi- 
nite number of infinitely small avalanches and the quasi-static 
process is infinitely divisible (a Levy process) as discussed at 
the end of Section IIIIF3I On the other hand, if one studies 
the original interface model ([T]) with smooth and short-ranged 
disorder, in the limit w = 0+ or in the limit of a small step 
in the force Sw, an avalanche has, strictly, an infinite dura- 
tion (diverging with some power of l/w or 1/ 5w). Indeed the 
starting point is a metastable state (zero force state) with one 
marginally unstable direction and the final state is generically 
a stable zero force state. Near both points the motion is very 
slow, so the duration is very large, but the associated displace- 
ment is negligible. One must thus focus on the part of the 
avalanche motion such that ii^ vq, or such that the interface 
has significantly moved by more than 5*0. This part of the mo- 
tion is universal and described by the ABBM model. It would 
be interesting to make this statement mathematically precise. 



27 



4. ABBM model: Connection between the instanton equation and 
the Fokker-Planck equation 

The Fokker Planck equation can be Laplace-transformed in 
A, or equivalently one can write the evolution equation (in the 
laboratory frame) for 

/■OO 

G{X,t):='^= dve^^P(v,t). (239) 
Jo 

Without specifying the initial conditions, the evolution equa- 
tion is 



The solution can be found in the form 



(240) 



(241) 



with Z{0, t) ^0 since G(A = 0, i) = 1. Then Z satisfies the 
equation 



9^ + 9Z 



(242) 



This equation admits a time-independent solution Z{X) = 
Z{\t) 



Z{X) = - X) 



(243) 



Hence we recover the result (IIOII) obtained via the MSR 
dynamical-action method. 

The connection to the instanton equation can be made as 
follows. The equation ( 1240b can be solved by the method of 
characteristics: Define a function X{t) which obeys the fol- 
lowing differential equation, 



AX{t) 
dt 



= X{t)-X\t). 



(244) 



Further define Z[t) := Z{X{t),t). Then, using Eq. ( |242] |. the 
total derivative is 



dZjt) 
dt 



= m- 



(245) 



Equation ( 1244b is exactly the instanton equation ( |99] l. if one 
identifies X{t) = u{t). For i < 0, it admits the solution 



A(t) 



An 



Ao + (l-Ao)e- 



(246) 



with boundary conditions X{—oo) — 0, and A(0) = Aq. In 
addition 



Z(t) 



A(t') dt' 



(247) 



Hence if we express .^(Aq) := Z{t = 0) as a function of 
Aq = A(0) we obtain precisely (12431 ). 



Eq. (1242b is solved for any initial condition Z{X, t = 0) 
ZoiX) as 



Z{X, t) = - ln(l - A + Ae"*) + Zq 



A 



A + (l- A)e' 



(248) 

Note that the argument of Zq is A(— 1)| Ao-s-a- Hence from 
Eq. (1239] ) we get 



G{X,i) = (l-A + Ae~*)-^Go 



A 



A + (1-A)e* 



Go(A) =e"^«(^) . 
This gives the decay to the steady state as 



u{t) = v{l - er"-) + e"'u(0) 



u{tY =w(l-e-*)2+2u(0)e-*(l-e-*) 



(249) 



(250) 



(251) 



(252) 



It is in agreement with the results of 118411 for a quench in the 
driving velocity. Note that for any t > Q ( |249b behaves as 
G(A,t) ^ A(t){-XY with A{t) = (1 - e-')--Go{-^). 
hence P(u, t) ^ A{t)u"~^ /T{v) and the current at the origin 
vanishes. 



5. Back to the Brownian force model 

Having recalled the properties of the ABBM model, which 
contains the information about the center of mass, we now 
reexamine the BFM which contains all spatial information. 

The Langevin equation (1232b for the ABBM model can be 
rewritten as 



(253) 



■qdtUt = V UtS,{t) + m {wt~ ut) , 



with ^{t)^{t') = 2a6{t — t') a Gaussian white noise. It de- 
scribes the original model (|5]) only if wt > 0. 

Similarly, the BFM can be defined focusing on the evolu- 
tion of the velocity, by the following Langevin equation in the 
laboratory frame: 



■qdtUxt = V Uxti{x, t) + fxt + (V;; - m )uxt , (254) 



W\thi{x,t)^{x' ,t') = 2cr5(i-i')(5(x-a;') uncorrected Gaus- 
sian white noises, with obvious generalization to an arbitrary 
elastic kernel g~^, . It does describe the motion in a stationary 
Brownian random-force landscape if (and only if) driving is 
monotonous f^t > for all times. However, from the discus- 
sion in Section ( IIIIEI) and in llsill . if one complements it with 
an initial condition 







(255) 



it does also describes the motion in the non-stationary Brow- 
nian random-force landscape 



F{x^ u)F{x' , u') = 2(7 min(w, u')6{x — x') 



(256) 
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for u,u' > with initial condition Uxt=o = 0. This setting 
has advantages since the landscape is defined by uncorrelated 
BMs which all start as F(0, x) = 0. This avoids the construc- 
tion of stationary BMs in a large box, as a limiting process. 
The catch is that in position theory it does not satisfy STS; 
this is seen on observables such as ds iixs, whose averages 
are time dependent. If one adds a large box, then these con- 
verge to the stationary BFM observables. 

If one focuses only on velocity observables, and forgets 
about the position theory, the BFM is uniquely defined by 
Eq. ( |2541 i. If one drives for some time, the memory of the 
initial joint distribution of velocities P[{uj.t=o}] is lost. In 
the case of a steady drive, fr^t = v, the system evolves to- 
wards a time-translationally invariant steady state, e.g. for the 
1 -time distribution, 

steady [{wj] , (257) 

which generaUzes (12341 ) and is more complicated to calculate; 
(it requires solving the instanton equation with a space depen- 
dent source). This steady-state measure for the full velocity- 
field identifies with the one of the elastic manifold in dimen- 
sion d > due, and for small v and to, as discussed in the 
previous sections. 

In addition, the BFM is an interesting model to study by 
itself. It can be solved in arbitrary space dimension d, and for 
arbitrary driving, from the general formula: 

^ g/^, i^J^t _ (258) 

which assumes (monotonous) driving from the far past, or for- 
mula ( 11861 ) for the initial condition (12551 ). More details can be 
found in 184(1 . including a formula for an arbitrary initial ve- 
locity distribution. 



H. Spatial fluctuations 

We can now use the full power of the tree theory, i.e. 
the BFM, and calculate space-dependent observables within 
mean-field theory. The space-dependent instanton equation 
allows to go beyond the ABBM model, which describes only 
the center of mass, and to compute spatial fluctuations. 

In addition, the results below are exact for the BFM in 
any space dimension d. Most results concern the BFM 
in the steady state, i.e. they are time-translational invariant, 
as discussed in the previous section. Time-dependent non- 
stationary generalizations are left for the future. 



1. General considerations 

Let us write for completeness the instanton equation for an 
arbitrary elastic kernel g~,j, 




ivdtSxx' - 9xx')^'^'t + f^^^t + ^xt^^ , (259) 



with a = — A'(0"'"). Below, we first perform our calculations 
using the local elasticity 

9xx' = + m^)Sxx' . (260) 

At the end we indicate how the formulae generalize to an ar- 
bitrary elastic kernel. 

Time-independent, but space-dependent solutions of the in- 
stanton equation with a (5-function source in space were stud- 
ied in Refs. iFtsL FtsI [86ll . and allowed to obtain the distribu- 
tion of local avalanche sizes. Finding solutions which are 
both time- and space-dependent is notably more difficult''' and 
must be left for future research. Here we analyze solutions 
which are "almost" space independent, i.e we choose 

Kt = At + Hxt , (261) 

where the spatially dependent part Hxt is small. The solution 
of ( 12591 ) can then be obtained in an expansion in powers of /i. 
We write here the two lowest orders: 

Uxt ^ u1 + + + 0{^i^) (262) 

(,^dt - m^)u°t + <7{u°f ^ -At (263) 

{r]dt + Vl-w? + 2(7u?)uit = -^lxt (264) 

(7y9t + V2 - to2 + 2nu'l)ul, + a{ul,Y = . (265) 

The solutions of Eq. ( 1263b have been discussed in section 
unci Since no general solution for all At exists, let us pro- 
ceed with a solution of Eqs. ( 1264b and ( 1265b . supposing we 
know {t( : 

wit = / Mx't' '^x't'.xt , (266) 

Jx' .t' 

ult = cT [ {ul,t,yRx't',xt . (267) 

We have introduced the dressed response kernel M.x't',xt, 
which will be a fundamental object in the remainder of this 
article. It is solution of the equation 

- +rn'^ 2(7u0] M,,t',.t = 5''{x - x')5{t - t') . 

(268) 

Note that since the instanton equation has the time-derivative 
reversed, we have reversed the order of the arguments in R, so 
that, as defined, it has the usual causal structure of a response 
function. Thus as noted in Eqs. ( 12661 ) and ( 12671 ) it "acts from 
the right", in contrast to the usual convention. 

It is easy to express R in Fourier space, i.e. M.x't',xt — 
/fc]Rfc.t',te*''('='^'=' with 

Rfc t, t, = ^(^^ _ 

V 

(269) 



Note the resemblance of the instanton equation with the KPP-Fisher equa- 
tion for front propagation. 
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First, this allows us to obtain the avalanche statistics in the 
small-velocity stationary state, working, as in Section (IIIICI ). 
to first order in i;. Integrating Eq. (12621 ) over space and time, 
we find 



Jxt 



2 r -d 



Z' = m^L 
Z^ 



^q=0,t',t I IJ^x't 



(270) 
(271) 
(272) 



t,t' 



<i=o,t',t / (uU,)^ . (273) 



To this order in fi we thus obtain averages of the velocity field 
containing two space-dependent velocities. Indeed, to first or- 
der in V, in the small-velocity limit, one can write^*^ 



= 2vm'^a 



^Xiti,x't'^X2t2,x't'^q=0,t' ,t ■ (274) 



(The factor of 2 comes from the fact that u^^ — 0(^^).) This 
is easier^' expressed in Fourier space^^. 



2vm, a 



= 2vm a 



^,ti,t'J»-9,t2,t'^g=0,t',t 




(275) 



Note that we have introduced a graphical notation that will be 
useful later, when calculating loop corrections. 
A nice feature is that the source, which couples to the center 
of mass, is still quite general. For A = 0, R reduces to the 
usual response function and we recover after integration over 
t and f 



1 



-^{rn^W)\tl-t2\ 



(276) 



Here and below we drop the factor v in the term v + iitm the exponential 
since it is subdominant at small v. 

Note that the first-order derivative does not yield any new infomiation: It 
confirms iixiti + v = v for A = 0. 

We use, for an arbitrary function A the short-hand notation Uqti ^q't2 ^ — 
(2TT)'^5'^{q + q')uqtj^U—qt2A for translationally invaiiant correla- 
tions. Hence Ma;itiMa:2t2^ = Xj e"'''"i~^^''"qtiM-9t2^ and, in- 
tegrating over x\,X2, one obtains the center of mass litjiiijA = 
L~'^UqtiU—qt2A\q—o, hcncc rccovcring Eq. l |54t . Everywhere = 

r d<^q 
J (27r)<i ■ 




FIG. 6: Disorer-averaged unfolding of an avalanche following a local 
step in the force at a:o = 4 and of duration T = 1, according to for- 
mula l |28U . Mean intermediate positions u{x,t = 0) + J^{uxt)oT 
are shown at t multiples of T/IO. The velocity Uxt > 0, so the 
motion is towards the top of the plot. For the sake of illustration we 
chose a random but fixed initial condition ^(a::, t = 0). 



This is a finite-momentum generalization of Eq. (|54] l; a factor 
of 2m^ has canceled. 

Next, using the results of Section IIIIDI we also obtain in- 
formation about avalanches following a small local step in the 
applied force at time to, i.e. Sf^t — Sfxd{t — to), 



u^^t'Sfx + OiSf), 



(277) 



where Sf^ ~ J^, 9xx'^^x', where we also work for v 0+, 
but due to the step the leading result is non- vanishing, though 
of order zero in v. 



2. Dressed response function, and space-dependent shape 
following a local force step 

Let us now apply our formulae to the case At = XS{t — t^), 
and pursue in dimensionless units. We recall the instanton 
solution 



A 



Ut = 



A + (1-A)e*3 



Z7^(i3 - 1) 



(278) 



It leads to the dressed response function for a single-time 6- 
function source 



1*3 



1 + A(e*''-*-^ - 1)0(^3 - U) 



l + A(e'»-*3-l)0(t3-t,) 



(279) 



Using Eq. ( 12771 ) to first order in fi, and Eq. ( 1266b . we obtain 
the (linear) response to a local step of the driving force at time 



^t2 



»*2 

^Xi,ti,Xo,to 



^fxo 



(280) 
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Taking A — )• ~oo we obtain the average local avalanche shape 
(i.e. the average velocity conditioned s.t. the avalanche starts 
at to and ends at t2) as 



02 



-Pduration(i2 — ^o) 5 fx 

4sinh(^2^)sinh(^i^) 



sinh(^2^) 



(282) 



(wi)q2 is the center-of-mass shape given in Eq. ( 11561 ). Thus 
the avalanche velocity spreads on average diffusively (for d = 
due) from the seed, i.e. the point where the kick was applied. 
It looks even simpler in Fourier space^^ 



^9*1/02 



^1, 



02 



Rq,ti-t„Sfk 
Rq=Q,ti~to^fq=0 

e"'^'(*^-*°)T§^(^^i)o2 (283) 

OJq=0 



On figure|6]we have drawn the mean advance of an avalanche 
following a local step in the force. 



3. 3-tiine, 2-space point correlation 

Let us now compute the 3-time correlation, in the steady 
state to lowest order in v, using the single-time source at 
and Eq. ( |279t . The t-integral in Eq. (12751 ) is easily performed, 
assuming that > t', 



1 + A(e-I*'-*^I -1) 



t<t' 



q=0,t',t 



1-A 



For ti < t2 < ts, the second integral over t' leads to 



UqtiU-qt^e 



XL''ut„ 



(284) 



(285) 



[AO 



1) + ly [A(e 



(A-l)4(q2 + l) 
X { 3, 2 (q' + 1) ; 2q2 + 3 



A- 1 



By analogy with the procedure used in Section IIII C 61 page 
[TtI taking now the limit A — > — 00 allows to select the contri- 
bution vq'-^ i2S{u3)'Pi2,3{ui,q,U2,-q) the 3-times joint dis- 
tribution P{ui^q, U2,-q, U3). The normalization uqg 12 should 
be the same as for zero momentum q = 0, since if a piece of 
the manifold is moving, the center of mass is also moving. It is 
equal to the probability that the avalanche starts before ti and 



The center-of-mass velocity iit and the velocity of the zero mode (q = 0) 
iiot are in our conventions related via L'^iit = uot- 



ends at t^. As in Section ITlIC 61 we determine it as q'^ ^2 — 

/p°°dsi /g°°ds2limA3^-oo<9Ai9A2^3|Ai=-si,A2 = -S2 = 

ln(z3i/2;2i), and we check that 9*3^3^12 = l/(e*2~*i ~ 1) = 
It^-ti dTPduration(T) Can be obtained also from the duration 
distribution. We can thus take dt^ limA^-00 of Eq. ( 12851 ) to 
obtain the conditional average 



^qti ^ — qt2 



)3 = (5*3 93,12 ) ^ '9*3 , 1™ UqtiU 



qti ^—qt2 



2 (e*2 „ g*3) g9^tl-(9'+l)*2-4t3 



2^tl+t3 I „2^t2+t3 



q e 



+ q'e'''+'' - {q' - l)e 



*l+*2 



X 2F1 (3, 2(g2 + 1); 2g2 + 3; eti-*3) | . (286) 

It is conditioned, s.t. the avalanche started before ti and ended 

at t3. For q = itreduces to {uMiUQt^)^ = 2(1 — e*^^*-^)(l — 
gti-*3) We can obtain the large- g asymptotics using the for- 
mula 



2F1 {a,b + x,c + X, z) 



= {i-zr^ 1 

This yields 



a(c — b) Z 



O ix-^) 



(287) 



(e*3 _ e*2) e9^(*i~*2)-t2-*3 
q2 (-gta _ g*i-) 

X (2e*^+*^ - e*i+*^ - e*i+*3) (288) 



Fixing ti and <3, the function (12861 ) decays monotonically to 
zero for ^2 ts- Depending on the value of q, it is either 
concave (small q) or convex (large q). 



4. 4-time, 2-space point velocity correlations and asymmetry ratio 

To compute the average at a given q, conditioned to both a 
starting time to and a final time for the avalanche, we need 
the more general observable, for to < ti < t2 < t^. 



(289) 



The calculation is more compUcated and done in appendix |F| 
by considering a source At = XoS{t — to) + X3S{t — t^) and its 
associated dressed response function. The full result for (12891 ) 
is displayed in Eq. ( IF8I ). An interesting observation is that at 
q ^ Oil is not invariant under time reversal, i.e. the simultane- 
ous changes to — > — <3, ti — >■ — ^2, ^2 — >■ —ti, ts —to^ and 
Ao ^ A3. This invariance is recovered only at g = 0. Hence 
at the level of the tree theory there is no way to tell the arrow 
of time by watching the center-of-mass velocity, but there is 
an arrow of time for modes with non-zero q. This can already 
be seen on the 4-time velocity correlation function obtained 
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FIG. 7: Plot of the asymmetry ratio A defined in equation ( |291| l. The 
different curves are for g'^ = (solid gray), = 0.2 (dotted red), 
— 1.703 (dashed blue), and q^ — 10 (dot-dahed, green). The 
maximum of A<Atx = -T/2 is attained for = 1.703 (dashed 
blue) The plot is for T = 1. 



started at to and ended at ^3. The times are chosen ordered as 
to < ii < ^2 < ^3. We compute it both for an avalanche 
(0 generated by a uniform small force step at time t = to, i.e. 
Sfxt = Sf6{t — to) with Sf ~ rn^Sw; (ii) for an avalanche in 
the stationary state to first order in v. The two protocols give 
the same result, as was explained in Section ITlIDI it is based 
on the identity ( 11821 ). We present the calculation of (/) in the 
main text; (//) is more involved and is presented in Appendix 

E 

Let us consider the following velocity average following a 
uniform force step at time to. 



U-XitiUx2t2C *^ - / 7 J Of 

Jxo 0f^Xiti0fix2t2 

= 2Swm^a [ R*,% x't'K\ x'fK-of- 

I Xiti^X t X2t2iX^ Q — ^1*1 

Jx'tt' 



from the expression ( 12891 ) by applying L ^'^9ao dx^ \ a,, = A3 =0 ■ 
The general result ( |F9t is bulky, so let us display it here for 
h = tz: 



r2d- ■■ ■■ ^ 2(2g^ + 3)e-^ 



2q^e 



2 -2(g" + l)ti-(g"+2)T 



(l + g2)(2 + g2)(i + 2<72)(3 + 2g2) 



14e 



-T 



(290) 



which is clearly not symmetric under ti — ti, although 
it is for (J = 0. Note that here we do not know when the 
avalanche starts and ends, we only know that the duration is 
larger than T. We define the asymmetry ratio of the 4-time 
velocity correlation as 



(291) 



U-T/2UqfiU-qflUT/2 

It is plotted on figure |7] Since the asymmetry ratio is larger 
at large q at the beginning of the avalanche, it implies that 
the local velocities in an avalanche are higher in the begin- 
ning of an avalanche than at the end. Stated differently, the 
avalanches are more compact at the beginning and the parts 
which move move more quickly. This is consistent with our 
physical intuition that an avalanche starts at some seed Xsccd, 
grows quickly around that point, while at the end it is spatially 
extended, but stops more uniformly. 



We have worked to linear order, i.e. up to terms or order 
0((5/2, w). In Fourier space, and dimensionless units, the lat- 
ter reads 



^qt\ U—qt2^ 



L''X 



= 2Sw I dt' 

to 



q,ti,t' q,t2,t' q=0,t',to ■ 

(292) 

The function R is given in Eq. ( 12791 ). and as written, we can 
drop the 6'-functions. Taking the limit A —00 enforces the 
center-of-mass velocity at the final time to be = 0, leading 
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"Second" shape of an avalanche at non-zero momentum 



We now obtain 



'■qti 



-9*2)0 3 ^-^^ shape fluctuation, or 



FIG. 8: Plot of the conditional average (uqnU-qti) given in 
Eq. l |294| ( for an avalanche starting at time —T/2, and ending at time 
T/2, in our dimensionless units. The different curves are for = 
(solid gray), — 0.5 (dotted red), q^ =2 (dashed blue), and g^ = 9 



second shape, at non-zero wave vector for an avalanche which (dot-dashed, green). The plot is for T ■- 
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FIG. 9: Plot of the (normalized) second shape S2(ti), i.e. the ra- 
tio of conditional second moments of the local velocity S2(ti) 



Jq 



Uq,tiU-q,ti 



f e- 

J n 



UqfiU-qfi 



for an avalanche 



starting at time —T/2, and ending at r/2, normalized s.t. S2 (0) ~ 1. 
The different curves are for a — >■ oo (solid gray, equivalent to the 
same-colored curve on figure |7), a = 1 (blue dashed), a = 0.1 
(green dot-dashed), and a — > (dotted, red), which approaches a 
parabola. Both limiting curves for a — s- and a — > oo are symmet- 
ric, while for generic values of a they are not. The reason why for 
a — ^ the curve becomes symmetric is due to a diverging symmet- 
ric contribution to <^J^ e"' " Uq_tjU_q,t^^, not due to a vanishing 
of the asymmetric part. 



to 



24 



2(e*i-*3 - l)2(e*2~*3 - i)2g(9^+i)(to-ti-t2-i-t3) 



= Sw 



e9^(to-t3)(e*o-*3 _ 1)2 
g(t„™t3)(2,^+i) 2; 2 - 2g2; j-^) 

(e*o-*3 - 1)2(1 -2q2) 
g(t,^t3)(2,^+i) (^1^ 2; 2 - 2q^; 

(e*i-*3 - 1)2(1 _ 2g2) 



(293) 



This is a joint expectation value conditioned to the event that 
the avalanche ends before ^3. 



As in Section IIII CI we obtain the conditional average s.t. 
the avalanche ends exactly at by taking a derivative dt^ of 
the above average ( |292t . and dividing by the total probability 
-Pduration(t3 — tQ)6f foi the avalanchc starting at to and to end 



Here we use the Kronecker .symbol = 1 or according to whether 
« = or not, i.e. the characteristic function for the event ii = 0, which is 
dimensionless. 



at ts, leading to 

2^U;(e*i-*3 _ l)2(e*^-*3 _ l)2g(g^-(-l)(2t3-tl-t2) 

p,(to-t3)(2,^+i) (1, 2; 2 - 2<z^ j-L^) 
(2^2 - l)(e*o-*3 - 1)2 
g(t,-t3)(29=+i) ^F, (1, 2; 2 - 2q2; 



(2^2 - l)(e*i-*3 _ 1)2 
2 2 



1 



0^2—^3 



coth 



a2q^{to-t3)+to+t3 
(gto _ 2*3) 2 



gio — ^3 — 1 
g2g2(ti-t3)+ti+t3 

(e*i — e*3) 2 



(294) 



The resulting function, for t2 ~ ti and T = 1 is plotted on fig- 
ure[8] One sees again that higher wave-vectors q are (slightly) 
skewed towards earlier times. 

It is interesting to perform the same calculation in real 
space. One can either Fourier transform the above result 
(which is not easy) or go back to Eq. ( 1292b and directly work 
in real space. Because of divergences indicated below, we 
need to compute the more general observable, smoothed on a 
small region of space (i.e. for close-by points xi, 0:2): 



2a\/TT 



(aq) 



Integrating over momentum directly in c? = 4 we obtain 



(295) 



{Uq,tiU-q^t^) e 



(aqf 



9 



87r2 



■ sinh' 



V 2 



sinh2(tL^)sinh2(^2^) 



sinh^( 
dt' 



2/t3-to 



sinh2(*a^)(a2-2<' + ti+t2)2 



(296) 



For ti < t2 we can set a = and obtain a finite result. How- 
ever, for equal times ti = t2, there is an ultraviolet divergence 
and the integral diverges like 1/a as a — ?• 0, hence we must 
keep a > 0. This allows to define a (normaUzed) second shape 
at time ti as the ratio 



S2(ii 



-q,tl/ 



-(aqf 



Iq {<t^U.qjJ e-i-ir 



ta + ts 



(297) 

This is the second shape normaUzed to unity for ti = the 
mid-time of the avalanche. The result is plotted on figure |9] 
where the integral over t' in Eq. ( 1296b was performed numer- 
ically. Note that upon normalization the limit a — > exists 
(even if both numerator and denominator diverge) and is a 
parabola. Another possibility to regularize the function is to 
chose ti <t2\ the role of the parameter is then replaced by 
the difference t2 — ti. 
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For the Brownian force model, the tree theory remains exact 
below due = 4, hence we can use the formula in any d. Upon 
integration over momentum in d < 4, the factor (a-^— 2t'+ii + 
t2)^^ is replaced by (a^ — 2t' + ti + 1-2)^'^^'^ ■ In dimensions 
(i < 2, the limit a — > can be taken. In smaller dimensions, 
the asymmetry is less pronounced. This is expected, since 
for d — we must recover the result for the particle, where 
V? = u^- The same holds true for LR elasticity. 



6. Arbitrary elastic kernel and non-local elasticity 

Finally we can now give the result for an arbitrary elas- 
tic kernel, g~^. Since we often use dimensionless units, we 
must first define g'^^ := gq=ogq^ = dq/tn^- Thus one has to 
substitute ^ g^^ — 1 in all above equations containing q 
expHcidy, e.g. Eqs ( |285] ), ( |286] l, ( |293] ). and ( |294l ). 

The equations where q has been integrated over, such as 
( 1296b . have to be recalculated. There the changes to be made 
can be condensed to a change of the integration measure over 
momentum. For the simplest form of a long-range elastic ker- 
nel this is explained in section HVGI 



So + Sdis, So given in Eq. (ISTT ) and 

5dis = -W UxtUxt'dtdt'A{v{t -t') + Ua;tt') ■ (301) 

Jxtt' 

We now rewrite this term with no approximations as 

Sdis = -(J UxtUxtiv + iixt) (302) 

Jxtt' 

+ \ ( UxtUxt'{v+Uxt){v + Uxt')^'Lg{'"{t-t')+Uxtt') ■ 
^ Jxtt' 

We have defined 

A{u) = ^a\u\ + A,eg(w) , (303) 

such that A"eg(u) is the second derivative of A(m) without the 
(5-function part; hence A"eg(0) = A"(0+), and Ai.cg(M) has 
a regular Taylor expansion in |u| around zero starting at order 
|up . Below we loosely denote A"(0) = A"(0+) since the 
right and left second derivatives coincide. 

B. Simplified model 



rv. LOOP CORRECTIONS 

Until now we found that the mean-field theory involves 
only the cusp parameter a = — A'(0+). As was the case 
for static avalanches ItsIi . the small dimensionless parame- 
ter which controls the importance of the loop corrections (and 
thus the deviations from mean field) is the second derivative 
of the (renormalized dimensionless) disorder correlator, i.e. 
using the same notations as in iFtsIi 



A = 
a :~ 

h ■■= 



m''-'*A"(0+) , 
e/2TO-^A"(0+) = -A"(0+) 

d'^g 1 

(27r)2 (l + g2)2 ■ 



(298) 
(299) 

(300) 



The parameter a is of order 0{e = due ^ d). Below we study 
first the 1-loop corrections using a simplified theory, which 
retains only a and A"(0+). This simplified theory streamlines 
the calculations, and allows to derive, at least heuristically, the 
result, which we then analyze. Finally we present a detailed 
derivation from first principles. Note that the presentation here 
focusses on standard short-range elasticity, i.e. due = 4. The 
generalization to LR elasticity is straightforward, so we only 
detail the main features in section lTV Gl and give more explict 
formulas in appendix U 



A. General framework 



The decomposition ( 1302b is exact. Now we make a sim- 
plification. We neglect the higher derivatives A'^")(0+) with 
n > 3. We will see below that this is sufficient to give the 1- 
loop result for the generating function almost completely, up 
to some subtleties that we discuss below. With this assump- 
tion, we have Sdis = '5^ + with 



osimp 

'-'dis — 



-cr / iLl^{v + iixt) 

xt 

+ l-A"{0) [ UxtUxt'{v + Uxt){v + Uxt') ■ (304) 

^ Jxtt' 

We now work with this "simplified" model, and discuss later 
on the effects of the neglected terms. 

The nice feature of this simplified model is that the new 
term can be written as an average over a fictitious (centered) 
Gaussian disorder 7]x with correlations 



{nxVx')r,=m^-^AS^{x-x') 



(305) 



where A is dimensionless, and we will choose later A 



-.d— 4 A" 



A"(0). With these definitions one can write 

G[A] = {G,[X])r, 



25 



(306) 



with 



G„[A] = J D[u]D[u]er^'^+^-''^^^'^''+"^'^ (307) 
Srf=So-(7 / ultiv + iixt)- / r]xUxt{v + Uxt) ■ (308) 



In order to compute the generating functions (|77l ) and ( |78] l 

beyond mean-field, let us start again with the dynamical action 25 ^ote that the noise r;, is unrelated to the friction r; despite the coincidence 

( ISTT i of the velocity theory, which we recall has the form S = in notations. 
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For each realization of rjx, the theory has the same features 
as the mean-field theory (|87] | of Section IIIIBI In particular, 
the total action (including the sources) is linear in the velocity 
field. Integrating over the latter, as in Section (IIIIBI ) one finds 



(309) 



The quantity u^^ is now solution of the (modified) instanton 
equation 

ivdt + V2 - m^)u':^ + a{u'jr - -A.,, - rj^u^J , (310) 

which has an additional "random-mass" term. Using this 
equation, Eq. ( |309t can be written as 



^ ^vL''Z„l\] 



Jxt 



xt 



-Xri 

HI I 



To lowest order in v we thus find 

i;=0+ 



Z[\] 



m 

17 



(311) 



(312) 



(313) 



As we discuss later, we will need to take ^4 < at the fixed 
point, hence the sign of the random term ( 13051 ) is not consis- 
tent with an additional real disorder Since all we want to do 
here is perturbation theory in A"(0), more precisely in the 
parameter a — 0(e) defined in Eq. ( 1299b . this is immaterial. 
It should be considered as a trick to simplify the perturbative 
calculations. 



C. Perturbative solution 

1. General equations and formal solution for arbitrary \xt 

For simplicity we switch from now on to dimensionless 
units, which amounts to setting 77 = m = cr = 1. We want to 
solve perturbatively in rjx the equation 



(314) 



We expand the solution in powers of rjx, denoting by u"^ the 
term of order 0{rf'), 



Uxi ^Uxt+U^t+Uxt+ - 

One must thus solve a hierarchy of equations 

1] = -A.* - 

'^xt \ ^xt 



(315) 



[dt + Vl^l + 2^0,] ul = -{ul^f - i^^ul, . (318) 



V2 - 1 + 2^0 1 u\ - - ^-'O 



(316) 
(317) 



The first line, for order zero, is the usual (mean-field) instan- 
ton equation ( |9l1 l. This perturbation problem is distinct, but 



similar, to the one studied in Section lTlI HI We introduce again 
the dressed response kernel (1268b . now in dimensionless vari- 
ables, 

{-dt - + 1 - 2ul^ W^.f^xt = S^'ix - x')S{t ~ t') . 

(319) 

It has the usual causal structure of a response function, and 
obeys a backward evolution equation. It allows to rewrite the 
solution of the system of equations ( 1316b to ( 1318b as 



rjx' u^i^iKx't' ,xt , 



(320) 



x' Jt'>t 



xt~ I I [{ult')^ +Vx'ul,t,]^x't'.xf (321) 

x' Jt'>t 



Consider now the average (13131 ) over rjx using ( 1305b . i.e. in our 
(dimensionless) units {rixi]y)n — AS'^{x — y). Since {ul.^)n = 
0, the lowest-order correction is given by the average of u^^, 



Z[X] = Zu-oclM+L-'' f 

Jxt 



^xt)n 



(322) 



Inserting Eq. (13201 ) into Eq. ( 13211) . and performing the average 
over 1], one finds 



xtin 



A 



t<tl<t2,t3 J X\x' 

X ^X>t'2,Xxt^x'ti,X\t'^XA_t\,Xt 

A I I u^, ^f^M.x't2,x'ti^x'ti,xt ■ (323) 

It<ti<t2 Jx' 



It admits the following graphical representation 





xt/V 



(324) 

The symbols are as follows: (;) a wiggly line represents u'^^, 
the mean field-solution; (//) a double solid line is a dressed 
response function K, advancing in time following the arrow 
(upwards), thus times are ordered from bottom to top. Note 



that for the choice Xt = XSIf ), one has {t° 
t > 0, hence the integrals only involve negative times. 
We now define the combination 



Ux't'^x't' ,xt , 



for 



(325) 



t'>t 



in terms of which one can rewrite 

^{y,x',t'f + '!?ix\x\t') 



{^xt)n 



^x't'.xt 



(326) 

In section lTV F 2l we shall show that there is an additional term. 
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2. Space-independent source, X^t = At 



D. 1-point velocity distribution 



We now pursue the calculation in the case of a spatially 
uniform source X^t = Xt, i.e. we study the center-of-mass 
velocity. Since then u'^f ~ is also uniform, we can express 
the solution of Eq. ( |319t - as in Eq. (|269b - in momentum 
space 

^UMM = e-^'^'+i^^^-^^^+^^n ^^<e{h - h) . (327) 



The same is possible for Eqs. (13251 ) and ( |326t . by defining 

<P{x',x,t) = /^$(fc,i)e^'=(^'-^) and 



'tl<t2 



Jt{k) = f [^k,h)^ + ^k,h] 
Jti>t 



(328) 
(329) 

lk=OM,t ■ (330) 



From Eq. (|313t we find that Z\X\ is then given by 



Z[\] = Z'^^^[\]+Aj^J{k) (331) 



J{k) = jJt[k) 



(332) 



As discussed below, some counter-terms are missing, and the 
correct formula is obtained by J{k) J{k) + Jct{k)- 

We now consider the space dimension to be d w 4, since 
we want to perform an e = 4 — c? expansion. Since A ^ e, 
it is sufficient to calculate J{k) in d = 4. In that case, 
we note that for any isotropic integral one can write (recalling 
A = -m'^-^i:^"{Q), and Eq. ( |299l )) 



d'^k _ a 



We used that 

A = i 
5*^ 2 



d(fc^) 



2, {k^f-^'^ 1 



(fc2 + 1)2 



(333) 



(334) 



where . . . denotes higher-order terms in e and Sd the unit- 
sphere area divided by (27r)'*. 



1. Generating function Z{\) and moments 

We now specify to \{x,t) = \6{t) to obtain the 1-point 
velocity distribution. 

Let us recall the solution of the instanton equation 



We found useful to define 



ei-t) 



1 



1 - K^ie 



lo-t 



e{-t) . (335) 



1-A ' 



(1-«)(1-A) = 1, 



(336) 



which we often use below as it simplifies the calculations. The 
relevant interval A € ]— oo, 1[ maps onto k S ]--oo, 1[ (with 
reversed boundaries). 

From the previous section we have 



Z{X) = Zo{X) + -5Z{\) 



(337) 



6Z{\):= ed{k^)[J{k,K)+J,,t{k,K)], (338) 
Jo 



where we denote Z'"'^ Zq and J{k) in (l332l ). ( l330l l by 
J{k^ k) to make the k dependence explicit. The calculation 
of SZ{X) then proceeds as follows. We need the dressed re- 
sponse only for ^2 < (since vP vanishes at positive times). 
It reads 



g-(fc^+l)(t2-tl)('gi2^ _ 1)2 

t,t2,ti = n -TT? -^(t^ - ^i) ■ (339) 



(e*iK-l)2 



This yields for ti < 



fc2^e*i + e*^'*! (fc2(l - k) - 1) + 1 - /c2 



A:2(fc2 „ l)(e*i«;- 1)2 
with$(fc,ti) =Oforti > 0. 



(340) 



t<ti<0 



^k,tifRk=o,tut 



t<ti<0 



[fc2(K - 1) + l] K-^^B^{k'^ + 1,0)+ fc2K + hl(l - k) 
fc2(fc2 _ 1) 

2fc2(fc2_l)2 {2 (2fc' + 1) - 1) + 1] ' n--'''B, (2fc2 + 1,0) 

-6 (fc2 + 1) [fc2(K - 1) + 1] n^^^'B^ (fc2 + 1, 0) 

+ ^2/4 [2fc2(K - 1) + K - 4] - 2(fc2 + 2) ln(l - ti)} . 



(341) 



(342) 
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We have introduced the incomplete beta function B^{a,h), 
defined as 



sion of the counter-terms: 



B^{a,b) := r<«-i(l-t)^-i 
Jo 



(343) 



and related to the hypergeometric function 2F1 via 



j3,K0)= ""^^^(^'°'«+^'") , (344) 



which can equivalently be used. Note that while Bi^{a, 0) has 
a branch cut for negative k, it is a spurious one since in our 
results only the combination K^°_BK(a, 0) appears, which is 
perfectly regular on the negative k axis. 
The final result for J{k^ k) is 



2fc2(fc2 _ 1)2 

X I - 4 [k^{K - 1) + 1] K-''^B^{k^ + 1, 0) 
+ [fc2(K-2)K-21n(l-K)] 
+ 2 [k^{K - 1) + 1]^ K^2fc"s«(2A:2 + 1,0) 

The special cases we need are of the form 

't\'' 1 



(345) 



K~''B^{l + x,0) ^ [ dt 
Jq 



kJ 1 - 1 



(346) 



Taylor expanding the denominator 1/(1 — t), and then inte- 
grating leads to a very useful series representation 

K-^S„(1 + X, 0) = V — — = K^K, 1, X + 1) . 

^0^ + ^ + 1 

(347) 

$ is known by Mathematica as the HurwitzLerchPhi function. 
Using the above series expansion, one can easily obtain the 
small- and large-fc behaviour: 



J{k,K) ^ -K+^K^ + 0{k'^) (348) 

K K + 21n(l-K) ^ / 1 \ 
J{k,n) = -- ■ ^^^"^^ 



Hence J'{k,K) is integrable with the measure k'^d{k'^) at 
small fc, but contains a quadratic and a logarithmic divergence 
at large fc. These will have to be cancelled by the counter- 
terms, leading to a finite result. We will show in section HV Fl 
that the exact expression for the counter-term is 

^ ^ (3 + fc^)^ + 21n(l-K) 



Using the series expansion ( 1347b . the integration over fc can 
be performed, keeping a large-fc cutoff in the intermediate ex- 
pressions. The final result is after simplifications, and inclu- 



SZ{X) 



kM{k'')[J{k,K)+J,tik,K)] 



= K^(l - hi 4) 



00 



n— 3 

(n-3)(n-2)2ln(n-2) 

61n(2) - 2ri(n + l)(ln(2) - 1) 
n-^(n + 1) 

{n - l)(n((n - 6)n + 2) + 6) ln(?i - 1) 



(351) 



(n2 



n?{n + 1) 
8n + 3) ln(n) 



2(?i + l) 



(352) 



Note that lim„^.2 a„ — 1 — 2 In 2, i.e. the first term 02 follows 
the same relation, if the coefficients are properly interpreted. 
For later convenience we set oi = 0. 

It is also possible to calculate the cumulants of the velocity 
directly in a perturbative expansion in the full disorder to 1- 
loop accuracy. This is performed in Appendix IM 1 1 up to the 
third cumulant. We have checked that this indeed agrees with 
our explicit series expansion up to order A"^. As the reader 
will see, the calculation of the appendix increases formidably 
in difficulty with each new order, while the present method 
allows to sum these diagrams much more efficiently. 

It is interesting to give the lowest moments. In dimension- 
less units they read, expanding ( 1351b in powers of A using 
(13521) and (l333l l. 



(353) 
(354) 
(355) 





~ V 














14 


-a{l 


-ln4)] + 0(1-2 ) ^ 






= V 


2 - 




+ 9 In 3- 131n4) 




i4 


= V 


6- 


a9 
^25 


(20- 132 In 2 + 69 In 3) + 0{v'^ 



(356) 



We recall the mean-field result = (p — 1)! which follows 
from Z{X) = — ln(l — A). The general formula for the mo- 
ments p > 2 is easily obtained as 



p 



i-iy 



(357) 



Let us recall that the small parameter a is related to the second 
derivative at the fixed point and equals (see (B12) of ItsII ): 

a = -A"(0+) = -i-^ + O(e2) 



1-C: 



i. + o(a 



with Ci = 
class, i.e. 



1/3 for the RF class, and (1 
2 



(358) 
(359) 

for the periodic 



(RF = nonpcriodic disorder) (360) 
(periodic disorder) . (361) 
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FIG. 10: Plot of 5Z{X) defined in equation dSSB . For real A the 
function (5Z(A) is defined for A £] — oo, 1 [, with a singularity at A = 
1. The Taylor series in k has convergence radius 1, which translates 
into convergence for A from — oo to 1/2. This is plotted with an 
upper bound on the series rimax ~ 2\ for i = 1, 2, 10, starting 
with the dotted curve for i = 1, dashed for i = 2, the remaining 
ones solid, and finally i = 10 (fat). This establishes convergence 
in that interval. The last 6 curves are indistinguishable, except at 
A 1/2. See figure[T3lfor the resulting function. Note that 5Z{X) 
can be obtained numerically from the alternative formula ( |H6t with 
excellent agreement. 



Hence for non-periodic depinning 

^ = w(l + 0.0858432e) 
^ = w(2 + 0.014924e) 
i4 = v{6 - OMUGAe) 



(362) 
(363) 
(364) 



More ambitiously, we will now determine the correction to the 
velocity distribution in an avalanche. 



2. From Z{X) to P[u): Distribution of velocities in an avalanche 



The series for SZ{X), defined in Eq. ( 1351b . as a series in 
K, has convergence radius 1, since a„/a„+i — > 1 at large n, 
equivalent to 3?(A) < i. This is demonstrated on Fig. [lOl 
The physical singularity however is outside of this interval, at 
A = 1. 

We now obtain the avalanche-size distribution. As ex- 
plained in section nil C II we have 

P{u) = (1 - p'v)6iu) + vp'V{u) + 0{v^) (365) 



with 



Z{X) = p' / du(e-^" - l)V{u) 



(366) 



We have obtained the expansion of Z{X) to order 0{a) in the 
form (13371 ). hence 

1 • a 

p'Viu) = -e-" + -SV{u) , (367) 
u 2 

/>oo 

6Z{X) = / du(e^" - 1) dViii) . (368) 
Jo 



In the case where 5Z{X) admits an inverse Laplace transform 
we can also write 



V{u) = Vmf{u) + -^SViu) . 



For ii > the inversion reads 

"-'^ dX 



5V{u) = 



27ri 



6Z{X)e 



-Xii 



(369) 



(370) 



the contour being closed to the right. Note that 

5Z{\ = 0) = 
5Z'{X = 0) = 



d?i(57'(u) = , (371) 
Auu5V{u) = Q. (372) 



To construct the probabiUty distribution we first note the in- 
verse Laplace transform 

LTI^^.k" = 5{u) - n iFi{l + n,2,~u) 

= (5(u)+c-"9,iL„(7i) (373) 

in terms of the hypergeometric function iFi, or equivalently 
the Laguerre-polynomial L„. For u > it can be found by 
rewriting the contour integral (which with our conventions 
must be closed to the right): 



?.oo 

d_ 

dv 
d_ 
dii 



dA 
27ri 



-A 



n 



-Xu 



1 - A^ 



da 



-\u-a{l-\) 



„n-l 



2TTi 



T{a) 



r(n) 



leading to ( I373I ). Thus we can now write the formal series 

oo 

6V{u) = a«e-"9,iL„(u) . (374) 

n=2 



Unfortunately, this series is divergent. 

This problem can be cured as follows: We will subtract 
from the series (13511 ) terms which can be summed analyti- 
cally, resulting in polylogarithmic functions, and their deriva- 
tives, and inverting the latter via a cut-integral. These terms 
are chosen to render the remaining sum (quasi-)convergent. 
To this aim, we note 

SZiX) = SZs,,iX) + <5Zcut(A) . (375) 

We start by Taylor-expanding a„ around n = oo, 

_l-41n(2n) 2 61n(2)-g 5-61ii(2) 



2n 
61n(2) 



101 
30 



O 



2 Inn 

n 



1 



OO , 



(376) 
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1 
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The coefficients a„ are what remains of a„ after subtracting 
their asymptotic behavior, 



2 ) — 



(382) 



Especially note that a\ becomes non-zero, even though a\ 
0; in fact, this coefficient grows rather quickly with j.^. 
while the other coefficients decay. 



(5Pser(u) 



oo 

E 

r(=l 



a„e ■"9i,L„(u) . 



(383) 



FIG. 1 1 : uSP{u) as given by Eq. l |380t . or, equivalently by Eq. ( |H7 



with bi = i — 2 In 2, 62 = 2, and an explicit formula for bj for 
j > 3 is given in Appendix iNl We recall that 02 = 1 — 2 In 2 
and that we set ai = Performing the summation, we 
obtain (if the series converges) the alternative representation 



(377) 



Lij (k) is the polylogarithm function, which is analytic on the 
complex plane with a cut on the real axis for k E [1,oo[, 
which maps on the same interval for A e [1, oo[ with reversed 
boundaries. It is along this cut that we have to integrate. The 
discontinuity there is given by 



lim hij{K + ie) — ljij{K — ie) 



0, 



K> 1 , 
K <1 . 

(378) 

Note that this also holds true for the derivative w.rt. j, and for 
j = 1, i.e. Lii(K) = — ln(l — k). 

Thus, the inverse Laplace transform (13701) becomes a com- 
pact and simple cut-integral 



27E + 21n(lnK) + ^ 



(379) 



However, this series also diverges. Therefore we choose j„ 
as a cutoff, by defining 



SVcutiu) = 



dA 



27E + 2 ln(ln k) + ^ ^ 



(380) 

bjilnny-^ 



(381) 



Although this may appear to impose an artificial constraint bj = 

it will be immaterial in what follows since we will use only a finite sum 
and add and subtract the same terms. 



Both expressions, SVcut (") and 6Vsci {u) can be obtained nu- 
merically with good precision, and seem to decay rapidly at 
large ii. One then checks that the sum of the two, for any ii in 
the bulk of the distribution, converges extremely well versus 
the result at j = jmax, e.g. for u = 1 excellent precision is 
already obtained for jmax=3- Of course, for a fixed jmax the 
sum over n in ( 1383b should be stopped at n not too large since 
it is an asymptotic series, which is ultimately divergent, but 
in practice the range of convergence (with respect to n,nax) is 
rather broad. 

Practical values are jmax = 15, and ( |383t can also be 
stopped at n = 15. With this choice, we find that the pre- 
cision is excellent and that all moments uPdP{u)du be- 
tween the fourth and 36th are at least given with a relative 
precision of 10^^, most even of 10^^°. jmax should not be 
taken too large, since otherwise this shifts too much weight 
into the moment di, leading to numerical problems (cance- 
lation of large terms.) As an example, for jmax = 15, one 
has 51 = -51.97, ^2 = 0.002976, £3 = 1.359 x 10"^ . . . , 
020 = 2.373 X 10~^^. There are no convergence problems at 
small or large ii. 

The final result for is 



'a5V{u) 
2P^ 



(384) 
(385) 



where we remind the value of the small parameter a from 
(13581 ). Note that the second formula (13851 ). while being equiv- 
alent to order e, has the property to resum the logarithmic be- 
havior at small v into the correct power-law behavior This is 
why we have chosen it in Fig. [T2l 



3. Small-velocity behaviour, the critical exponent a 

Let us now obtain the small-velocity asymptotics of P{u) 
and extract the a-priori new critical exponent a. It is controlled 
by the asymptotics of 5Z{X) at large negative A, i.e. A — > 
—00. This corresponds to the behaviour for n ^ 1^ of the 
series ( I351t . It is determined by the leading behaviour of a„ 
at large n, i.e. from the leading term a„ ~ —2\n{n)/n of 
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FIG. 12: Left: V{u)u'^ as a function of ii^ for RF-disorder; e = i.e. MF (blaclc dasfied), e = 1 (blue, dotted), e = 2 (green, dasii-dotted) 
and e = 3 (red, solid). Right: log-log plot of T'{u) as a function of ii, for the same values of e. For both plots resummation formula ( I385t was 
used. We note that since only a appears as a parameter, at this order e = 3 RF and e = 2 RP are indistinguishable, see Eqs. l l360t . i36U . 



( 1376b . Resumming with this term alone, we obtain 

5Z{\) = ^a„K" « -2^ ^^k" = 2daUa{K)\^^^ 
n=2 n=2 ^ 

,2/ 



(386) 



= -ln'(l-A) + ... . 
This yields for A ^ — oo 

Z{X) = ZoiX) + ^SZiX) 



In(l-A) [l + -ln(l-A) + 



(387) 



It is easy to see that this is consistent with a modified critical 
behaviour at small velocities 



SZ(X) 
2.0 r 



(388) 
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FIG. 13: (Color online) SZ{\), as obtained by ( [3511 (thick red line), 
a reexpansion in A (dashed thick green line) and numerical integra- 
tion of Eq. l l368t . using Eq. J380t (thin black line). All functions 
agree in their respective area of convergence. 



To show this, we start from the trial probability 

1 



triaUU) 



■,1+x 



(389) 



for which the associated Z{X) can be computed exactly via a 
Laplace transform, using ( 13661 ). Expanding the result in small 
X yields 



Atrial (A) — 



dii- 



1 



,1+3: 



"(e^" - 1) 



= -ln(l - A) 



-7E + ln(l - A) 



A) 



(390) 



the first term is Z*'°°(A) and the second one the correction. 
Comparing the behavior at large negative A of Eqs. ( 1387b and 
( 13901 ), we can thus identify x ^ a, consistent with Eq. ( 13881 ). 
Note that multiplying (13891 ) by a prefactor ^ I + 0{x) 
or changing the exponential to e~"[i+o(2:)] produces only ^ 
xln(A) terms, subdominant w.rt. the ln^(l — A) at A — — oo. 

Let us now discuss our results for the small-velocity expo- 
nent. Using ( 1388b . together with ( 1360b and ( 13611 ). we find 



1 



a = 1 - 



nonpcriodic 
periodic . 



(391) 
(392) 



Our predictions for the change of a are thus quite large, and 
tend to reduce the exponent. A naive extrapolation to d = 1, 
e = 3 (depinning of a line) would suggest a sa 1/3 signifi- 
cantly reduced from the mean-field value aMF = 1- Prelimi- 
nary numerical results indi cate that the exponent a may even 
be negative in d = 1 flOOll . A 2-loop calculation (or higher) 
would settle the question from an analytical point of view. 

We can compare the above formula to the one for the dy- 
namical exponent to one loop 



(393) 
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Hence we could also write 



Indeed that would imply 



z-l + 0(e2) , 



(394) 



which holds for both periodic and non-periodic systems. 
Again it would be interesting to obtain the higher-loop cor- 
rections, since we did not find any general argument why they 
would be absent. 

Finally the small-u behaviour can be studied more system- 
atically. This is done in Appendix |0] where we obtain the 
amplitude at small u > 0, 



P{u) « 



C = 1 



2(47E + 6i) 



(395) 



where bi = i — 21n2 as defined above. This yields C ~ 
1—0.71 1 284a in good agreement with our numerical Laplace- 
inversion. In principle this amplitude is universal and can be 
measured. 



4. The behavior of SZ{\) for A — !> 1, and tail ofV{u) at large u 

The behavior of Z{\) in the limit of A — > 1, which controls 
the tail of V{u) for ii oo, is obtained in Appendix IgI The 
final result is 



Z(A)=:-ln(l-A) + -(5Z(A) 
5Z{X) 



8(1-A)2[ln(l-A)]2 



(396) 
(397) 



To obtain the tail of V{u), one needs to inverse Laplace tran- 
form Z{\). Before doing so, let us point out that this form is 
incompatible with the naive expectation of a stretched expo- 
nential at large velocity. 



C" 



-Bii" 



(398) 



with C" = B = a' = (5 = 1 in mean field (e = 0). While it 
would be hard to extract £?, C and a', we could extract 5 as 
follows. Expanding near (5 = 1, we find 

'^5V{u) = -((5- l)e-"ln?i + 0(((5- 1)2) . (399) 

This is equivalent to 

f ^^(A) ^ {5 - i)Mi_^)±I^ + oiiS - If) . (400) 



Clearly, this is not of the form (I397l l. Noting s := 1 — A, we 
claim that Eq. (13971 ) is equivalent to 



5V{u) ~ -e-"uV(w) 



(401) 



at large ii, where f{u) has a Laplace transform /(s) := 
d?i/(u)e~''" which behaves at small s as 



/(s) = /(O) 



In s 



(402) 



8s2(lns)2 



(403) 



for small s, which is exactly the result ( 1397b . It is then easy to 
guess that 



f{u) 



d__l_ 
dii In ii 



TTTr^ (404) 



at large u, for ii > uq. Indeed, the contribution for ii > uq 
reads 

1 

Inu 
1 



, s e 



Uo 



dw e 



UqS 
1 

Ins 



In K7 — In ^ 



(405) 



In the partial integration from the first to the second line we 
have dropped a term (1 — e^^"°)/ Iniio, which is of order s. 
In the last step, we have used that for s 0, first In w — In s w 
— In s, and second uqs — > 0. 

For the velocity distribution at large ii, we thus finally ob- 
tain 



-Tr)l — loop 

it— foo 



U 



2 



+ 0(a2) 

+ 0{a^) (406) 



We remind that a < 0, which has motivated us to write the re- 
sult in an exponentiated form. Other forms are however possi- 
ble, such as corrections to the pre-exponential only. The form 
( 14061 ) renders the tail stronger decaying; it is plotted on fig- 
ure[T4] In all cases, given the smallness of the correction, this 
tail will be hard to see in numerical simulations. With the 
help of Eq. (13801 ), we have been able to evaluate 5V{u) up to 
ii w 100, while the alternative representation (IH7b works up 
to M ~ 10. For these values of ii, the tail-behavior (14061 ) is not 
yet reached. 



J. Alternative approach: Integrating over momentum first 



Our result for 5Z{\) given in Eq. ( 1351b . was a compact se- 
ries expansion from which one first had to extract the asymp- 
totic behavior at large A, before being able to perform the in- 
verse Laplace transform. A complementary approach, per- 
formed in detail in appendix|Hl is to start from Eq. (1330b . cal- 
culate $(fc,i) as given in Eq. ( 1340b . and first integrate over 
t and fc, the final result for 6Z{\), given in (|H6b . is now an 
integral over ti. (Recall that above in Eqs. ( 1341b and ( 13421 ) 
we integrated first over t, and ti leaving the k integral for the 
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\n(uP(u)) 




FIG. 14: The exponentiated version ( I406t of the function 
^ii^-'oD ("") ~ ~i (^-S- disorder in d = 1) (red dashed 

line), compared to the mean-field result (i.e. e = 0, blue, solid line). 



end). We did not succeed in performing the final integral over 
ti analytically, although it is easy to compute numerically. It 
confirms the above results for SZ{X). The advantage of this 
method is that the inverse-Laplace transform can be done ex- 
plicitly, yielding a (relatively complicated) integral represen- 
tation (as integral over ti) of 5V{u) given in ( |H7l i. It confirms 
all statements made above, including the asymptotic behavior 
for small and large ii. 

Note that in Eqs. (13411 ) and (13421 ) one can interpret t as the 
time of a kick (infinitesimal step in the force), or starting time 
of the avalanche, while time zero is the measurement time. 
The time ti < is an intermediate time, which must be inte- 
grated over the duration of the avalanche. Hence, if we instead 
integrate over k and then ti G [t, 0] at fixed t we obtain the 
joint probability that u(0) = ii and the avalanche started at 
t. Although it is a straightforward generalization we will not 
give this result here. 



E. Recovering the avalanche-size distribution to one-loop 

As discussed in section UlIFl to recover the avalanche-size 
distribution, one can use a source constant in time Xxt = A 
during a large time window T. The avalanche-size generat- 
ing function, noted here ^^(A), is obtained from the dynamic 
generating function studied here via Z[ A] = TZs{\). In prac- 
tice it amounts to suppressing the final time integral in the ex- 
pression for Z[}]. 

For a source constant in time, the solution of the unper- 
turbed instanton equation (77^ = 0) is 



iP = Zl = Zl{\)^-{l-y/T^). (407) 

The dressed response kernel then becomes 

^k.t^M - e-(^''+i-2^°)(*^-*i)0(t2 - h) , (408) 
which is simply the bare response up to the replacement 



^ ^2 _ 2Zg(A). The formula (13281 ) then gives 
ci>(fc, t,) = Zl I ^kMM = ,2^f^ ^^0 ■ (409) 

Following the steps in Section irVC2l this leads to 

Zs = Zl + Zl + ... (410) 

4 = (^i*^, = 4 / ^'^''^ (411) 
ei2 Jk 

Jtik) = j^^^Jtik). (412) 
The coefficient A ~ ~, and we have defined 

eh 

which is time independent. Graphically Eq. (14121 ) can be writ- 
ten as in (1324b . 




(414) 



replacing the external wiggly lines of Eq. (1324b by the fac- 
tors Zg. Note that we have recovered Eq. (152) of ll74ll for 
the statics, up to the two counter-terms discussed below. For 
pedagogical purposes, we want to make further contact with 
the self-consistent equation obtained in I73il . To this aim we 
rewrite Eq. (14121 ) as 

4 / J,'^(fc) = 4(l-2Z°) 

= (Z5-Z0 + ...)(1-Zg-Zs + ...) 

= Zs- {Zsf {Z% - {Z%f\ 

= Zs- {Zsf - A . (415) 

Note that by going from the first to the second line, we have 
added in each parenthesis a subdominant term. From the sec- 
ond to the third line, we have regrouped the terms, and finally 
from the third to the fourth line we used the exact relation 
- [Z^sT = ^- Eq. ( 14151 ) can thus be written as 

Zs = \+{Zsf + ^ I Jt^ik). (416) 




where we recall the graphical interpretation of each term. 
(The amputated lower response gave the factor of 1/(1 — 
2Zg).) Comparison with formula (151) in [7^] shows that 
one recovers the result of the static calculation, provided (/) 
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one replaces in J^{k) the tree generating function by Zs, 
which does not make a difference at this order; (//) one adds 
to Eq. ( |416t two counter-terms, discussed below. 



Zs = \+ [Zsf + 



a 

7h 



(fc2 



1) 



/j2 + 1 



(417) 
(418) 



In the statics these counter-terms appeared naturally by using 
everywhere the improved action. The first one comes from the 
renormalization of A(u), thus all parameters which appear are 
renormalized ones. The second also appeared naturally in the 
statics from the definitions used there, while here it comes 
as a correction from using the (over)simpUfied model, as is 
explained below. 



F. Counter-terms and corrections to the simplified theory 

1. Counter-terms from renormalization 



In II73I1 the static avalanche- size distribution was computed 
using the improved action, i.e. in terms of the renormal- 
ized disorder A(u), which automatically includes the counter- 
terms for the renormalization of the disorder In the dynamics, 
there is an additional operator which is marginal at d = due 
and corresponds to the friction term in the dynamical action. 
Computing from the start in terms of the renormalized fric- 
tion T] is possible, but less convenient, hence here we perform 
the calculation first in terms of the bare disorder Ao(w) and 
the bare friction rjo, and then reexpress at the end the result 
in terms of the renormalized disorder and friction. This yields 
an expUcit derivation of the counter- terms. 

We start from the bare action given in Eqs. (ISTI ) ff. 



S = 



+Aq(0+) / UxtUxt{v + iixt) 



+ iAo(0+) / UxtUxt'iv + u^,t){v + iixt') 

^ J xtt' 

Here the subscript zero denotes bare quantities. 

The effective action to one loop, F = 5 + 8S, reads 

55 = Ao(0+)/ u^t{v ^ u^t')^^^t'y-xt) 

J xtt' 

-2A'o(0+)A^'(0+)x 



XitiXtt' 



UxtUxt'Ux^tA^x^t^Ux^t^Uxt'Uxt) + 



xitixtx't' 



(419) 



where here averages (...) are w.rt. 5o. Corrections to A"(0+) 
were omitted, since they do not matter to this order From F 



we can now identify the renormalized parameters. The second 
term leads to A'(0+) = A^(0+) + 5A'(0+), with 

1 



(5A'(0+) = -3A^,(0+)A[,'(0+) 



(420) 



This is the correct FRG equation for A'(0+) Ii73i1 . The first 
term gives 



+AS(0+) 



1^ k'^ + 



Uxt 



UxtUxt 



(421) 



-A'o'(0+),7o 



1 



Ik (fc^ -r I Jxt 

The last term gives, in agreement with iItsIi . the renormalized 

rj = riQ + Sr], 

1 



UxtdtUxt + 



<5,, = -a;,'(o+),,o 



(422) 



Reexpressing Z{\) instead in bare parameters as a function of 
renormalized ones, defines the counterterms as 

Z(A; r;o, Aq) = Z(A; 77, A) + Z=*(A; 77, A) . (423) 

Using that 

Z{\- ?7, A) = Zt,ee(A; 77, A) + |JZ(A) , (424) 

where a ^ A"(0), and given by Eq. (12991 ), we only need 
to expand Ztroc to first order in the differences 5 A and Srj. 
Eq. ( 11051 ) allows to restore units, 

^.cc(A; 77, A) = ^^Z,., ( Zh^] . (425) 



-A'(0^ 



77777' 



HereZtrco(A) = — ln(l — A) and we remember that A' (0+) < 
0. To compute the rh.s of Eq. (14231 ) we substitute 77 — > 770 = 
77 — 5rj, A Aq = A — (5A, expand to linear order in the 
differences, and in the final result we replace, to this order, 
bare parameters by renormaUzed ones. This gives 

^=.(A;„,A)=r^^-^2)^ (426, 



A'(0+) 77 J -A'(0+) 



~AA'(0+) 

— —ri — 



We now switch back to dimensionless units, setting 77 — > 1, 
777 1 and -A'(0+) 1. Using (1 - A)(l - k) = 1, we 
then find 

Z^'(A; 7;, A) ^ -2A"(0+) ^^j^ [ln(l - n) + n)] . 

(427) 

Comparing Eq. ( |423] l with Eqs. ( |337] |. ( |299] | and ( |333] |. we 
finally obtain 

Note that to derive this counter-term we have used that m = 
nif) i.e. that the mass is not corrected, a property that we now 
discuss in detail. 



J^^ik,K) = 2 K + ln(l-K) 
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2. Corrections to the simplified theory 

Let us examine more closely the effective action derived 
from the simplified theory, i.e. the first two terms in Eq. (14211 ). 
We see that there appears an apparent correction to m^, ob- 
tained from the second line of Eq. (14211 ). 



(S^in^pTO^ = A^'(0+) 



1 



fc2 + 



(429) 



However we know from the STS symmetry that the mass can- 
not be corrected. The reason for this artifact is subtle. Let 
us go back to the exact theory (1302b . When computing the 
effective action, there is an additional term 



6S ^ Uxt / AtlRtix,t'x{Uxt+v){Uxt' +v) 

Jt'<t Jf 

X A;;'g(z;(t - t') + u,t - u,,t,) , (430) 

not present in the approximation AJ4g(ti) = A"(0). Although 
it contains a third derivative (which to this order is not sup- 
posed to matter), it gives a correction. To see this, we recog- 
nize that 

{uxt + v){uxt' + w)Aj."g(u(t - t') + U^t - U^t') 
= -dtdt'K,g{v{t - <') + u,t - u,t') . (431) 
Inserting this relation into Eq. (1430b . we obtain 

SS ^ - [ Uxt [ dTR{T, a; = 0) 
Jt'<t Jo 

xdtdt'A'^^g{v{t - t') + u^t - u^t') . (432) 

We now integrate by part w.rt. t': there is no boundary term 
at t' = t (since the r-integral then is zero); and there is no 
boundary term at t' ~ — oo, since then 9tAJ.j,g(v(t — t') + 
Uxt — Uxt' ) = 0. Thus only the upper bound of the r-integral 
contributes, and gives 

<55 = - y UxtR{t-t',x = Q)dtKogHt-'t')+Uxt-Uxt') ■ 
t'<t 

(433) 



Thus we arrive at 

5S ^ ~ I Uxt{v + iixt) 



XI I dTe-^^''' + ^'>A"{vT + Ua:t-Ux^t-r) (434) 
k 

In the limit of small v, the term A[!^g{v{t — t') + Uxt — Uxt') 
can be approximated by A"(0+), thus 



6S = -A"(0+) / uxtiv + iixt) 1^ . (435) 



Thus there is an additional term 



<5addm2 = -A^'(0+) 



1 



(436) 



This term cancels the spurious mass correction. Note that in 
another derivation, given in appendixlVl both terms ( 1429b and 
( 14361 ) appear. 

Two observations are in order: First of all, one can rewrite 
the two terms graphically as 





(437) 



While the first one naturally arose in the velocity theory, it is 
the second one which we derived above. Their crucial differ- 
ence is where the field ii is sitting in time, as iit at the same 
time t as the response field ut, or as Uf at the earlier time t' . 
Thus there is no correction to m? due to this cancellation, also 
kno wn a s the mounting property (and frequently used, see e.g. 

iaiiol). 

Second, we have used that A" (u) decays to for u — > oo, 
i.e. short range disorder 

Finally, the additional loop correction (1436b must be added 
to our calculation based until now only on the simplified the- 
ory. It can be interpreted as an additional "counter-term" to 
subtract ( |429b . To calculate it let us consider how this addi- 
tional term (14361 ) contributes to Z{\). Indeed, it changes Eq. 
(12681) to 



''Ult] ul 



= -(u-xt) - VxUxt + (JaddTO Wxt ■ 

This is equivalent to an addition to Z{X), equal to 



Su'^xt = -<5addTO^ / U^^IKk=o,tut 



(438) 



t<ti<0 



In terms of J^{k, k) it reads 



1 



(439) 



(440) 



Both counter-terms together give, as already used in Eq. (13501 ). 

^ (3 + fc^)^ + 21n(l-^) 

(1 + F)2 ■ ^ ' 

3. First-principle calculation in u-theory 

We note that the two terms in Eq. (14371 ) naturally appear 
together in calculations based on the position field u{x, t), in- 
stead of the velocity u{x, t), see e.g. Eq. (3.22) and Fig. 9 on 
page 13 of ||48ll . The question thus arises whether one could 
construct the field theory directly for the position field instead 
of the velocity field, and whether this would give directly the 
combination ( 1437b . As we show in appendix|V] both answers 
are "yes". 
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G. Distribution of velocities for long-ranged elasticity 



Then the resuh for 6Z{X) is 



Ahhough some systems with long-range elasticity are stud- 
ied at their upper critical dimension (usually interfaces with 
due = 2), some require an e expansion around d^c- This is the 
case for instance for the contact-line or fracture fronts {d ~ 1, 
fi ~ 1, due ~ 2), i.e. e = 1. We now indicate how the one- 
loop calculations of the previous sections can be extended to 
these cases. 

It turns out that the details of the velocity distribution de- 
pend on the precise form of the elasticity kernel at large scales. 
This was already the case for the statics, and in iFtsIi we estab- 
lished a general formula for the avalanche size to one loop as 
a function of the elastic kernel. This formula was applied in 
in the case of the contact line. 

Although we sketch below the calculation for an arbitrary 
kernel, for simplicity we will concentrate on a kernel of the 
form 



6(g) =.9-1 =c(g2+^2p/2 ^ rn^^cfi-^. 



(442) 



we set c = 1 by a choice of units. The upper critical dimension 
due = 27 is identified by the large-g divergence of 



l2= 9n = Cd.,^fJ. 



,1 



(443) 



Here e = due — d > 0, and Cd,-y = e/2 where I2 = Jg{q^ + 
2^7/2 jjjg rescaled disorder parameter is defined by 



(444) 



At the fixed point, it reaches, in the limit of small m (small 
II), the same value as before, independent of 7, 

1 



a 



-A*"(0) 



,(e- + 0(6^) 



(445) 



Note that the avalanche size becomes 



m-^A'{0+) = M"27A'(0+) = (e/2)^iA'(0+)/i-''+'^ 

(446) 

and we refer to iFtsIi for more details. We now use dimen- 
sionless units meaning that we express x in units of l/i-i, 
time in units of t„i = rjm/m^, and all velocities in units of 
Vm = rn'^Sm/Tm (or v,rt = L"'^ Sjn / Tm)- In these dimen- 
sionless units the result for the center-of-mass velocity does 
not change at the tree level, i.e. for mean-field. We will write 
the 1-loop result for Z{X), or V{u), in the form 



Z(A) = ZMF(A)+a-(5Z(A), 



V{u) 



a-5V{u) 
d 



(447) 
(448) 



inserting the factor of 2 /d for later convenience. For SR elas- 
ticity d = due — 4, and one recovers the previous definition. 

The calculation of Section|lV]is easily extended to an arbi- 
trary kernel gk- All we have to do is to replace (fc^ + 1) by 
gZ ^ . Let us define, from formulas (1345b and ( 13501 ) 



f{y) [j{k,K) + J-\k,K)\^ 



(449) 



5Z[X)^Y-^ 



' dk^k''-^f[l/g,) 



(450) 



For the choice g^7^ = [k"^ + 1)'*'''^ on which we focus from 
now on, the calculation can be brought in a form very similar 
to the case 7 = 1 as follows: 



5Z{\) 



dPk 



_}_dSd 

e/2 4 

e/2 27 7i 



(451) 



Taking the integral to the critical dimension d = due, and 
using that due = 27 and that lime-i.o e/2 = Sd^^ for any 7, we 
arrive at 



6Z{X) 



dyy^'^-'iy^'^-lY'^-'fiv) (452) 



The two cases of most interest are short-ranged elasticity (7 = 
2, dc = 4), and long-ranged elasticity of the contact-line or 
fracture front (7 = 1, dc = 2). For these cases, Eq. (14521 ) 
reduces (after a shift from y to x + 1) to 



5Z{\) 



5Z{\) 



SR 

d=i 
LR 

d=2 



dxxf{x + l) + 0{e) 



(453) 



da; {x + + 1) + 0{() (454) 



Hence the two calculations are very similar For short-ranged 
elasticity, the results where given above. For long-ranged elas- 
ticity (7 ~ 1, dc = 2), we have plotted the resulting functions 
for 5Z{\) and 5V{u) on figures[T5]and[T6] More details about 
the calculation and the results are presented in appendix |T] In 
particular we find that the exponent of the small-velocity be- 
havior changes to 



a = 1 + 2a + Oie^) = 1 - -(e - C) + ©(e^) . (455) 



V. FIRST-PRINCIPLE CALCULATION OF GENERATING 
FUNCTIONS TO ONE LOOP IN THE POSITION THEORY 

A. General framework 

Let us now go back to the more conventional formulation 
of pinned elastic systems formulated in the "position theory", 
i.e. Uxt rather than in the velocity variable iixt- 

Let us go back to the original equation of motion in the 
laboratory frame 

/ Klx't'U^'t' ^ f 9^l'Wx't + F{uxux) (456) 

J x' .t' J x' 

R^L't'^Sn'i6,,.'Vdt'-g-^,). (457) 
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FIG. 15: SZ{X) for LR elasticity (7 = 1, = 2). 




FIG. 16: iiSViu) for LR elasticity (7 = 1/2, dc = 2). The form of 
the curve is slightly different from the SR case, e.g. it crosses zero 
for ii slightly larger. 



We want to compute an arbitrary generating function in the 
position theory 



G[fi, w] = e/xt A'xtu, 



(458) 



It can be written as an expectation value with respect to the 
dynamical action S[u, u], 

(459) 

This dynamical path integral is normalized to unity, 
J 2?[u]2?[u]e^'^["'"l = 1. The dynamical action, now for the 
displacement u, and a different response field u instead of u 
is 



S — So + Sd 
5o = 



(460) 
(461) 



Sdis ^ -7; [ UxtUxt"^{uxt - Uxt') ■ (462) 



Note that here we have chosen to consider w as a source and 
not included it in S, although this is a matter of choice. (Dis- 
order independent) initial conditions are easily specified con- 
sidering the path integral with fixed endpoints and convolving 
with the normalized initial distribution P[{ux,t=t„}]- Non- 
zero temperature leads to an additional term — 77T / ^ u^^ in 
So. 

To obtain an exact formula for the observable G[fi, w], we 
need to consider the effective action r[M,u], associated to 
S[u,u], defined in the usual way as a Legendre transform. 



Uxtgxx''^-^'t+ / f^xtUxt ■ (463) 



T[u,u] + W[^i,w] 

Knowledge of F allows to obtain our observable as 

G[fi,w] = eJ'x*Mx.<r+L.Cr9:i,-x'*-rK-,«-n 



(464) 



in terms of the solutions u^j"' and u^j"' of the "exact" saddle- 



point equations 

ST , 



ST 



[U, U\ = fixt , 



OUxt OUxt 

These solutions are such that 



I g-l.wx't. (465) 

J x' 



- u.w 



[Uxt, 



fl.W 



Qxx' 



5W 

SWx't 



"^xt — \U'xt/fj.,w 



sw 
s^J-xt ' 



(466) 
(467) 



and thus u^^"" vanishes when = 0. There are other inter- 
esting properties. The covariance of the action under the STS 
transformationist Uxt + <Px, Wxt Wxt+gxx'4'x' implies 
that G[iJ,,w + g4>] = e'^'^G[fi,w], hence taking a derivative 
w.rt. w one finds the property 



— 1 ^u.-w 

9xx'<'t = / ^^xt 



(468) 



Note that because of the saddle point equation, in the deriva- 
tive 



J x' 



9xx'K:^ (469) 



one can differentiate only the explicit dependence on w. 

The effective action can be computed in a loop expansion 
as follows. Consider U :~ {u,u) a shorthand notation for the 
fields. Then for an action of the form 

S[U] = So[U] + Sdis[U] (470) 

the associated effective action can be computed as 

m = So[U] - ln(e~^-[^+*^l)i;Pi . (471) 

Here {■■■)^,^ indicates that averages over SU should be per- 
formed using the action So and that one keeps only graphs 
which are 1 -particle irreducible w.rt. the vertex S^is- Hence 
these diagrams are sums of 1-loop diagrams. 
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B. Tree calculation 

It is easy to see that, if one allows only for tree diagrams, 
one has 

r'™<=[C/] = To + So[U] + 5dis[C/] = To + S[U] , (472) 

since the only IPI tree diagram is the vertex itself. We have 
defined Fq = |tr hi Sq which is just a constant. 
This leads to the tree approximation of G[^,w], 



(473) 

where in this Section the u^j"', u'^'t" are solution of the saddle- 
point equation ( |465t with the replacement F — )• 5, and will 
also be denoted by u^^'"'''^'^", in the following. As is 

well known, this is the sum of all tree diagrams in perturbation 
theory of the non-linear part i.e. 5'dis. It leads to the mean-field 
theory, as discussed below. 

Note that because of the saddle-point equation, in the 
derivative 

Jt,x' 

one can differentiate only the explicit dependence on w. 
Choosing e.g. Wxt ~ vt one obtains 



vt] 



=0+ 



= / .9^^ 



(475) 



u) = 0] = 1, which is not neces- 
sarily true, except if the system is prepared in the Middleton 
state, which we now assume. 



using that A'{u) is an odd function. In the absence of disorder 
the solution is w = R^fi and u — C ■ fi + R- g^^ -w with C = 
2r7oTi?^i? and one checks that G'[Ai]"''='= = e'^fl g-'-^'+^M C M^ 
as expected. Taking a time derivative of the first equation, one 
notes the structure 

{R-^ + ±)-u = g-'^ - w (479) 
{R-^''^ + t)-u = fi (480) 

^xt,x't' ^ Sxx'^W Uxt"A'{Uxt ~ u^t") (481) 



The scalar product "•" denotes integration over the common 
space and time arguments. We can now compute ( 14731 ) by 
substituting the solution of ( 1476b : again using ( 1476b it can be 
simplified into the two equivalent forms 



^fl-U-— ^U-B. ^ -u— ^u-g ^ -W 



(482) 



2. Expansion at small driving w — 0'^ 

The solution of the above saddle-point equations can be ex- 
panded in powers of Wxt, assuming that fxt ~ (tj^ — Wx)wxt 
is a monotonous function of time for each x. We find 



u^I = + u 



xt — "-x ^ ^xt ^ ■■■ 1 ^xt — '^xt ^ '^xt 



(483) 



From Middleton's theorem we know that we should look for a 
solution of the saddle-point equation such that w^J*" — u'^'^' > 
for t — t' > Q, hence should satisfy this property. 

a. Lowest order: At lowest order, i.e. Wxt = 0+, the 
first saddle-point equation leads, using (14781 ). to the quasi- 
static solution-^ 



1. Tree saddle-point equations 

Let us now specialize to = + ni^. To tree level we 
need to solve the following saddle-point equations: 



VodtUxt + {m^ - Vl){uxt - Wxt) 



Uxt'A{Uxt - Uxf) = 



(476) 



{~mdt " + m )Uxt - J UxtUxt' A' {Uxt - Uxf) 

= l^xt ■ (477) 

Its solution is called u^'*", u'''"' only when needed, otherwise 
u,u. At non-zero temperature there would be an additional 
term -2T]oTuxt on the rh.s. of Eq. ( 1476b . Note that u^-™ 
vanishes for ii = 0. We now consider sources fixt which 
vanish at t = ±oo, hence we also assume that Uxt vanishes at 
t = ±oo. Note also that Uxt Uxt + 4'{^)iWxt — ^ Wxt + 
(— + m^)(t>{x) is a symmetry of the equations (STS). We 
further have the remarkable property 



{m' -Vt) Uxt 



l^xt 



(478) 



<. = A(0)K-V2)-2 / 



(484) 



while the second saddle-point equation leads to the "instanton 
equation" for u, 

{-riodt -Vl+ m^)ult + au^t J ult,sgn{t - t') = ^ixt ■ 

(485) 

where here and below we denote 



a := -A'(0+) , 



(486) 



and we use 



A'{Uxt - Uxt') = -crsgn(i - t') + lS."{Q){Uxt - Uxt') 

+0{{uxt - Uxt')"") ■ (487) 



(ii) 



Note that we expect that there are other solutions corresponding to a non- 
steady state, e.g. solutions with other prescribed boundary conditions. 
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b. Next order: To first order in Wxt one finds 

Jx'.t' 



(488) 



A"(o) {r^)-^ + j: 



X / ul,,,ul,,^{uU, - ui,tj . (489) 



We have defined 



Stt' / sgn(i - - sgn(< - t'){i'^t 



We also used that u],^ = 0. Note that 



. 



(490) 
(491) 



(492) 



3. Case Hxt = and connection to the velocity theory 



In the velocity theory one is interested in observables ( 145 8t 
such that 



f^xt = 



fJ-xt 



'dtXx 



(493) 



where \xt vanishes sit = ±oo. Then Eq. ( 1478b implies that 

Uxt=^ , Uxt = -dtUxt , (494) 



where Uxt vanishes at i = ±oo. Note that at the level of the 
MSR action one can rewrite 



^xtRxlx-t'^^'t' = / ^xtRxlx't'^x't 



(495) 



The saddle-point equations in the velocity theory then read, 
after some integrations by part: 



{R-^ +t)-u^g-^ - 4] = f 
(i?-^-^ + t)du = dX 

'^xt,x't' = ^xx'^tt' 



(496) 
(497) 
(498) 



2aUxt+ J Uxt"Uxt"^reg(.'Uxt - Uxt") 



To lowest order in w, i.e. for w = 0+ we obtain 



which is exactly the instanton equation ( |9T1 ), recovered here 
from first principles. In section ITlIB I we have obtained it by 
neglecting higher derivatives than the first of A(u); we see 
here that the contribution of these derivatives indeed vanishes 
if one looks at tree diagrams for w ^ 0. They do not vanish 
however to higher orders in w, or at non-zero velocity. 

We now go beyond the tree calculation and consider one- 
loop corrections. 







-A, 



(499) 
(500) 



C. 1-loop calculation 

Now we compute r[U] by including all tree and one-loop 
diagrams. It is then easy to see that 



^tree+l-loopj^j = S[U] + [U] , (501) 

ri [U] = itr In 5" [U] - ^tr In [U] , (502) 

and we assume T^[U] to be small. Let us denote A := 
{g~^w, fi). The saddle-point equations are thus 



hence U = U^""" + 0{T^). To compute 

Q ^ ^Au-s[u]-r'[u] 



(503) 
(504) 



(505) 



we can thus consider as an explicit perturbation and to the 
same accuracy, i.e. neglecting terms of order (F^)^, 

G = G"-'=e-i^'[^"°°l . (506) 
Going back to our explicit notations, we thus need to compute 

G[/i, w] = G*>--[/., ^]g-rM^-— . (507) 

D. Explicit calculation 

From now on, we focus on velocity observables, i.e. the 
case 



f^xt — , fJ,xt — ~dtXxt 



(508) 



for which (14931 ) and ( 14941 ) hold, and will be used extensively 
below. One thus has that Fi = for w = 0^. In this section 
U = {u,u) denotes = m^j'"''''™), and all 

derivatives are taken at the tree saddle point. 

To compute Z{X) we need to expand to first order in if. 
The small-w dependence of denoted U here, can be 

obtained from ( 1503b : 



(509) 
(510) 



We need to compute 



1 



Fi = -tr(ln5"[C/]) -tr(lni?-i) 



2 
1 

liS")-pS:p,iS")-}-g-'-w + 0{w') . (511) 



For now, we ignore the quadratic substraction. Greek indices 
denote either u or w and all contractions are implicit. 



48 



At this stage this is still general enough to treat a non- 
uniform \xt- However for simplicity we will now focus on 
the case of a uniform X^t ~ At, i.e. on center-of-mass ob- 
servables. The saddle-point solution is then uniform and we 
denote Uxt = = —dtUt- is independent of w. 

We need first S", the matrix of second derivatives. It is 
computed in Appendix U for general U, then specified for 
^tico ggjjgj-^j ^ Here we need it only in the case (15081 ). 
and for a uniform A, hence we can use ii^ = and it sim- 
plifies further into 



S,~,,, — R 



T\-l 



(512) 
(513) 
(514) 
(515) 



The "self-energy" S is defined in ( 14901 ), and reads 

'^xt,x't' — ^xx'^W (516) 

Stt' = (y[5tt' I sgn(i - - sgn(f - t')ul] 

= -a[2<5um°-sgn(t-i')at'4] • 

Note that 



(517) 



^tf = 



(518) 



The first component is actually A(0), but can never appear 
for velocity observables; hence we dropped it. To pursue, we 
define the dressed response 



7^ = (i?"^ + E)-i 



In Fourier 



{Tlk)t 



kW 



(519) 



(520) 



with (i?^, ^)tt' = Rktf '^^is dressed response is related to the 
one defined in ( |269t and ( 1327b , 

^kw 9{t - t') e-('='+i)(*-*')+2/.' ^521) 

Namely one has 

nut' ~ idt)-^Rktt'dt' , (522) 

where the ss means that it is true up to a zero mode. The 
correct identity, proven in Appendix iKl reads 

nktt'fpt' = J {dty^Rkwdf {(j)t' - 0-oo) + \ 

(523) 

upon acting on a test function (j)t . This implies the following 
property 



dtTlkt 



^kW Ot 



(524) 



used extensively below. The above relations arise because we 
are working in the position theory in a case where we compute 
velocity observables. 

We now need the inverse second-derivative matrix. One can 
first invert the 2x2 block structure 



(525) 



where the inversions on the rh.s. refers only to the space and 
time dependence. Given that in addition one has S'^^^ = 0, 
there are only three distinct terms in the sum (151 lb of order 
0{w), and which we denote 



/'c"\ — Ic'" /'c"\ — 1 — 1 



(526) 



We now specify to a uniform w^t = Wt ■ The third derivative 
tensor is computed in Appendix |j2l It is important to note 
that iS"„ and all components of S'" are local in space, i.e. 
^xt x't' XI ti = 5xx'xiS^[, j^. We can then make the momen- 
tum structure more explicit, using the above second-derivative 
matrix, and write 



5Ti ^ 5V\ 



(1) 



5T 



(2) 



5T 



(3) 



(527) 



Jk 

= I [nl] \S^uu ■ nl ■ 51 • 7^o • H . 

Jk 

All three terms are matrices in the time indices only, i.e. 

[5l]^^,=A"(0)u?4, (528) 

and [SZuh'u = [S'iUtUt'- Note that we can define 



n-ott'Wt' 



(529) 



and replace it above since it appears on the right in all three 
terms ( 1527b . 

We now specify to the choice of most interest for us here, 
namely the driving at small but finite constant velocity Wf = 
vt' . In that case TZt is not a well-behaved expression, since it 
may contain an additive term in the position of the parabola. 
Fortunately, in the calculation below, using dSM only the fol- 
lowing combination will appear: 



dtUt 



If Jt'<t 
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In particular, 



lim dfTZt = V 



R 



0,4*2 



(531) 



t2<t 



It is shown in Appendix that the third term vanishes, 



(3) 







(532) 



Hence we only need to compute two contributions. 

Substituting ( |528t into ( 15271 ) we compute the first term. 



^ ■'k,t,t',ti,t 



1 



^ Jk,t,t' ,ti,t2 



Jk,t,ti,t2 
Jk,t'<t<Q 



(533) 



To obtain the second line we have integrated by part over ti 
and t2 and used (15241 ). No boundary terms are generated since 
Mt vanishes at i ~ ±oo. We used that 9t9f' (7?.t —7?.t')sgn(t — 
t') = -dtdt'Tlt'Sgnit - t') = -2dt''Jlt'S{t - t'), i.e. there is 
a factor of 2, not 4. This is the first term obtained in Eq. (13261 ). 
Graphically, this can be written as 



Wt5 




(534) 



Only the first term is non-zero. 
For 6T'f'\ we find 



(535) 



Jk,t,t' 

= -m2A"(0) / Uktfdtu'linv - Tit) 

Jkd,t' 
Jk,t,t' 
Jk,t,t' 

We have used that TZtt ~ 0, absence of boundary terms 
[TZku'u"tiTZt' - 7^t)]^Z^^ = and dtUt. = dfUt = 0. 
We have also employed (15241 ). 



Now we can use ( 15231 ) with = 1 which gives J^, TZkw 
J^, Rf.^^, = (1 + fc2)-i and obtain, using ( l530l l and ( |53T] ) 



k Jt2<t'<0 
1 



fc fc^ + 1 Jt2<t 



lot't2'^(k,t') 



t ^0**2 



(536) 



The first term is exactly the term proportional to the single 
<i>(fc, t) in our previous calculation (1326b . The last term can be 
calculated, recalling the definition k : 



A 
1-A 



dtu°Rau2 



m'vA"{0) 



Graphically, this can be written as 



(537) 



5rf ' = 



We can now put all together and obtain 




(538) 



(539) 



Z{X) = Z*'''='=(A) - lim — , 

which coincides with the result (13321 ). (13301 ) apart from the 
additional contribution A n ( 1 + k'^)"^. This contribution, 
equivalent to ( 14391 ) and (14401 ). exactly cancels the 0{k) in 
Z{X) to one loop, as it should and automatically removes the 
quadratic divergence. It is thus exactly the quadratic counter- 
term. While in Section lTV F 21 it came via some manipulations 
on the seemingly vanishing term A'" (v{t—t')+Uxt — Uxt'), in 
the present calculation it appears automatically, and is related 
to the zero mode of the velocity theory. 



VI. CONCLUSION 

In this article we presented in detail the novel tools and 
methods which allow to calculate the statistics of velocities in 
an avalanche for the prototypical model of an elastic interface 
driven in a random environment. It is the extension to the dy- 
namics of our work on static avalanches, and the quasi-statics 
reveals to be closely connected, albeit different, from the stat- 
ics. The dynamical observables are much richer as we aim to 
calculate many-time correlations. The problem of how to de- 
fine an avalanche, and the steady state measure for avalanche 
statistics, is addressed and allows to make progress. At the 
same time connections to avalanches following a kick, or non- 
stationary avalanches are discussed. The Middleton theorem, 
which allows to order all realized configurations in time, plays 
a crucial role at all stages of the derivation. 

Our construction starts by identifying the correct mean- 
field theory, valid in space dimensions d > d^c- We discover 
that it is given, up to renormalization of a few parameters, by 
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a simple tree theory, itself equivalent to a non-linear instanton 
equation. This tree theory is interesting in itself. For the cen- 
ter of mass of the interface it exactly reproduces the ABBM 
model; it settles an important question concerning the validity 
of the ABBM model, introduced before as a toy model. The 
full space-time statistics of the velocity field is found to be 
given by the Brownian force model (BFM). This model is ex- 
actly solvable, reducing the problem to solving a space-time- 
dependent instanton equation. Our methods allow to obtain 
a host of new results for the probability distributions at sev- 
eral times and a number of results at non-zero wave vector q, 
which go beyond the ABBM model. A salient result is the 
time asymmetry of the avalanche shape, which, within mean 
field, manifests itself at the local level (or at non-zero q) but 
not for the center of mass. The universality of our results is 
discussed and quantified. 

Continuing to 1-loop order, we obtain the distribution of 
instantaneous velocities in an avalanche for an elastic mani- 
fold, as e.g. a magnetic domain wall, driven through disorder 
These results are new, and have never been addressed before. 
They are the basis for further work on the avalanche duration 
and shape, beyond mean-field theory. 

Many of the results of the present article can be confronted 
to experiments, and for this purpose we have extended them 
to long-range elastic kernels which are ubiquitous in nature. 
There are numerous experimental systems at their upper crit- 
ical dimensions (e.g. magnets) and non-zero q observables 
have not been measured and discussed previously. For other 
classes of systems below their upper critical dimension, the 
techniques introduced here provide a novel and at present the 
only way to attack them. 

Let us list a few important prospects for the future. Since 
we now know how to describe the space-time structure of 
avalanches within the mean field theory, using the Brown- 
ian force model, it would be interesting to develop analyti- 
cal and numerical techniques to solve its evolution, and solve 
the space-time dependent instanton equation beyond what has 
been done here. This should yield a detailed description of 
the space-time processes involved in an avalanche, and shed 
light on their physics. Avalanches have similarities as well as 
differences with branching processes, and the spatial shape of 
an avalanche is an important observable for experiments. We 
have voluntarily focused on the small driving-velocity limit, 
since at present the FRG is better controlled in that limit, 
but an important challenge is to understand the finite-v be- 
havior, and in particular whether the w-dependent avalanche 
exponents present in the ABBM model survive beyond mean- 
field theory. Other more far-reaching issues are to treat non- 
monotonous driving, hysteresis and to extend the theory for 
systems which do not obey in an obvious way Middleton's 
theorem. 
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Appendix A: Laplace inversion for a time window 



We give here the inverse Laplace transform (1211b in a series 
representation. By inspection we find that for any finite T the 
LT has simple poles on the negative real axis at s = s„ < 
— 1/4, n = 1, ... the closest one to zero crosses over from 
si(T) = -l/T at small T to si = -1/4 at large T. Since all 
s„ < —1/4 we can write s = — . Noting x = tant/j the 
poles are solutions of — V'n = -j tan ipn — 7i7r/2. The function 
s„(r) is better represented as a function of s„, 

4[^ - arctan(V-l - 4s„)] , ^ 

T = ^ , ^ ^ , ^ ^ s = Sn{T) (Al) 

represented in Figure [TT] Now we can compute the residues 
and using the equation satisfied by the poles we find, amaz- 
ingly, that they are all simply all equal to 1/T. Hence 



oo 

^ ' TU ^ 



-\Sr.(T)\U 



The small-T behavior of the poles is 

ki(T)l = ^ + ^ + o(r), 



\Sn{T)\ 



+0{T) . 



1 



T V4 7r2(?i+l)2 



(A2) 



(A3) 



(A4) 



We are grateful to Alexander Dobrinevski for numerous 
useful remarks. We thank Andrei Fedorenko, Alejandro 
Kolton and Alberto Rosso for stimulating discussions. This 



Hence at small T we get 



'U/T 



TU 



(A5) 
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consistent with the velocity distribution, as discussed in the 
text. For large T the poles behave as 



K{T)\ 



T2 



T3 



rp4 



o 



1 

2^5 



To leading order at large T one can keep only the first two 
terms, and approximate the sum by an integral, which repro- 
duces the correct asymptotic result 



dn ( 



1 



(A6) 

equal to the avalanche-size distribution as discussed in the 
main text. 



Let us thus discuss here the dynamical part of F, and con- 
sider the dynamical coefficients D„, as examples of irrelevant 
operators. For concreteness we restrict to SR elasticity with 
due = 4. The perturbative correction to the inverse response 
function reads, to lowest order in the disorder (see e.g. [il02il ) 



dt 



jqJa V 
which leads in the limit of tj 0+ to 



-(1 -e~'"')A"(0+) , 
(B4) 



= -77/2A"(0+) , 
SD^ = (-l)"77"7„+iA"(0) 

rA 



In 



(B5) 
(B6) 

(B7) 



Appendix B: Irrelevant operators and response function 

The effective action of the position theory in the laboratory 
frame can be written in an expansion in powers of the response 
field M as 

°° 1 f 

p=l P- •'^iM 

The term p ~ 1 expanded to linear order, T^^^ [u] ~ 
TZ^^^, t^t'^x't' + 0{v?) defines the exact inverse response 
function. Expanding the latter in time derivatives defines the 
renormalized dynamical parameters, more conveniently ex- 
pressed in the frequency domain. 



A is an UV cutoff. For d < 6, to which we restrict, we have 

limA^^oo In = m'*-2nj^ ^jtjj ^ J^(;,2 ^ lyn.^ ^gjj. 

defined for n > 3. We define /a = {AttY/^T{2 - |) as the 
analytical continuation to any d, with I2 = m''^^/2 for d < 4; 
the integral I2 becomes UV divergent for d > 4. 

One now defines the dimensionless inverse response func- 
tion at (7 = 0, with times scaled using the characteristic time 



00 

n=2 



(B8) 
(B9) 

(BIO) 



T^qlo,^ + Vi^ + Y. Dniiojr ■ (B2) 



n=2 



Similarly, in the limit v = 0+ the local time-persistent part of 
the term p = 2 defines the renormalized second cumulant of 
the disorder A(m), 



1™ / ^ii.t,ii^,t, [{uxt = ut}] = L''A{ut - Uf) 



(B3) 



Similar definitions hold for the p-th disorder cumulant C'(p) 
from the term or order p in T. All renormalized quantities 
depend implicitly on m. 

The main point is that near d = due and in the limit m ^ 0, 
the only relevant terms, i.e. operators in F are and A(u), 
irrespective of the details of the bare model. For d = due ~ 
e > 0, all other pieces of F are irrelevant, i.e. higher orders in 
e. For d = d„f. they are higher powers in 1/ ln(A/m). 

In Refs. SItIIIII, this property was discussed in detail 
for the disorder-part of F, for instance that the dimensionless 
(i.e. rescaled by the appropriate power of m) higher cumulants 
(JP = 0{eP) forp > 3, and similarly, that the non-local part of 
the second disorder cumulant is 0{e^). Since the local second 
cumulant A = 0(e), it implies that the renormalized disorder 
V ~ 0{y/e) is local and gaussian, and that all other disorder 
operators are irrelevant. 



The Dn are dimensionless. Let us now discuss the two rele- 
vant cases: 

(i) d < 4: 

Using -ma„/„+i = (2n + 2 - d)/„+im''-2"-2, 
Eqs. (iBSt-dBT]) lead to the RG equation, up to 0{e'^), 



- mdm7] = -riA"{0+) 

-mdrnDn = -2{n - l)bn 



(BID 



+ (-l)"(2n + 2-d)^A"(0+) , (B12) 



using the rescaled disorder (1281 ). Since for d < 4 the behavior 
of A(7i) is universal for small to, so are the behavior of r] and 
of the coefficients I)„. The first equation gives 77 ^ rn?~'', 
i.e. Tjn ^ with 

^ = 2-A*"(0+) = 2-i— ^e + 0(e2). (B13) 

The exponent z is the dynamical exponent, with z < 2. In the 
second equation we can use iFtsIi 

4+1 ^ F(n + 1 - d/2) d^4 1 
e/2 2F(n+ l)F(3-d/2) 2n(n - 1) ' 
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Hence for d < d^c the scaled dynamical coefficients converge 
as m ^ to universal fixed point values, given to lowest order 
in e by 



(-1)" i-^r l^.^oie^) 



-A*"(0+) = 



2n(n-l) ' ' 2n(n-l) 3 
For d = due = 4 one finds analogously 



(B15) 



(-1)" A*"(0+) (-l)"47r2 1-Ci 



2n(n- 1) ln(A/TO) n(n - 1) 3 ln(A/TO) 



(B16) 

Hence at the upper critical dimension the dimensionless coef- 
ficients Dn decay to zero at small m, thus the model is faith- 
fully described by the BFM and ABBM mean-field equations 
of motion, with (only two) parameters, r,„ and cr,„ . The be- 
havior is universal, and largely independent of details of the 
bare model. For d < d^c the model is not described by mean- 
field theory, but by a new universal fixed point which is stud- 
ied in SectionHV] We can obtain the inverse response function 
for d = due — e by inserting ( IB15t into (IB 8b and summing 
over n > 2, 

fiy) = l + y+ ^{{y + 1) ln{y + 1) - 2/)A*"(0+) . (B17) 

Thus the final result for the inverse response function to one 
loop, i.e. 0(e) accuracy is 



1 + iu;T„,-y + Oie^) . (B18) 



We used the result (IB13I ) for the dynamical exponent z. The 
behavior for large ujTm ^ 1, i.e. in the limit of small mass m, 
is TZqlj ^ (ioj)^/- as expected from scaling. This provides a 
derivation of the dynamical exponent at finite frequency. 

(ii) d > 4 

The F RG flow of the disorder for this case was discussed 
in ifiol (Appendix H) and III (Appendix B). There are two 
phases: (a) if the (smooth) bare disorder is small, it remains 
smooth under coarse graining, i.e. there is no metastability, 
no cusp, and no avalanches, (b) if the bare disorder is larger 
than a threshold, A(u) acquires a cusp, but flows back to 
zero as A(m) ^ (^)'^^^v4(m), where A{u) is non-universal, 
and equivalently, A(u) ^ A'^^'^A{u) is non-universal. Al- 
ternatively, if one considers a non-smooth and weak bare 
disorder (i.e. with a cusp in Aq{u)), then perturbation the- 
ory converges, schematically A — Aq ^ /20(Aq) where 
I2 ^ A'^~'^ — ni'^~^ since I2 is now UV convergent and dom- 
inated by the UV cutoff (see iFtsIi for details). 

Since the rescaled disorder A flows to zero as m — > 0, 
the FRG equations dBl II ) shows that 77 converges to a non- 
zero value as m — > 0, hence z = 2. The value of 
jyp, = rjo exp(- ^A;^(0+)) obtained from (IbTTT ) is non- 
universal, since the flow of the disorder is itself non-universal. 
The coefficients Dn, on the other hand, using dBl II) converge 
to zero as 



Dn 



i-iy 



2n{n- 1) 



- (-) A"(0+) ; (B19) 



hence for m ^ the model is well described by the ABBM 
model with constant but non-universal parameters ?/ and a. 



Appendix C: A differential equation for Z(A) 

We give a very general argument of how to calculate Z{X), 
without calculating the instanton. This method works for all 
first-order instanton equations. 

The instanton equation away from the source reads 



dtiiit) = uit) - u{tf =: f{u{i)) 



(CI) 



where we have allowed for a possible generalization to an ar- 
bitrary function /(u). To obtain Z{\), one has to integrate its 
solution 



Z{\) 
u{t{X)) 



t(A) 



dtu{t) 



A 



(C2) 
(C3) 



Note that the translational zero-mode in time of u{t) is not 
fixed in (ICll ). but by the condition ( IC3b . Compared to the 
standard solution, there is an arbitrary change in the time of 
measurement. Taking a derivative w.rt. A of the last two equa- 
tions yields 



dZ{\) 
dA 

dtu{t) 



t=t{\) dA 
Combining these two equations yields 



dZ(A) u{t) 



dA 



dtu{t) 



uit) 



t=t(A) 



(C4) 
(C5) 

(C6) 



t=t(A) 



where in the last step we used the instanton equation ( ICll ). 
Using (IC3b we find the simple result 



dZjX) 
dA 



A 



/(A) 



(C7) 



If /(u) ^ u — u^, the case usually considered, we arrive at 



dZ(A) 
dA 



1 



1-A 



The solution is 



Z(A) = - ln(l - A) 



(C8) 



(C9) 



where the integration constant is fixed by demanding that 

Z(0) = 0. 



Appendix D: More details on the ABBM model 

In this appendix we use dimensionless units. Let us rewrite 
Eq. d233]) as 

dtQ = -d,J (Dl) 
J(v,t)=-(9v(vg)-(v-«)Q), (D2) 
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where J(v, t) is the current of probability. The equation for 
the eigen-modes is 

-sQ = d,[d,{yQ) + (v - v)Q) . (D3) 

Let us first discuss the case w > 0. The general solution is 

Qiy) - v-'-ie'^ [CiL;1+^(v) + C2U{-s,v,sf)\ , (D4) 

given in terms of the Laguerre polynomials and confluent hy- 
pergeometric functions. The Laguerre polynomials can only 
have s = n = 0, 1, 2, .. since for different values they do 
not decay fast enough at infinity For these integer values of s 
the two solutions become linearly dependent. These Laguerre 
solutions for all s = n have the peculiarity that the current 
vanishes at the origin, i.e J(v = 0, <) = 0, more precisely 
J(v = 0, t) ~ v" at small v for all rt > L In addition the cur- 
rent vanishes everywhere for rt = 0. For the hypergeometric 
solution the current is J(v = 0, t) = y^^z^j- 

In their work 10 lU ABBM retained the solution with zero 
current at the origin, hence the solution which vanishes for 
V — > oo, 

g„(v) =v''-ie-^L,r""i(v) , s = n = 0,1,2,... . 

(D5) 

They thus obtained the normalized propagator |j2l|3|], 

^ Tiv + n) 

a formula valid for u > 0. The term n = is (5o(v) = 
y^~^e~^ /T{y) and integrates over v > to unity, the oth- 
ers to zero. Hence dv (5(v, vi, = 1. Since the cur- 
rent vanishes at the origin for all times (i.e the total proba- 
bility for V > remains unity), for large times the probability 
reaches the stationary state which has zero current everywhere 
g(v,vi,t) Qo(v). 

Let us now consider w = 0+. One then finds that (i) the La- 
guerre polynomials must again be of integer order to behave 
well at infinity (one has L(^^(v) = 1, Lj^^(v) = — v, and so 
on), (ii) The Laguerre solution corresponding to n = be- 
haves as e~"/t;, hence is not normalizable. (iii) The Laguerre 
solutions for rt = 1, 2, .. have a non-zero current at the origin. 
(iv) The hypergeometric solution does not behave well at the 
origin ^ \/v unless s is positive and integer, in which case it 
again becomes identical to the Laguerre solutions. The only 
possible solution for the propagator thus seems to be 

oo 

g„=o(v,vi,t) = v-ie-^^ni-i(v)L-i(vi)e-"* , (D7) 

ri=l 

which is the limit of ( |D6b for w = 0+, where the term n = 
has dropped because its prefactor l/r(u) vanishes. 

On the other hand, inspired by our result from the text, we 
found that there is another expression for the propagator at 
w = 0+, namely 

Q(v2,vi,t) = Q(v2,Vi,t)+5(v2)e-^^ (D8) 
(5(v2,vi,i) = Vie 1 e - . 



We recall z ~l—e * and that Q satisfies (IDll l with as v ^ Vi, 

Q(v,vi,t)p.vi ' ' p.J(v2-vi). (D9) 

V47ri(viV2)3'^ 

We have checked with Mathematica that JJ^ dv Q{y, Vi , t) = 
1, and that the (5-function piece in ( 1236b is crucial for this prob- 
ability conservation. 

It turns out that the two expressions (|D7| i and ( |236b co- 
incide for V > 0, i.e. Q(v,Vi,t) — Qv=o{^,^i,t) as we 
have checked numerically with excellent accuracy (the con- 
vergence of the sum over n is very good). However the 6 
function in ( 1236b is not reproduced. Hence the terms n > 1 
now have a finite integral over v. This integral does not add 
up to L Somehow the n = term is replaced, for v = by 

a delta function, multiplied by the factor e "'-i . This factor 
takes into consideration the absorption at zero, which is now 
present. 

Other boundary conditions at v > 0, such as absorbing 
ones, can be studied, which we leave for the future. 



Appendix E: Checks of the 3-time formula for MF (ABBM) 

We now want to check the 3-times correlation. We use the 
formula 

poo 

/ dve-"Ii{2ay^)Ii{2bV^) ^ Ii{2ab)e'''+''\ (El) 
Jo 

which yields: 

= . / — e e = ^ ^ 

V "3 z" 7 

X h {2 ^^^/' ./^^ . (E2) 

Herer7 = i + ^-1,7 = 7-A2 and 1-z" = (l-z)(l-z'). 
For A2 = we find 

J dv2Q{v3,V2,z')Q{v2,vi,z) = Q{v3,vi,z") , (E3) 

as expected for a propagator. Other useful identities are 

dv3e^'^'Q{v3,V2,z') = e^'^^-^') J^e"'^^!™^ - 1^ 

(E4) 

roo p — vi 1 / i_ \ 

/ di;ie^i"iQ(w2,?;i,z) = — e"^ (e^^Mi-Ai, _ 1] 

Jo Vi V2 \ > 

(E5) 
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This allows to obtain 



„AiUi+A2t>2 + A3t'3 



X Q{V3,V2,Z')Q{V2,VI,Z)- 



fl ,V2 ,V3 >0 

= ln(l — Xiz" — \iz' + \\\2zz) 

+ ln(l - \2Z - Xsz" + X2X3ZZ') - ln(z" - X2ZZ') 

— In ^1 - Ai - A2 - A3 + A1A22; + A1A32;" + A2A3Z' 
- A1A2A3ZZ') . (E6) 

We recognize that the last logarithm is Z3. Taking the three 
derivatives dx^ d\,^ gets rid of the other terms, and shows 
that 

VlV2V3Q{v3,V2,z')Q{v2,Vi,Z = LTIJ^.^^,.9ai9a,9a3^3 ■ 

(E7) 

Since the latter expression is also the inverse LT of 
9i23^i^2W3'P(wi, V2, W3), and since neither function contains 
a S function, we obtain 



9123^ ("1:^2,^3) = Q{V3,V2,Z')Q{V2,V1,Z)- . (E8) 

Vl 



This shows that the 3-time velocity probability can be writ- 
ten as a product of 2-time propagators, i.e. the 3-time velocity 
probability at tree level (i.e. in the ABBM model) is Marko- 
vian. 



Appendix F: Spatial correlations in the tree theory 



Here we give further calculational details concerning Sec- 
tion nil HI in particular we work in the steady state to lowest 
order in v and in dimensionless units. The results are exact 
for the tree theory, i.e. the BFM in any d, or for SR disorder 
in the mean-field theory. For the 3-point function to first order 
in V, computed in the text, let us indicate the following inte- 
gral formula, useful to generate a series expansion in q (with 
h<t2< 0): 



1 - A 



1 + A(e*i-1) 1 + A(e*^-1) 



o2(l+g")t' 



[1 + A(e*'-1)]' 



(Fl) 



Let us now detail the calculation of the 4-time correlation function, from which we will also extract the avalanche shape in the 
stationary state. Consider the source Af = Ao(5(t — io) + X-i8(t — t^) for to < < ^2 < ^3- In dimensionless units, this gives 



1 



1 + i^e*3-t 



1 _ A3Aoe'o-(l-A3)(l-Ao)e'3 
^ (l + An)A3e*o+An(l-A3)et3 ^ 



e{t < to) 



The dressed response function has six sectors. We indicate only those needed: 

(1 - Ao)(l - A3) + Ao(l - A3)e'''-*o - AoAae*--*^' + (1 + Ao)A3e*'>-*^ 



(1 - Ao)(l - A3) + Ao(l - A3)e*<'-*o - AoA3e*o-t3 + (1 + Ao)A3et<--*3 

1 + A3(e*''"*-^ -1) 

(1 - Ao)(l - A3) + Ao(l - A3)e*»-*» - AoA3e*o-*3 + (1 + Ao)A3e*-'-*3 
1 + A3(e*''-*-^ -1) 



, ta < tb < to < ts 
, ta < to < tb < t^ 



l + A3(e*"-*3 -1) 



to < ta < tb < ts ■ 



We now use 



(l + A3(e/-*^-l))(l + Ao(e*"-*'-l)) 



t<t' 



^q=0,t',t 



for t' > to 



(1- Ao)(l- A3)-AoA3e*o-*3 
^ , Ao(l-A3)(e*'-*"-l) + e*'-*3A3(l + Ao) 

(l-Ao)(l-A3)-AoA3e*o-*3 f°'^^<^« 



We must split the integral into two parts. 



gAoL'^iioy 



\3L''U3 _ 2v 



t'<to<ti Jto<t'<t 



q,ti,t' q,t2,t' 



(F2) 

(F3) 

(F4) 
(F5) 



(F6) 



(F7) 
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The result is 



((A3 - 1) e*^ - A3e*0 ^ ((A3 - 1) 6*=' - Xse*-) ' 
((Ao-l)(A3-l)e*3-AoA3e*o)4 

^ 2^1 I J, +'^) e*oAoA3-e*3(Ao-l)(A3-l) 



9^ + 1 



3(-g^-i)(ti+t,)-3t3 ((^^ _ 1) _ 1) e«3 _ AoA3e*°) ■ 
(A3-l)3(2g4 + 392 + i) 



X I (Ao - 1) (2g2 + 1) [(_e2(.^+i)*o) (^3, 2 (q^ + l) ; 2?^ + 3; '^j^^ 



+e'i<''+'h ( 3, 2 + l) ; 2g2 + 3; ^^J^ ) ] 



e2(9^+i)to (3,2g2 + l;2(g2 + i) 



e*^-*='A3 
A3 - 1 



+2Ao {q^ + 1) 



2qh,+to+U 2F1 3,2q2 + l;2(g2 + l) 



e*»-*3A3 
A3-I 



(F8) 



We checked that for g = this expression yields c^Aaf^Ai ^4(Ao, Ai, A2, A3)|a2=Ai=o and that for Aq = it yields ( 12851 ). This 
expression is not invariant by time reversal i.e. by simultaneous changes to — t3,ii —t2, t2 —ti, —to Aq O A3. It 
is invariant however, at g = 0. The non-invariance by time reversal can already be seen on the 4-point function, taking dxg dx^ : 



-2d 



2g-(q^+2)T-(q2 + l)(2ti+3t2) 



(l + g2)(2 + g2)(l + 2g2)(3 + 2g2) 

X 2 (2q6 + + 13g2 + 6) e((39'+2)*i + ( V+3)t. + (9^ + l)T) + 4 (^2 ^ 3) q2g(?^+2)ti+2(,^ + l)t. + i 

+ 4 + 2) 52g(,^+l)ti + (2,-+3)t. + f _ 3 (2q2 ^ ^2g(9^ + l)(U+2t,) ^ + 7^2 ^ g-j ^(,^ + l)(3ti+2t.+T) 

This function is not symmetric by ti — > — i2 and f2 ^ ii- 

If we take the limit Aq, A3 — > —00 we obtain SugUqtiU-qt^Su^ which we do not reproduce here. One can check that the first 
hypergeometric term yields zero, although the limit is quite delicate. Taking —dtgdt^ and dividing by the duration distribution 
we find our final result; 



(F9) 



(e*3 - 1) (e*3 - e*i) (e*^ - e*^ ) e(-'''"i)*i-('''+i)*=-3*^ 
^ 2g4 + 3q2 + 1 

(1 - 2q2) e*i+*^ + (1 - 292) e*i+2*3 + (1 _ 2g2) e*^+2*3 + (2g2 _ 3) ^t,+u_+ts ^ ^3^2 + 1) 6^*3 
+ (292 + 1) e*i+*^ + (2g2 + 1) e*2+*-^ - (2q2 + 3) 6^*3 

(2g2 + 1) (-e^ii"+^h) 2F1 (3, 2 (g2 + l) ; 2q2 + 3; e*i-*3) + (2^2 + 1) (3, 2 (g2 + 1) ; 2^2 + 3; e-*3) 



2 (g2 + 1) e*^' (e(^«'+i)*i 2i^i (3, 2g2 + l; 2 (q2 + 1) ; e*i-*3-) _ ^Fi (3, 2q^ + 1; 2 (g2 + l) 



(FIO) 



where we have set to = for simplicity; the general case is obtained setting ti ti — to, i = 1, 2, 3. 

I 
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Appendix G: Behaviour of the 1-loop correction SZ{X) near 
A = 1 

Here we indicate how we extract the behavior of 5Z{X) near 
A = 1. We recall our result 



6Z{\) = ^a„K" 



(Gl) 



n=2 



with a„ given in Eq. ( 1352b . and repeated here 

_ (?i-3)(7i-2)2lii(n-2) 
°" " 2^12 

61n(2) - 2n(n + l)(lii(2) - 1) 
n2(ri + 1) 

{n - l)(n((n - 6)n + 2) + 6) ln(n - 1) 



rj2(n + 1) 



(n2 - 8n + 3) ln(n) 
^ 2{n+l) 
a2 = lim a„ = 1 — In 4 . 



(G2) 
(G3) 



From the relation (1 — A)(l — k) = 1, in order to get Z{X) 
in the limit of A — )- 1, which controls the tail of for 
M — > oo, we need this expression for k — oo. However, 
the series expansion has a convergence radius in k of only 1, 
equivalent to A < 1/2. A first thing one can do, is to re- 
express this series in A: 



5Z{\) = a„K" 



OO 

p=2 



nAf 



The formula for the coefficients Cp is 

p 

Cp = {p~ l)!y]a„ 



(-ir 



(G4) 



(G5) 



The convergence radius of 6Z{X), as a series of A, is 1. While 
this is useful for intermediate values of A, it does not allow to 
study the singularity for A — > 1. In order to analyze the latter, 
we now derive an expansion of 6Z{X) in powers of — We 
start with 



The contour starts at oo + iS, goes to 3 + iS passes left of 3 
and then goes to oo — iS, for any < S < 1, see figure [18] 
The formula uses the residue theorem, and that the residue of 
— r at integer n is (—1)". Two remarks are in order: a„ 

sm(7rn) fe V y 

has three different branch-cut singularities, starting at n = 2, 
n = 1, and n = 0, and going to rt = — oo. Singling out 
the term 02 avoids crossing the branch-cut starting at n = 2, 
which would not be allowed. Second, one could try to move 
the explicit factor of (—1)" into (— k)". This does not work, 
for two reasons: First of all, tt/ sin(n7r), when prolonged to 
the complex plane, converges exponentially fast. This would 
not be the case for tt cot(n7r), to be used to produce the non- 
alternating sign. Worse, k", for negative k, when prolonged to 
the complex plane, actually diverges in the lower half-plane. 
This is why we use the formula as is. 

Having an integral representation for dZ{X), we can now 
prolong analytically for k — 00, by deforming the contour 
of integration to C2, which starts at —00 + iS, goes to 2 + iS, 
then passes at the right of 2, and finally goes from 2 — id to 
—00 — iS; see again figure [TSl This gives 



SZ{X) 



dn 
27ri 



sin(7r7i) 



a2 

n-2 



-\2 



(G7) 

Note that while the integral representation (IG6b is convergent 
for — 1 < K < 0, the representation (IG7b is valid for — 00 < 
K < — 1; the smaller k, the better the convergence. We have 
checked the integral representation (IG7b for k = —8, i.e. A = 
8/9 numerically. Then both ( IG7I ) and the A-series (IG4b give 
SZ(8/9) = 8.17538, with a relative error of 10"^. Therefore 
trusting our integral representation, we can now analyze it for 
large negative k. Then it will be dominated by the contribution 
at the beginning of the cut singularity of a„, which starts at 
n = 2, see the first term of (IG2I ). and the corresponding plot 
[Ts] Therefore for large negative k, the integral (IG7b is given 
by 



SZ{X) 



dn -(ri-2)ln(n-2) 



2ni 



O 



[ln(-K)]3 



(G8) 



<5Z(A)=^a„(-l)"(-^)" 

n=2 

dn TT 



2Tii sin(7r7i) 



a^i-KY (G6) 
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FIG. 18: The complex-n plane with the contours Ci of Eq. iG6t and 
C2 of Eq. \G1\ . The branch cut starting at n = 2 is indicated. 



One can obtain more subleading terms by expanding a„ to 
higher powers in (n — 2). Doing this, we find 



dZ{X) = 



1 



1 



21 + 27r^ 



[H-k)Y 2[\i-i{-K)f 16[ln(-K)]4 
15 + 47r2 585 + 210tt^ + IAtt^ 



4[ln(-K)]5 
757r2 



45- 



48[ln(-K)]6 
1 

" 2 ) M-n)Y ^ 



(G9) 



We can test this series against the integral (IG8I ): We find for 
K = -10^° that 5Z = 2.887 x 10^^ with a relative error of 
10-4. For K = -10^"", we find 5Z ^ 2.400 x lO^^"*, with 
a relative error of 10^^. For k — —10^™°, we find 6Z = 
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2.361 X 10^^^^, with a relative error of 10~^ (probably from 
the numerical integration). 

Expressed in terms of A, our final result is given in Eq. (I397t 
of the main text. 



Appendix H: An alternative approach to express the 1-loop 
contributions SZ{X) and SP(u). 

Here we calculate the 1-loop correction SZ{X) by first in- 
tegrating over momentum. More precisely we start from Eq. 
( 1330b . calculate $(fc,t) as given in Eq. ( 1340b . but instead of 
Eq. ( 1341b and ( 13421 ) we first integrate over t, and then k, leav- 
ing the <i-integral for the end. In order to be able to per- 
form the fc-integration, we have to introduce counter-terms 
right away. The term involving ^{k,t), with the necessary 
counter-term J^l\k, k, ti) becomes 



(fc2)dfc2 



J^l\k,K,ti)+ I dt$(fc,t)]R(fc,ti) 



K2e2*iEi(-ti) 



KEi (ti) 



Ke*i (k + Ke*i (2ti - 1) - h) 



J^P{k,K,h) 

Ke*i (2Ke*i - 1) 



1) 



(fc2 + 1)2 (^e*i - 1) + 



1 



fc2 + l 



(HI) 



(H2) 



The second term involving <l>(fc,<)2, with the necessary 



counter-term [k, k, ti) becomes 



(/fc2)d/fc2 



j!:p{k,K,ti)+ I dt^{k,t)'''R{k,ti) 



(^,>_l)3L^e^Ei(-2ti)(A^ + 2A^ti-2) 
- 2Ke'iEi {-ti) (Ke'i + kc^Hi - 6*^ - l) 



k2 - 2K2e*i + 2Ke*i - In (ti) - 7e - 1 + ln(2) 



J^^'{k,K,h) 



(fc2 + l)2(Ke*i -1) 



(H3) 
(H4) 



Several checks are in order: First, the two counter-terms, 
when integrated over ti reproduce the one given earlier in 
Eq. ( I3501) . 



ti<0 



ji'\k,K,tl)+J^^\k,K,h) 

K(3 + fc2) + 21n(l-K) 



(fc2 + l)2 ■ ^"^^ 

Second, both J^^^K,tl) and J^'^\K,ti) have a finite Umit 
for ti 0. This is why the last term in Eq. (IH2I) was added, 
even though the fc-integral would have been convergent with- 
out the term at fixed ti. 

We thus have found an integral-representation for 5Z{X) as 
defined in Eq. ( 13511 ), with the same counter-terms. 



5Z{X)= J^^\nM) + ■ (H6) 

Jti<0 



The two contributions were given in Eqs. (IHll) and ( IH3I ). 

We now note that all terms in Eq. (IH6b are algebraic func- 
tions of K, and thus of A. Hence the inverse-Laplace transform 
is possible. Replacing ti by t to alleviate the notations, this 
becomes 



5V{u 



f2{t) 



hit) 



Jt<o 



e-"/i(0 + e"i-'=* 



f2it) 



h{t)u , fi{t)u^ 



(e*-l)4 (e*-l)5 (c*-l)6_ 
= e\2t + 3)Ei(-i) - e\2t + l)Ei(-2t) + Ei(i) + e' ^2 - - e"' + i + : 



(2e* - 8e^* + 12e'^*) i + e* - 4e^* + Be"^ 



be 



Ei(-2t) 



- 2e' + 8e^' - 12e^' t - 3e^ 



lOe^ 



7e3* + 6e" 



Ei(-<) 



ln(V2) + 7E 



(2e 



2* ' e^n 



13e' - 126^* + 4e 



.3t 



3e* - 3e 



2t I „3t 



e^' - 1 



h{t)^ (8e 



2e^* ) t - e^* + 4e-^ 



e 4 

De 



Ei(-2t) 



2e^ 



-6 

4t 



) t + 2e2* - 2e^* + 6e''* Ei(-i) 



ln(</2) + 7E 



„3t 



o4t 



2e^ 



Ei(-2t) 



6e* - 9e2 



4e 



3t 



Ei(-t) + Y - ^ 



2t 



„3t 



„3t 



ln(t/2) + 7E 



(H7) 
(H8) 



(H9) 

(HIO) 
(Hll) 
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This is a closed expression for 5'P{u). We can now check all 
our statements made in the main text. First of all, we repro- 
duce the plot on figure [TT] 

For the small- it behavior, we remark that the integral (IH7l i 
is dominated by the terms proportional to c^^/i'^^'^ \ in the 
limit of small t. The leading contribution comes from expand- 
ing f2{t) for small t, and reads 



with 



t<o 



„M hit 
e t — 



_^\n(u)+j^^ (H12) 



Note that /^{t) and f4{t) could also contribute at the same 
order, but they have no term proportional to hit, thus they 
only correct the subleading term ^ l/ii leading to the final 
result 



, , lii(w) + 27E + T - In 2 
SV{u) = -2—^ ^ + 0{\nu). (H13) 



To obtain a systematic expansion one rescales t ^ iit and in- 
tegrates term by term in t the series expansion at small ii. This 
confirms the predictions given in Eq. (I388l l for the exponent 
a, and for the constant C in Eq. (|395t . 



Appendix I: Long-ranged elasticity 7 = 1 

In this appendix, we calculate all relevant quantities for LR- 
elasticity 7 = 1, dc = 2, with the kernel defined in the main 
text. We found in Eqs. ( 1447b and ( 14541 ) that 



LR/ 



(II) 



5Z'^'\\)= dx(a; + l)/(a- + l) + 0(e) , (12) 
Jo 



where f{x) is defined in the text, in other words, the calcula- 
tion is identical to the short-range case, except that when inte- 
grating over fc, we have to replace / d(/c^)fc^ by J d(fc^)(l + 
k^). This replacement can be performed before or after the 
time integral. 



a. First method 

In this method, we first integrate over t leading to formulas 
( 1345b and (13501 ): then we integrate over k with the modified 
measure. The series expansion is then given by 



(13) 



LR 



ln(2) 

2 (n2 + n - 6) hi(2) 



(14) 



">2~ n2(n + l) 

{n - 4)(n - 3)(n - 2) ln(n - 2) 
2^^2 

(2 - n)(6 + 2n - In^ + n^) \xi{n - 1) 
^ 7i2(n + l) 

(n- l)(?i2-9n + 2)ln(7i) 
2n(n + 1) ■ 
For n ^ 00, the leading behavior is 



_ -21ii(n)-f ~21n(2) 



0{n- 



(15) 



(16) 



Comparing with Eqs. ( 1376b and ( 13861 ) shows that 

,5Z^R(A) = -ln2(l-A) + ... forA->-cx). (17) 

Thus 

= -ln(l-A)[l + ahi(l-A) + ...] . (18) 

This is consistent with a modified critical behavior at small 
velocities, 

V]:^^ {ii) 4 ' a = 1 + 2a + 0(e2) . (19) 

The behavior for k —00 (i.e. A 1) now reads 



O 



,.2 



(110) 



41n(-K) V[ln(-K)]2 
This implies a different tail than in the SR case. 

b. Second method 



We find, analogously to Eqs. ( IHlb and ( IH3I ). the integral 
representation 



5ZLR(A)= / J^^\kM) ^ J^^K^^M) 
The contributing terms are 



(111) 



(112) 



X [2Ke*i (2tiEi(-ti) - 1) + 2k - 27Eti - 2ti In (-ti)] 

Ei (-ti) (-2K2e2*i - 4K2e2*iti + 4Ke*i + 4Ke2*i) 
+ Ei (-2ti) (k^c^*! + 4K2e2*iii - 4Ke2*i) + 2^^ 
- 4^26*1 + 2K^e^^^ - (7E + In (-ii)) (2Ke*i + l) 
+ ln(2) (2K + 2ti + 1)1 . (113) 
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Inverse-Laplace transforming Eq. (Ill II ) yields an integral representation for 'P'^^(m), 



(114) 



l-c* 



t<0 
t 



/2^«(^) , , frit) 



e\At + l)Ei(-2<) + Ae\t + l)Ei(-i) + 



(c*-l)4 (c*-l)5 (e*-l)6_ 

(e-*-l)(2e*t-2f + e*) 



(4e* - 16e2* + 24e^*)t + e* - 46^ 



3t 



-4e* + 16e^' - 246-^' t - 4e' 



+2e-* + 30e* - 32e- 



2t 



12e 

„2t 



3t 



|-6e 

12e 
3e* - 3e 



12e^ 



2e^ 



2t I „3t 



t 

Ei(-2t) 
126'*' Ei(-0 

i + (e* - 4e2* - 6e3*) [in(-i) + 7e] 



e^* + e^*) lii(2) + (46^* + 2e^*) nn(2) - 12 



(16e 
-12e* 



3t . 2t 



2t 



4e"')t - e 



46^* - 12e4* 



(46^* - 16e^')< + 2e^' + 4e''' + 126*' Ei(-t) 



„3t 



2t , . 3t 



„2t I Q„3n 



3t 



p3t \ 

^-2e4*jEi(-2t) 

„3t 



Ei(-2i) 
In(-t) + 7E 

2e'H + e^* + 2e4*) Ei(-0 + e* - 26^* + e^* 



„4t 



(115) 
(116) 



(117) 



+ (76^* + 26^*) ln(2) + (26^* + 4e'''*) t ln(2) - 2 (118) 



^jln(2)-|e^*ln(-i) 



e3*iln(2) 
I 



(119) 



The analysis of the small ii behavior gives 5V^^{u) 
_2 inu^ hence is consistent with the above result ( Il9l) . 



Appendix J: Second-order derivatives S" and third order 
derivatives 5"' 

1. Second-derivative matrix 

We give here the matrix of second derivatives of the action: 

5xx' -UxfStt' / Uxt^A"{Uxt - Uxti) 

(Jl) 



M = = 0, and ii. = iiP . Hence 



9" 



9" 



+ UxtUxt'A [Uxt - Uxt') 

6tt' ivodf +711'^) 

-Sxx' Stt' i ilxti^'iUxt - Uxti) 

L Jtl 

-A'{Uxt ~ Uxt')Uxt' 

Swi-Vodt'-yl+m^) 

-Sxx' Stt' / UxtiA'{Uxt - Uxti) 

+ A'{Uxt - Uxt')Uxt 
-Sxx'A{Uxt - Uxt') 



(J2) 



(J3) 
(J4) 



We will need it at the tree saddle point and to lowest order in 
w, i.e. for w — 0'^, where according to the previous section 



9" 



6xx'A"{0) 



{R ^ + T,)xt,x't' 



T\-l 



-Sxx'A{0) 



t^tt' I U^ti + UxtUxt' 
ii 



xt.x't' 



(J5) 



2. Third-derivative tensor 



In the text we need the third derivative tensor only at the 
tree saddle point with = 0+. It can be obtained from dMb 



MxtMa;'t'Mxiti Xiti 



5xx'A'{Q+){ul,^ul,,)Bgn{t~t') 



ti 



"xt'":c't'"xiti 2:1*1 



(J6) 



(J7) 



= <5,,,A"(0) Sw I ul,^{ul,-ulJ~ulA<t-uit,) 

(J8) 



'-'MxtMx't'^xiti 2:1*1 



4.'A'(0+) Stt' / ultsgn{t ~ t2) - sgn{t - t')u 



'\ -.1 

xt' 



*2 
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Consider now the uniform case fi^t = fJ-t and u1.j. = u^. 
T^^nS'^L't'.xuti = S^^'^iKt\U with: 

[S'lLh't, = <j{6tu " St'u )sgn(t - t') (J9) 
= -A"(0)(<5,,-,, / < 

[^'■a'uultt'ti = [S'i^'uuhtit' 



Appendix K: Dressed response functions for velocity 
observables in the position theory 



We note that with notation ( 



{R- 



)fct (pt. 



(Kl) 



Jt' 

Hence for a smooth function (ht : 



(K2) 



= dt<j}t + k^<Pt+<j}t+ J sgn(t-t')w°'/'t 

Hence apart from a zero-mode in time, 

iR-' + ^)tt' = (at)-'e'^' (a, + A;^ + l)e-^^' . 

(K3) 

The zero-mode can be treated as follows. Consider the con- 
stant vector (j)t = (f>^oo = const. Then because of ( 15181 ) one 
has 



J^{R-^ + E)u'(f>t' = {k^ + !)<!>- 



(K4) 



This implies that the vector (pt = (p-oo is an eigenvector of 
R~^ + S, with eigenvalue fc^ + 1. Hence one also has 



This yields 



1 



(K5) 



1 



fc2 + l 



(K6) 



Using the definition of Rfctf given in Eqs. ( 13271 ) and ( 1521b . we 
can rewrite Eq. (IK6b to get the fundamental equations 



(K7) 
(K8) 



nttt'^t' = / (at)-'Mfctt'at. {(Pf - 0-00) 
1 



fc2 + 1 

dtTlkw = ^ktt'df ■ 



The subtraction of 0_oo from cpf is noted for clarity reasons 
only. A first corollary is 

tr(ln(i?-i + E)tt-) =tr(lni?-i) . (K9) 

Similarly one finds that 



it' 



ft 

(KIO) 



Appendix L: Third diagram SFi 

(3) 

We now turn to the third contribution 5Tl 



(3) 



(3) 



Ti-kt'tS'- 



t4 



(Ll) 



Using that TZktt = 0, and exchanging t and t' we get 

Jtt't2t4 

= WfTA"(0) / at[7^fe«'Sgn(t-<')] 

"'tt't2t4t5 

X Rot2t"?2<IRot4t5 , (L2) 
where we have used (1524b and ( 15301 ). Now we use that 



dt[nktt'^gn{t-t')] 

= ^ktt'dt'Sgn{t -t')+ Uku'dtsgnit - t') 
^ ^2S{t-t')[Rktt' -Uku'] ■ 

Since TZktt = Kfett = we find 

st[^^ = . 

Graphically, this can be written as 



(L3) 



(L4) 




t' ' <2 U 



(L5) 



These terms are zero: The first term is the response at equal 
times. The second term, when viewed in standard diagram- 
matics can be mounted, moving one arrow-head from t2 to <4, 
or vice versa. So it is expected to be zero anyway. 
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Appendix M: 1-loop expansion for the lowest cumulants 

1. Expansion in A of ^( A) 

Let us first recall the result for the one loop contribution to 
Z{X) to all orders in k derived via perturbation of the instan- 
ton equation and displayed in Eq. ( 1345b . Here we reexpress it 
as a function of A and display it up to to order 4 in A, 



ZiX) 
Jik,X) 



1 



k 

2(3 + fc2) A2 



(Ml) 



1 



A + 



P (l + /c2)(2 + fc2) 2 
2(108 + 128fc2 + Alk'^ + 6k^) 



A^ 



(l + fc2)(2 + fc2)(3 + F)(3 + 2fc2) 3! 
6(16 + 13P + 2fc4)(45 + 22fc2 + 4,k^ 



) 



(1 + fc2)(2 + fc2)(3 + /!c2)(4 + fc2)(3 + 2A:2) 4! 

(M2) 



+0(A5) . 
The counter-term has the expression 

A 2 (fc2 + 2) a2 



J'\k,X)^- 



fc2 + l (fc2 + i)2 2! 

2 (3fc2 + 7) A3 12 (2fc2 + 5) A^ 



(fc2 + l)^ 3! {k^ + iy 



4! 



(M3) 
O(A^). 



As requested for a counter-term, in the sum J^{k, A) + 
J''^*(A;, A), the terms proportional to l/fc2 and 1/fc^ at large 
k cancel, and one is left with 



J{k,X)+J 

-A2 + 



•(fc,A) 
A3 5A4 



OA'' 



1 

fc6 



O 



(M4) 



2. Diagrammatic calculation of the lowest-order cumulants 



We recall from Section HITBl that 



n=l 



\n 

uf—=vZ{X}+0{v') 



(M5) 



The cumulants, or equivalently the moments, were computed 
at tree level up to n = 5 (and arbitrary times), in Section |III] 
i.e. using only the local cubic vertex. Here we compute the 
1-loop correction to this result, at equal times, and show how 
the result ( I345l l, after re-expansion in A, is recovered. The di- 
agrammatic rules are those of the simplified theory, which has 

(i) a cubic, local-in-time vertex proportional to (T = — A'(0"''); 

(ii) a non-local-in-time quartic vertex proportional to A"(0), 
which comes from the (simplified) interaction 



'^d™^ ^ ^'^ ^xtUxt{v + Uxt) 
Jxt 

+ iA"(0) / UxtUxt'{v + Uxt){v + Uxt') ■ (M6) 

^ Jxt 



Due to the quartic vertex, 1-loop diagrams are now possible 
in contrast to the cubic theory, which has only tree diagrams . 
Since we use dimensionless units below, we set cr 1 and 
A"(0) — )- —A. Note that we have written the action ( IM6I ) in 
the co-moving frame to make apparent the v terms, but the 
calculation can also be made in the laboratory frame; then one 
must remember that u has an average v. 

Let us first discuss the two lowest orders and their diagram- 
matic representation. 

To order A there is a single diagram 



Ut 



1 



1 +fc2 



(M7) 



This term involves the vertex A" (0) represented by the dashed 
lines. It is also the usual representation of the disorder vertex 
A(ii) and identifies to it whenever there are 2 entering legs. 
Since all our contribution are 0{v) the v has been chosen in 
the lowest v + ii field, which will be the case in all diagrams 
written in this section. Propagators with arrows are bare re- 
sponse functions, l/(fc2 + l)e~'^ ■* in Fourier External 
arrows are in the same number as n in m" to match the exter- 
nal ii fields. External legs are at zero momentum (since we 
compute center-of-mass velocity moments) but internal ones 
carry momentum, to be integrated over (1-loop diagrams). 
To order A2 (two outgoing lines), there are 4 contributions: 



2(3 + /e2) 



„ (l + F)(2 + fc2) 
v{2h + 2/2 + 4/3 + 4/4) 



(M8) 



Do 




h = 



2(l + fc2)(2 + fc2) 



1 1 



^2 = 7: 




2 1 



1 



(M9) 



1 



2(l + fc2)(2 + fc2) 
(MIO) 



2* 



Da 




1 1 



^4 = T 



42 + /c2 



(Mil) 



We see that both the cubic and the quartic vertices appear 
One can check that the sum of these terms with their indicated 
weights reproduces (IM2b . 
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At third order, one has 



1 2i, 3 



it^ — V 



V 



2(108 + 128fc2 + 47fc4 + 6fc6) 
'k (l + A:2)(2 + fc2)(3 + P)(3 + 2/fc2) 

16 5 8 21 



A:2 + 2 fc2 + 3 2fc2 + 3 P + i 
11 

vJ2d^I^■ (M12) 

i=l 



This comes from 1 1 diagrams: 




Ti 



6(1 + 



T. = 



6(l + /c2)(2 + fc2)(3 + fc2) 




5 + A:2 



6(l + /c2)(2 + fc2)(3 + fc2) 



T4 = 




6(l + A:2)(2 + fc2) 

U 2* 3* 




di = 12 



(Ml 3) 



do = 12 



(M14) 



= 12 



(Ml 5) 



dA = 12 



(Ml 6) 



= 24 



Tb = 



/e = 



T7 



/7 = 




7 + 4A:2 



6(1 + fc2)(2 + fc2)(3 + p)(3 + 2fc2) 

U 2* 3* 




1 



3(3 + fc2)(3 + 2fc2) 




19 + 5fc2 



36(2 + /c2)(3 + fc2) 

U 2A A3 



10 




'10 



6(2 + fc2)(3 + fc2) 

A3 



Til = 



hi = 




18(3 + /c2) 



dfi = 24 



(Ml 8) 



dr = 12 



(Ml 9) 



= 12 



(M20) 



dq = 12 



(M21) 



dio = 24 



(M22) 



dii = 24 



(M23) 



3(l + A:2)(2 + fc2)(3 + fc2) 



(Ml 7) 



This calculation illustrates how the complexity increases 
formidably with the order, and how powerful the algebraic 
method developed in section|lV]is in summing these contribu- 
tions. 
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Appendix N: Series expansion of the a„ 



Appendix P: Adiabatically switcliing on of the disorder 



The bj defined in the text can be obtained, for j > 3 as 

_ -16 + 2*= X 7 - fc[10 - 2^^ + 3fc(fc + 1)] 
^ k{k + l){k - l){k - 2) 

+6$L(-l,l,fc-2) , (Nl) 

where ^hio,, b, c) is the Lerch-$ function. 

Appendix O: Small- velocity behaviour 

Let us discuss in more detail the expansion of 5V{u) at 
small ii, looking also at subdominant terms. Denoting s := 
—A and thus k = s/(l + s), we can expand at large s, 

Li,(K) = -ln(l-K).^i^i^^ , (Ol) 

for j = 1, 2, where is analytic around n — 1 and 

0,(1) =C(i). Hence 



1 + s 



= ln(l 



ln(l + i)] 



m 



C(j) 



oo , 
SP 



We also have 
Inn 



n V 1 + s 



-(Ins)^ +lns 27) 



-K + 0[- 



(02) 



Hence we find for large s 



5Z{\ 



In^ 



-(Ins)' + (27E + 6i)ln.s + (&2 - 1) 

s 

-CI j + analytic . 



We have the following Laplace transforms: 

LTii^^, = i (In s)^ + 7e In s + analytic 
u 2 

1 



LTu^s - = - In s 
u 

LTii_j.su" Inw = — n! 



In; 



Bn 



(03) 

(04) 
(05) 
(06) 



for n = 0, 1, where in the first two lines the Laplace Trans- 
form is defined via the correctly subtracted formula. We can 
surmise that 



+K + 0{u) . 



I + ^u{b^ - 1 + Oiii)^ 



(07) 



In this appendix, we recuperate the missing terms of the 
velocity theory, as discussed in section lrVF2l It is suggestive 
from the discussion in that section, that these terms could be 
boundary terms, lost in a partial integration in time. Since the 
theory is causal, the time in question is t ^ — oo; physically 
it is related to the preparation of the system: Remind, that we 
crucially use that we are in the Middleton state. 

In order to be on the safe side, we could switch on the dis- 
order adiabatically slowly, which will suppress any boundary 
terms at time t = — oo, since there is no disorder at that time. 

Let us start from the equation of motion for the velocity (for 
short-ranged elasticity, and a source Wxt constant in space) 

(dt - + l)u,t = dt [F{vt + u^u x)gt] + m'Swt (PI) 

We have added an adiabatic factor gt which can e.g. be chosen 
as 



gt = e 



-St 



with (5-^0. 



(P2) 



Note that the exact form is not crucial, but this particular 
choice will simplify some of the ensuing calculations, since 
gt = gt-t-gt'- This gives 



S — ~Sq — Sdis 



Jxt 



(P3) 
(P4) 



-5dis = / u^tUxt'dtdf [A{v{t - t') + Uxt - Uxf )gtgt' 



_c(o) -5(2) 

dis dis dis 



(P5) 



-s. 



(0) 



UxtUxt'gtgt'dtdf A{v(t - t') + Uxt - Uxt') 

(P6) 



-'^dii = / UxtUxt'gtgt'dt'A{v(t~t') + Uxt-Uxt') (P7) 



UxtUxt'gtgt'A{v{t - t') + Uxt - Uxt') (P8) 



We now study corrections to S^is, which may intervene in our 
generating function e'^''^"''. Noting that all diagrams contain 
response-functions which decay in time at least exponentially 
fast, or more precisely faster as 



(P9) 



we have two types of diagrams for our new perturbation ex- 
pansion (for the case of interest 5 — 0): 

(i) Connected Diagrams: The disorder vertex at time t is 
attached to t = via a string of response functions; 
then we can make the replacement gt 1, and gt — > 0. 
Especially this reproduces all diagrams of the velocity 
theory. Only the vertex 5*^"^ contributes. E.g. all dia- 
grams given in appendix lM 2l are of this form. 
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(ii) Disconnected Diagrams: If the disorder at time t is not 
attached to t = via a string of response functions, then 
the integral over gt may produce a factor of j^gt — 1, 
even though gt ~ (5. As a consequence, —t is of order 
1/(5, and all response functions connected via a string of 
response functions to t may have both time-arguments 
at very large negative times, and thus are to be evaluated 
in the flat background = (since (ut) — ?■ for t — > 
—oo.) (For an example see below). 

We now discuss the leading-order correction. It comes from a 
term with one dtgt, i-e. from — iS^J] : 



UxtUxt'gt'gtdtA{v{t - t') + Uxt - Uxf) 

UxtUxt'gt'gtiv + Uxt)A' i / dT[v + Uxr] 
xtt' \Jt' 

(PIO) 



In order to conform to the rules discussed above, ut must 
somehow be connected to t = 0, whereas u-t' may not. This 
gives the only possible diagram 



S. 



(1) 



(pii) 



The times are <' < r ^ t < 0, where only t — t will become 
very large, ~ |. Therefore we can set gf = gt'-Tgr-tgt ~ 
gt-T, and the ensuing integral J^^^ gr-t = 1- The dotted 
line indicates this factor of J^^^gr-t- Furthermore, since 
both times t and t' are very negative, the response function 
Rkrt' ^ Rkrt' ■ This gives 



At leading order we now have to replace the remaining fields 
by their expectations; we also drop the term of 0{S): 




0tt2 Ut 



(P13) 



This is exactly the additional term found in Eq. (14361 ), or in 
the more rigorous derivation in Eq. (1537b . 

(2) 

We also note that S^^^ can not contribute (at least at leading 
order), since we need to gain 2 free time integrals. That im- 
plies that both response-fields must be contracted inside the 
interaction, which is impossible due to causality. However 
there will be a contribution at 2-loop order 

Further we note that, in spirit, the above derivation is 
similar to the one given in section IIVF2I In both cases, 
it was important that the second derivative of the disorder 
A"{v{t — t') + Ut — Ut'), decays, as a function of the time- 
distance t — t', to 0, which allows for a partial integration 
(eating up the time derivative gt). 



S, 



(1) 



t'f--- 



•k Jt'<T<t 

A"(0+) 



{V + Uxt)Rkrt'gr^tA"{0+)Uxt + 0{S) 

j{v + uxt)uxt + 0{5) (P12) 
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